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RESUMO 

Uma Contribuição ao Estudo de Controiadores Robustos 

O controlador preditivo baseado em modelo (MPC) tem sido aplieado eom sueesso na 
indústria, eom espeeial destaque na indústria petroquimiea. A earaeterístiea básiea dos 
algoritmos preditivos baseados em modelos é a formulação de um problema de otimização, 
para o eáleulo da seqüêneia de ações de eontrole que minimizem uma fúnção de desempenho 
num horizonte temporal eom a melhor informação disponível a eada instante, e submetida a 
restrições de operação e ao modelo da planta. A presença de ineertezas no modelo pode 
resultar em perda de desempenho da planta eontrolada baseada em eontroladores eom modelo 
nominal ou até mesmo na instabilidade da malha feehada. A elasse de eontroladores 
preditivos robustos (RMPC) baseia-se na eonsideração explíeita das ineertezas. Este trabalho 
apresenta dois algoritmos de RMPC. Um algoritmo que adieiona ao problema de eontrole 
robusto, eom ineerteza invariante no tempo na forma de otimização mín-máx, um eonjunto de 
restrições eom garantia de estabilidade robusta (RMPC-MMR). E outro algoritmo que 
desenvolve o RMPC usando desigualdades de matrizes lineares (EMI) e ineerteza na forma de 
politopo (LMI-RMPCR). Os algoritmos são avaliados, através de simulações eomputaeionais, 
aplieados em diversos sistemas típieos dos proeessos industriais, eomparando-se as respostas 
dos mesmos eom as respostas de eontroladores elássieos. Foram implementados 
experimentalmente os algoritmos de eontrole elássieo IHMPC {Infinite Horizon Model 
Predictive Control) e MPC-RS {Model Predictive Control with Reference System) em uma 
planta piloto de eontrole de nível e pH, eonstruída na Unidade Experimental da Faeuldade de 
Engenharia Químiea da Universidade Federal de Uberlândia. Os resultados de simulação dos 
eontroladores robustos foram avaliados eom as repostas de simulação e experimentais dos 
eontroladores elássieos implementados na planta. Os resultados indieam a aplieabilidade e as 
limitações dos eontroladores propostos. 


Palavras-Chave: Controle de Proeessos, Ineerteza, EMI, MPC, RMPC. 





ABSTRACT 


A Contribution for Studing of Robust Control 

The model based predictive controller (MPC) has been successfully applied in industry, with 
particular emphasis on the petrochemical industry. The basic feature of MPC algorithms is 
based on the formulation of an optimization problem, for calculating the sequence of control 
moves that minimize a performance function in a time horizon with the best information 
available at every moment, and submitted to restrictions, operation and the plant model. The 
presence of uncertainties in the model can result in loss of performance of the model based 
control System or even instability in the closed loop. The class of robust predictive controllers 
(RMPC) is based on the explicit consideration of uncertainties. This research presents two 
RMPC algorithms. One of them adds a set of restrictions with robust stability guarantee 
(RMPC-MMR) to the problem of robust control, with uncertainty in time variant as min-max 
optimization. And the other develops RMPC using Linear Matrix Inequalities (LMI) and 
uncertainty as a polytope (LMI-RMPCR). The algorithms are evaluated through Computer 
simulations, and are applied to a range of industrial process typical Systems. The responses are 
compared with the responses of traditional controllers. The classical algorithms IHMPC 
(Infinite Horizon Model Predictive Control) and MPC-RS (Model Predictive Control with 
Reference System) are also assessed experimentally by the application in a pilot plant for 
controlling levei and pH, built in the Experimental Unit of the Chemical Engineering College, 
at Federal University of Uberlândia. The results of RMPC simulation were evaluated with the 
simulation and experimental responses of the classical controllers implemented in the plant. 
The results indicate the applicability and limitations of the proposed controllers. 

Keywords: EMI, MPC, Process Control, RMPC, Uncertainty. 
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1 INTRODUÇÃO 


O controle preditivo baseado em modelos (MPC) tornou-se uma importante 
estratégia de eontrole para aplieações industriais, devido prineipalmente ao sueesso de suas 
implementações. Os diversos algoritmos de eontrole preditivo diferem pela forma do modelo 
usado para representar o proeesso, o ruido e pela função eusto a ser minimizada. Existem 
muitas aplieações bem sueedidas de eontrole preditivo, não apenas na indústria químiea eomo 
em outras áreas (QfN; BADGWELL, 2003). 

O eoneeito de eontrole preditivo robusto (RMPC) está assoeiado eom a utilização 
explíeita de um modelo para as ineertezas na formulação do eáleulo da lei de eontrole. Os 
algoritmos de eontrole preditivo robusto diferem pela forma de representação das ineertezas e 
pela restrição de estabilidade, além das diferenças eitadas para eontroladores MPC 
eonveneionais (MAYNE et al, 2000). 

1.1 Justificativas 


O eontrole de proeessos eom restrição é um dos mais desafiadores problemas na 
teoria de eontrole. A exigêneia se toma ainda maior quando se trata do projeto de eontrole 
robusto de sistemas. O problema de restrição robusta não é apenas signifieativo, mas um 
instigante problema na teoria de eontrole. 
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Durante as décadas passadas, uma sólida fundamentação teórica de MPC 
possibilitou a aplicação prática dessa estratégia com estabilidade garantida e certo grau de 
robustez. Pode-se afirmar que hoje existe um entendimento bastante apurado sobre o projeto 
de MPC linear robusto, mas ainda há espaço para pesquisas em versões desse controlador 
aplicáveis a problemas não lineares, que exploram eficientemente o conhecimento de 
incertezas. Como obstáculos aos algoritmos de MPC robusto, pode-se citar a alta demanda 
computacional e a pequena aplicabilidade decorrente da grande exigência computacional 
associada a esses algoritmos. 


1.2 Delimitações e Objetivos deste Trabalho 


A meta desse estudo foi propor algoritmos de controle preditivo robusto por 
construção que mantenham a estabilidade robusta e o desempenho do sistema em patamares 
aceitáveis. 

Dois algoritmos foram propostos: um que utiliza incerteza na entrada invariante 
no tempo com otimização mín-máx e restrição de robustez (RMPC-MMR), e outro que utiliza 
incerteza na forma de politopo com o problema de otimização na forma de LMl {Linear 
Matrix Inequalities) e restrição de robustez (LMl-RMPCR). Os algoritmos foram avaliados 
pela simulação de diversos modelos representativos de processos industriais, que exploram de 
alguma forma características de estabilidade, robustez ou desempenho dos controladores. 
Foram implementados controladores clássicos (IHMPC, RMPC-MMTV, RMPC-MM, MPC- 
RS, LMl-RMPC, MPC-ML), cujas formulações serão apresentadas em detalhe no Capítulo 2, 
a fim de se obter uma base comparativa. 
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Foram realizadas também as implementações experimentais dos eontroladores 
MFC: IHMPC e MPC-RS na planta piloto de eontrole de nivel e pH, eonstruida na Unidade 
Experimental da Faeuldade de Engenharia Quimiea da Universidade Federal de Uberlândia 
(UFU). O eontrole MPC-RS foi eseolhido para esse trabalho, eonsiderando-se que o Núeleo 
de Controle e Otimização de Proeessos da UFU já possuía experiêneia nessa elasse de 
eontroladores na unidade experimental. Os eontroladores robustos desenvolvidos são da elasse 
do algoritmo de eontrole de horizonte infinito, exemplifieado pelo IHMPC. Desta forma, os 
resultados de simulação entre os algoritmos desenvolvidos e o IHMPC forneeem uma boa 
eomparação para fins de aplieabilidade. Os resultados de simulação e experimentais dos 
eontroladores implementados na planta foram eomparados eom os resultados simulados dos 
eontroladores propostos. As eonsiderações deste trabalho indieam a aplieabilidade e 
limitações dos eontroladores robustos. 

Este trabalho buseou o treinamento em pesquisa, o desenvolvimento de algoritmos 
de eontrole MPC e RMPC e também a redução da distâneia entre a teoria e a prátiea 
industrial. Como eontribuições, além dos desenvolvimentos teórieos dos eontroladores e sua 
implementação real, pode-se eitar a divulgação e diseussão dos sistemas de eontrole pela 
publieação e partieipação em diversos eongressos e fóruns de diseussão téeniea e eientífiea. 


1.3 Resumo do Trabalho 


O Capítulo 2 apresenta a revisão bibliográfiea aeerea dos seguintes assuntos: 
MPC, ineerteza e robustez, RMPC e LMI, e ainda apresenta os algoritmos elássieos utilizados 
para o desenvolvimento de novas estruturas e eomo base eomparativa dos resultados. O 
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Capítulo 3 trata do desenvolvimento de dois algoritmos de eontrole robusto baseados na 
adição de restrição de robustez ao problema de otimização mín-máx e ineerteza na entrada 
invariante no tempo (RMPC-MMR) e ineertezas na forma de politopos, eom restrição de 
robustez e problema de otimização na forma de LMI (LMI-RMPCR). O Capítulo 4 mostra a 
simulação em diversos proeessos industriais. No Capítulo 5 é deserito o sistema experimental. 
No Capítulo 6 são apresentados os resultados de malha aberta da parte experimental. No 
Capítulo 7 são mostrados os resultados de simulação e experimentais de malha feehada no 
sistema experimental. No Capítulo 8 apresentam-se as eonelusões juntamente eom as 
sugestões de trabalhos futuros. Eneerrando o trabalho, eneontram-se as referêneias usadas na 
elaboração deste trabalho. No Apêndiee é possível ver as eurvas de ealibração, que fazem a 
relação entre as medidas da unidade experimental, eomo por exemplo, vazão das soluções, 
eom a tensão dos transmissores. Pode-se também visualizar as telas dos programas em 
Labview® que foram implementadas para a operação da unidade experimental em malha 
aberta e em malha feehada. 


1.4 Publicações 


Durante a elaboração deste trabalho de doutorado foram realizadas duas 
publieações em periódieos, treze publieações em anais de eongressos, sendo duas delas de 
natureza profissional e quatro resultantes de publieações do trabalho de mestrado em eontrole 
preditivo. Houve também a partieipação em vinte e eineo eventos de natureza eientífiea. 
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Os artigos listados abaixo tratam do desenvolvimento dos algoritmos de eontrole 
preditivo robusto e sua aplieação em diversos proeessos, eomo por exemplo, os modelos de 
eolunas de flotação, máquinas síneronas e modelos de pequenas eentrais hidrelétrieas e outros 
sistemas lineares. Os resultados desses estudos não se eneontram deseritos nesta tese, mas os 
mesmos eneontram-se alinhados eom as eonelusões desse trabalho. 


Artigos publieados em periódieos: 


SILVA, C. H. F., HENRIQUE, H. M., GOMES, L. G, OLIVEIRA-LOPES, L. C. Controle 
Preditivo Robusto Baseado em Inequações Matrieiais Lineares Aplieado a Máquinas 
Sineronas. Revista Ciências Exatas, v.I3, n° I, p.I2I-I29, 2007. artigo seleeionado para 
publieação na revista dentro dos anais do 1° Congresso Internaeional de Cooperação Indústria- 
Universidade 2005. 


SILVA, C. H. F., HENRIQUE, H. M., OLIVEIRA-LOPES, L. C., GOMES, L. R. Controle 
Preditivo Robusto e sua Aplieação no Modelo de Máquinas Síneronas de Pequenas Centrais 
Hidrelétrieas. Revista Controle & Instrumentação, v.ll4, p.72-77, 2006. Artigo 
selecionado dentro dos anais do IV Seminário Nacional de Controle e Automação, realizado 
em 2006. 


Artigos publicados em anais de congresso: 


SILVA, C. H. F., GOMES, L. G, HENRIQUE, H. M., OLIVEIRA-LOPES, L. C. Controle 
Preditivo Robusto: Um Estudo Comparativo. XVI Congresso Brasileiro de Engenharia 
Química, Santos, 2006. 


SILVA, C. H. F., GOMES, L. G, HENRIQUE, H. M. Controle Preditivo Aplicado À Máquina 
Síncrona. IV Conferência de Estudos em Engenharia Elétrica, Uberlândia, 2005. 


SILVA, C. H. F., GOMES, L. G, HENRIQUE, H. M., OLIVEIRA-LOPES, L. C. Controle 
Preditivo Robusto Baseado em Inequações Matriciais Lineares Aplicado a Máquinas 
Síncronas. I Congresso Internacional de Cooperação Universidade-Indústria, Ubatuba, 
2005. 
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SILVA, C. H. F., OLIVEIRA-LOPES, L. C., HENRIQUE, H. M., GOMES, E. G Controle 
Preditivo Robusto de Pequenas Centrais Hidrelétrieas. IV Seminário Nacional de Controle e 
Automação, Salvador, 2005. 


SIEVA, C. H. F., OEIVEIRA-LOPES, E. C., HENRIQUE, H. M. CONSTRAINED ROBUST 
Model Predictive Control for Linear Systems. XV Congresso Brasileiro de Engenharia 
Química, Curitiba, 2004. 


SILVA, C. H. F., OLIVEIRA-LOPES, E. C., HENRIQUE, H. M., GOMES, E. G Controle 
Preditivo Robusto Aplieado a Pequenas Centrais Hidrelétrieas. IV Simpósio Nacional de 
Pequenas e Médias Centrais Hidrelétricas, Reeife, 2004. 


SIEVA, C. H. F., HENRIQUE, H. M., OLIVEIRA-LOPES, E. C. Controle Preditivo Robusto 
de um Sistema de Flotação. XXX Encontro Nacional de Meios Porosos, Uberlândia, 2004. 
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2 REVISÃO BIBLIOGRÁFICA 


2.1 Controle Preditivo Baseado em Modelos - MPC 


A teoria de eontrole preditivo baseado em modelos teve origem no final da déeada 
de 1970 e desenvolveu-se de forma eonsiderável desde então (CAMACHO; BORDONS, 
1999; MAYNE et ai, 2000; NIKOLAOU, 2001; MORARI, 1993; SAFFER II e DOYFE III, 
2004, RAWFINGS, 1999). 

Os primeiros trabalhos foram Model Predictive Heuristic Control (MPHC), 
posteriormente eonheeido eomo Model Algorithmic Control - MAC), formulado por Riehalet 
et al. (1978), e Dynamic Matrix Control (DMC), formulado por Cutler e Ramaker em 1980, 
levando o eontrole preditivo a se tornar bastante eonheeido na indústria petroquímiea 
(CAMACHO; BORDONS, 1999). Esses eontroladores estão relaeionados eom o problema de 
eontrole ótimo de tempo mínimo e programação linear, inieialmente apresentada por Zadeh e 
Whalen em 1962 (GARCIA; PRETT; MORARI, 1989). 

O eontrole MPC também pode ser formulado no eontexto do eontrole moderno. O 
uso de modelos na forma de espaço de estados, em detrimento do uso de outros modelos, 
durante a déeada de 1990, foi responsável por um amadureeimento substaneial na teoria de 
eontrole preditivo (EÔFBERG, 2003). O uso dessa formulação permite não apenas a aplieação 
de teoremas eonheeidos da teoria de espaço de estado e sistemas lineares, eomo também 
faeilita generalizações para easos mais eomplexos. Essa perspeetiva leva a um simples ajuste 
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de regras para estabilidade e robustez (RICKER, 1990). O eontrolador MFC, nesse easo, pode 
ser interpretado eorno sendo um eompensador baseado em observador de estado, e sua 
estabilidade, desempenho e robustez são determinados pelos pólos do observador, fixados 
diretamente por ajuste de parâmetros, e pelos pólos do regulador, determinados pelos 
horizontes e pelos pesos de desempenho. 

O desempenho dinâmieo de um sistema eontrolado é grandemente determinado 
pela loealização dos pólos em malha feehada. Para se obter um bom desempenho é neeessário 
aloeá-los em uma região adequada. Para sistemas lineares sem restrição, estudos do problema 
de eoloeação de pólos ajustam a região eireular para aloeação dos pólos usando realimentação 
de estado e também das saídas. Para o problema eom restrição, entretanto, a eoloeação de 
pólos situa-se dentro de uma pequena região eireular, aearretando a redução substaneial da 
viabilidade da região inieial (ZHANG E XI, 2008). Wang e Rawlings (2004a e 2004b) 
apresentam um exeelente tutorial para a teoria de MPC indieando problemas e direções para 
essa importante téeniea de pesquisa aeadêmiea e de aplieação industrial. 

A primeira geração de MPC, eom algoritmos eomo DMC e IDCOM 
{Identification and Command), teve grande impaeto no eontrole de proeessos industriais e 
serviu para definir o paradigma inieial do MPC industrial. 

O algoritmo QDMC {Quadractic Dynamic Matrix Control) desenvolvido por 
Gareia e Morshedi (1986) pode ser eonsiderado representante da segunda geração de 
teenologia MPC, earaeterizando-se por desenvolver sistematieamente uma forma de 
programar restrições de entrada e saída. Essa formulação aperfeiçoa o MPC na medida em que 
posieiona o problema de eontrole eomo um problema de otimização quadrátiea (QP), euja 
solução pode ser obtida, efieientemente, por eódigos QP disponíveis. O ajuste do eontrolador 
para operação na presença de restrições pode ser muito difieil, mesmo quando o modelo do 
proeesso é perfeito. Muito esforço tem sido direeionado para a simulação da malha feehada a 
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fim de comissionar o controlador. Deve-se reconhecer, entretanto, que a simulação de muitas 
combinações possíveis da ação da restrição é de difícil execução e possui reflexos na 
aplicação prática. 

Desta forma, parece sensata a utilização de algoritmos que sejam construídos de 
forma inerentemente estável para todos os casos. Na teoria de controle tais estruturas são 
conhecidas como algoritmos nominalmente estáveis (QIN; BADGWELL, 2003). Como 
exemplo de formulação de estabilidade garantida, surgiu o algoritmo de MFC com horizonte 
infinito (IHMPC - Infmite Horizon Model Predictive Control), apresentado por Muske e 
Rawlings (1993a), que foi adotado pelos acadêmicos como o paradigma dos anos de 1990, 
constituindo a terceira geração de MFC e ainda considerando a estabilidade como conceito 
central nos estudos de sistemas dinâmicos (NIKOLAOU, 2001; MORARI; LEE, 1999). Essa 
geração se distingue devido aos vários níveis de restrições desenvolvidas por vários 
mecanismos. A solução viável é procurada e direcionada para variações na planta em tempo 
real, fornecendo um rico conjunto de soluções para realimentação. Permite-se, deste modo, 
uma larga faixa de processos dinâmicos (estáveis, integradores, instáveis) e especificações de 
controle. 

A quarta geração de MFC contempla preocupações como a consideração da 
incerteza e o controle robusto; otimização em múltiplos níveis direcionada para a priorização 
dos objetivos de controle; melhoria na tecnologia de identificação, dentre outras (QIN; 
BADGWELL, 1996 e QIN; BADGWELL, 2003). 

MFC não é um novo método de projeto de controle. Ele essencialmente resolve 
problemas de controle ótimo padrão. O que o diferencia é o fato de que, no MFC, o problema 
de controle ótimo requerido possui horizonte finito, em contraste com o horizonte infinito, 
usualmente empregado nos controladores ótimos lineares H 2 e //«. 
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Figura 2.1 - Genealogia aproximada para algoritmos MFC linear - adaptado de Qin e 
Badgwell (2003). 

A Figura 2.1 mostra de forma aproximada a genealogia dos algoritmos MFC 
industriais. A Tabela 2.1 apresenta uma síntese da teenologia industrial sobre MFC (QIN; 
BADGWELL, 2003; CAMACHO; BORDONS, 1999). 



TABELA 2.1: MFC LINEAR INDUSTRIAL - TECNOLOGIA DISFONÍVEL. 


Companhia 

Nome do Produto 

Descrição 

Adersa 

IDCOM 

Identificação e Comando 


HIECON 

Controle de Restrição Hierárquico 


FFC 

Controle Freditivo Funcional 


GUDE 

Facote de Identificação 

Aspen Teeh 

DMC-plus 

Facote de Controle por Matriz Dinâmica 


DMC-plus model 

Facote de Identificação 

Honeywell 

FCT 

Tecnologia de Controle Freditivo 


RMFCT 

Controle Freditivo Baseado em Modelos Robusto 

Shell Global 

SMOC-II 

Controle Shell Multivariável Otimizante 

Invensys 

Connoisseur 

Facote de Identificação e Controle 


Fonte: Qin e Badgwell (2003) 
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O MPC resolve o problema de eontrole ótimo em malha aberta em tempo real, 
eonsiderando o estado eorrente do sistema a ser eontrolado eomo mais que uma politiea de 
realimentação off-line. É neeessário que esse problema de otimização seja resolvido dentro do 
tempo de amostragem (MAYNE et al, 2000). Bemporad (2006) apresenta uma sintese sobre 
os eontroladores preditivos, eom tendêneias e ferramentas eomputaeionais. 

Na teenologia de MPC, áreas eomo estabilidade e grau de ótimo têm atingido 
sufieiente maturidade. A razão de ser do eontrole MPC é sua habilidade para lidar eom 
problemas de eontrole em que a eomputação da lei de eontrole off-line seja difíeil, embora 
outras oaraeterístioas, eomo sua eapaeidade para eontrolar plantas multivariáveis, tenham sido 
julgadas inieialmente eomo oaraeterístioas mais importantes (MAYNE et al, 2000). Cada 
aplioação impõe suas restrições. Os atuadores, por exemplo, possuem limites para aplioação 
de força. A habilidade em manusear as restrições de entrada, de saída e de estado interno pode 
ser atribuída eomo a mais signifioativa oontribuição para o suoesso na aplioação industrial 
(HENSON; SEBORG, 1997; DE OLP/EIRA; MORARI, 1994; LEE, 1997). 

Na ausênoia de erros de modelagem e perturbações, o eontrole MPC usando 
matrizes de pesos adequados e horizonte de predição sufioientemente longo pode ser visto 
eomo uma implementação exata da politiea de eontrole ótimo em realimentação (L0VAAS et 
al, Inpress). 

Outras vantagens a serem oonsideradas são: os limites de segurança; a 
possibilidade de obtenção de vantagens no emprego dessa téeniea em plantas não lineares sem 
restrição, bem eomo em plantas variantes no tempo; a eonsideração dos ruídos e dos efeitos 
reais do proeesso na formulação de restrições a serem atendidas quando da operação da 
planta; o melhor desempenho do eontrole; redução dos requerimentos de manutenção, e 
melhoria da flexibilidade e da agilidade (EATON; RAWLINGS, 1992). 
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Considerando a aplicação em processos MIMO {Multiple Input Multiple Output), 
0 MPC lida diretamente com a coordenação e balanceamento das interações entre as entradas 
e as saídas e também com as restrições associadas (ZHANG e XI, 2008). Mas apesar de ser 
realmente efetivo em situações apropriadas, ele ainda possui limitações que geralmente não 
são mencionadas pelos vendedores, como por exemplo, dificuldades de operação, alto custo 
de manutenção (uma vez que exige mão-de-obra especializada) e perda de flexibilidade, que 
podem resultar na fragilidade do controlador. Essas limitações não estão ligadas apenas ao 
algoritmo propriamente dito, mas também à justificativa desses controladores, à aplicabilidade 
e à manutenção (HUGO, 1996). 

Versões industriais de MPC não são estabilizantes automaticamente, necessitando 
de escolhas e considerações de forma a se obter propriedades de estabilidade. O sucesso 
industrial do MPC levou à investigação teórica da estabilidade. A ferramenta de análise de 
estabilidade por excelência é a teoria de estabilidade de Lyapunov. A estabilidade teve que ser 
direcionada restritivamente para sistemas lineares, configurados como MPC sem restrição. Os 
resultados estabeleceram a existência de controle de horizonte finito de tal forma que o MPC 
seja estabilizante. Hoje, o valor da função custo de controle é usado amplamente como uma 
função natural de Lyapunov para análise de estabilidade de MPC (MAYNE et al, 2000). 

A revolução de controle da década de 1960 foi grandemente influenciada pela 
teoria de Hamilton-Jacobi-Bellman ou programação dinâmica, cujo teor providencia 
condições suficientes para o grau de ótimo, um procedimento construtivo para se determinar 
um controlador ótimo em realimentação na forma u = Kx, e motiva a computação de 
algoritmos para a determinação de controle ótimo em malha aberta para um dado estado 
inicial, u'’(k; x), onde k é o instante de amostragem. O grau de ótimo não implica em 
estabilidade, mas sob certas condições, controladores ótimos de horizonte infinito são 


estabilizantes. 



15 


Nesse eontexto, é eonveniente eitar o MPC sub-ótimo. Essa elasse é eonstituída 
por eontroladores eujo grau de ótimo não é exigido, e sim a viabilidade da solução. Tal 
proposta não defende o sub-ótimo, mas permite que resultados satisfatórios de eontrole sejam 
eneontrados quando o ponto ótimo se torna impratieável. 

Existe ainda uma distâneia entre os resultados observados e a teoria (HENSON; 
SEBORG, 1997). Resultados experimentais, em termos de aplieação e simulação, eontinuam 
indieando um bom desempenho do MPC em uma larga elasse de problemas, mas a sintonia do 
algoritmo para obter estabilidade ainda é usualmente neeessária. Resultados teórieos requerem 
horizontes infinitos ou restrição terminal para assegurar a estabilidade, o que representa 
eondições que podem ser fortes demais para a implementação prátiea. Como alternativas para 
essas aproximações, é possível garantir a estabilidade eom o uso de uma função de penalidade 
terminal, eonforme é deserito a seguir. 

A viabilidade do problema de otimização a eada instante de tempo k deve ser 
garantida. Tipieamente, eonsidera-se a viabilidade no tempo k = 0 e se eseolhe a função eusto 
e a restrição de estabilidade de tal forma que a viabilidade é preservada nos passos dos 
instantes seguintes (BEMPORAD; MORARI, 1999). 

As téenieas populares para “forçar” a estabilidade podem ser divididas em duas 
elasses prineipais: 

• A primeira elasse usa os valores da função eusto dado por V(k) = j(s,Xj^,Np,Nj^), para 

atingir por minimização a seqüêneia de ações de eontrole ealeulada em 
H = i^j, em eada instante k, eomo uma função de Eyapunov. Np é o 


horizonte de predição e N^ é o horizonte de eontrole. 
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• A segunda classe requer explicitamente que o estado x(t) seja encolhido segundo alguma 
norma. É de grande utilidade ter um simples número que seja capaz de ter informações 
sobre a dimensão global de um vetor, uma matriz, um sinal ou um sistema (SKOGESTAD 
e POSTLETHWAITE, 1997). Para atender a esse propósito são usadas as funções 





chamadas de norma. A norma mais comum é a norma Euclidiana 

llHl2 = Jsk 
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1 V J 
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Outras normas bastante utilizadas são a norma-1 




e a norma-co 
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As principais técnicas são sucintamente descritas a seguir 


(BEMPORAD; MORARI, 1999; MAYNE et al, 2000; BADGWELL, 1997a; RALHAN; 

BADGWELL, 2000): 

• Restrição terminal: a restrição de estabilidade é dada pela Eq. 2.1. Essa forma gera 
uma seqüência de ações de controle na forma da Eq. 2.2. E interessante notar que 
no cálculo da ação de controle no instante k, após o horizonte de controle (N^), a 
ação de controle deve ser igual a zero, neste caso considerando-se o problema 
regulador. Para o problema servo haverá um deslocamento da origem do sistema. A 
Eq. 2.3 mostra a forma da seqüência de controle sem esta restrição. Esta seqüência 
é viável no instante k + 1. Por conseqüência, a Eq. 2.3 é a função de Lyapunov do 
sistema, uma vez que se trata de uma função não-crescente. 


Xk+N^ik=0 (2-1) 

'^1 ~ l*^k+l|k’-"’^k+N^|k(2-2) 

V(k + l)<j(s„x„N,,Nj<j(E,x„N,,Nj=V(k) 


( 23 ) 
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O principal empecilho para o uso de restrição terminal é que o esforço de eontrole 
requerido para levar o estado para a origem pode ser grande, espeeialmente eom o 
horizonte de predição reduzido. Desta forma, a viabilidade se toma erítioa. Além 
do mais, o desempenho pode ser afetado negativamente devido à restrição terminal 
artifieial. O valor da função eusto assoeiado ao problema de eontrole ótimo finito 
se aproxima do problema de eontrole ótimo infinito quando o horizonte se 
aproxima do infinito. O interesse da aplieação da restrição terminal de igualdade 
em MFC reside também no fato de que ela pode, ainda, ser empregada para 
estabilizar sistemas que não podem ser estabilizados por eontroladores em 
realimentação eontinuos (MAYNE et al, 2000). 

• Conjunto Terminal Invariante: a idéia dessa téeniea é relaxar a restrição terminal 
em um eonjunto de restrições, eomo apresenta a Eq. 2.4. Novamente, a prova de 
estabilidade é baseada no argumento de Lyapunov. A proposta do eontrolador 
preditivo é levar o estado para em um tempo finito. Dentro de um eontrolador 
estabilizante é empregado. Essa forma de eontrolador é eonheeida eomo modo 
dual. 




(2.4) 


• Horizonte de Predição de Saída Infinito: para sistema assintotieamente estável, 
nenhuma restrição de estabilidade é neeessária se o horizonte de predição for 
infinito (MEÍSKE; RAWLINGS, 1993a). A prova é novamente baseada de maneira 
similar ao argumento de Lyapunov. Marjanovie e Lennox (2004) investigam a 
formulação de MFC baseada em horizonte infinito e mostram que o eonjunto 
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terminal de eondições de inelusão que é assoeiado eom a lei de eontrole induz às 
limitações no sistema de eontrole. 

• Matriz de Pesos Terminal: a estabilidade pode ser garantida sem a adição de 
restrição de estabilidade pela eseolha de uma matriz de pesos terminal (QN), 
eonforme a Eq. 2.5, eomo solução da equação de Rieatti, dada pela Eq. 2.6. Nesse 
easo, é mostrada a forma para Np = N^. 


mmJ(S,x,,Np,N^) = xí:,^QNx^^ +Z.=Ô~ ^LikQ^kuik + «L|kR« 
A^QNA - QN + A^QNb(b^QNB + r) ' B^QNA + Q = 0 


k+i|k 


(2.5) 

( 2 . 6 ) 


• Restrição de eontração: mais que eonfiar no eusto ótimo V(k) eomo uma função de 
Lyapunov, a idéia é requerer explieitamente que o estado reduza segundo alguma 
norma, eomo apresentado na Eq 2.7: 


^k+l|k 


< a Xj, ,a < 1 


(2.7) 


• Função Custo Terminal: modifiea-se a função eusto para garantir a estabilidade da 
malha feehada pela adição de um eusto terminal, 'E(x), que é não-trivial. Para o 
easo sem restrição terminal, Xj pertenee ao eonjunto de números reais de ordem n. 
O resultado dessa modifieação é equivalente a ter horizonte de predição infinito e 
horizonte de eontrole finito, podendo-se demonstrar que a estabilidade é garantida 
se 0 horizonte de predição for muito maior que o horizonte de eontrole. 
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• Custo Terminal e Conjunto de Restrição Terminal: os controladores preditivos 
mais recentes pertencem a essa categoria. Existem boas razões para incorporar o 
custo terminal e um conjunto de restrição terminal no problema de controle. De 
maneira ideal, o custo terminal ^(x) deve ser o valor da função de horizonte 
infinito Joc**. Se esse for o caso, então o valor da função objetivo no horizonte de 
controle (N^) é igual à função objetivo de horizonte infinito como se segue: JNm*^= 
Joo*'. A otimização em tempo real seria desnecessária e com conhecidas vantagens 
do horizonte infinito, tal como a estabilidade e a robustez seriam automaticamente 
atingidas. O conjunto de restrição terminal Xj é escolhido como pertencente ao 
conjunto de máximas saídas admissíveis. 

A complexidade do problema de otimização depende da escolha do índice de 
desempenho e da restrição de estabilidade. Quando o horizonte de controle é infinito ou a 
restrição de estabilidade apresenta uma forma adequada, o problema de otimização pode ser 
formulado como uma otimização quadrática. 

Essa formulação é considerada a forma padrão e sua popularidade se deve á 
conveniência matemática. Quando a norma-2 é usada, os algoritmos de programação que 
utilizam cone de segunda ordem podem ser adotados convenientemente. Alternativamente, 
pode-se obter formulações com otimização linear (LP) por utilização de norma-1 ou norma-oo. 
Um dos principais fatores que motivam a formulação de problemas MPC na forma de LP é a 
percepção de que essas formulações podem ser resolvidas mais eficientemente. E certamente 
verdadeiro que códigos modernos podem ser eficientemente manuseados para uma larga 
classe de problemas, entretanto, também é verdade que ambos, LP e QP convexo, podem ser 
moldados como problemas de complementaridade linear e resolvidos usando o método dos 
pontos interiores pela formulação do problema na forma de LMI (VANDENBERGUE et al. 


2002 ). 
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Teorema 2,1 Suponha que as seguintes eonsiderações se mantenham para um sistema linear, 
um eontrolador nominal K(x), um domínio de estado terminal Xj e um peso de estado 
terminal 'P (x), expressos a seguir pelas Eq. 2.8-2.13. 


Xk+i = Axk + Buk 

(2.8) 

0 e XtÇ X 

(2.9) 

Ax + BK(x) e Xj, V x e Xj 

(2.10) 

T^CO) = 0, T' (x) > 0, V X 0 

(2.11) 

T(Ax + BK(x)) - 'F(x) < -x^Qx - K^(x)RK(x), V x e Xj 

(2.12) 

K(x) e U, V X e Xt 

(2.13) 


Então, admitindo a viabilidade de estado inieial, um eontrolador MFC usando o seguinte 
problema de otimização apresentado pelas Eq. 2.14-2.15 garantirá estabilidade assintótiea. 

mmX£B'~^^k+i|kQ^k+i|k + »í+i|k^“k+i|k + 'í'(Xk+Nmík) (2-14) 

Sujeito à, 

«k+i|k ^ U, x^+iik e X, Xk+Njnlk ^ (2.15) 

Prova: (MAYNE et al, 2000). 

Apesar da teoria do MFC apresentar eonsiderável maturidade, profissionais 
ligados à prátiea industrial argumentam que os limites de desempenho e aplieabilidade do 
MFC na prátiea não se devem a defieiêneias no algoritmo, mas a questões ligadas a 
difieuldades de modelagem, adequação de sensores e robustez insufieiente para falhas 


(MORARI; LEE, 1999; QIN e BADGWELL, 2003). 
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Pesquisas em robustez têm levado ao entendimento do problema e às limitações 
inerentes ao MPC eonveneional (em malha aberta). Enquanto o problema tem sido estudado e 
bem entendido, os resultados de pesquisas são eoneeituais e funeionam, a prineípio, mas são 
eomplexos demais para serem empregados, prineipalmente quando se tratam de eontroladores 
MPC mín-máx em malha feehada. No momento existem várias formulações de MPC robusto. 
Enquanto os esquemas existentes eontribuem para um aumento no entendimento eles são, 
eom raras exeeções, eomputaeionalmente intratáveis na aplieação prátiea. Desta forma, 
pereebe-se elaramente a neeessidade de desenvolvimento de estratégias de MPC robusto 
viáveis sob o ponto de vista da implementação e que reduzam a distaneia entre a teoria e a 
prátiea (FfNDEISEN; ALLGÔWER, 2001). Avanços em robustez, estimação de estado e 
MPC adaptativo podem eontribuir para a evolução de MPC eom realimentação das saídas 
(MAYNE et al, 2000). O mais importante dessa diseussão é a possibilidade de 
implementação em tempo real, o que é um fato esseneial para o MPC. 

A Figura 2.2 mostra as áreas de interfaee eom o eontrole MPC eujos 
desenvolvimentos afetam na melhoria e evolução dos eontroladores preditivos (RICHALET, 
1993; RICHARDS, 2005). Além dessas referêneias bibliográfieas, as idéias transmitidas pelo 
Prof. Dr. Karl Johan Âstrõm {Lund Institute of Technology - Suéeia) durante palestra proferida 
em 28 de novembro de 2007, por oeasião do aniversário de 80 anos da Universidade Federal 
de Mina Gerais, em Belo Horizonte, eontribuíram para a eonsolidação dos eoneeitos da atual 
situação do eontrole de proeessos e para a formação das idéias de futuro indieadas na Figura 


2 . 2 . 
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Figura 2.2 - MFC e áreas de interfaee. 
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2.1.1 Formulação Clássica - IHMPC (Infiniíe Horizon Model Predictive 
Control) 

De maneira geral, os diferentes algoritmos de MPC são formulados no instante 
presente k, eonsiderando um horizonte de predição Np e um horizonte de eontrole N^, eomo 
ilustra a Figura 2.3. Usando a resposta do modelo do proeesso, as variáveis manipuladas são 
preditas. A ação da variável manipulada é seleeionada de tal forma que a resposta predita 
possui eertas oaraeteristieas desejadas de estabilidade e desempenho. Apenas a primeira ação é 
implementada e, no instante k+1, é feita uma atualização das medidas da planta, além de um 
novo eáleulo para as variáveis manipuladas e saídas no horizonte de predição. Após o 
horizonte de eontrole, eonsidera-se variação nula para a ação de eontrole. Para os estudos 
deste trabalho é feita a eonsideração de que o horizonte de eontrole é igual ao horizonte de 
predição, ou seja, Np = N^. Os desenvolvimentos em eontrole robusto serão feitos em 
algoritmos da família do MPC de horizonte infinito. 



k 

Horizonte de controle 

◄-► 

^_ Horizonte de predição 

Figura 2.3 - Horizonte de ações, eentro da teoria de MPC. 
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Um algoritmo básico é descrito a seguir: 

1. Obter 0 novo estado x por medição e/ou estimação; 

2. Resolver o problema de otimização^ J(k) 

3. Apliear apenas o primeiro movimento de eontrole Uk+o|k; 

4. k k+1, retoma para 1. 

Considerando a sua importâneia, por se tratar de um algoritmo representativo de 
estabilidade nominal garantida, bem eomo pela simplieidade (MORARI, 1993), o algoritmo 
apresentado por Muske e Rawlings (1993 a) será deserito a seguir, eonstituindo o ponto de 
partida para o estudo da estabilidade robusta. Será apresentada apenas a formulação para 
sistemas estáveis. Para sistemas instáveis, sugere-se a utilização da versão apresentada em 
Muske e Rawlings (1993b). 

O MPC eom horizonte de eontrole infinito (IHMPC - Infmite Horizon Model 
Predictive Control) é definido de forma que o movimento de eontrole Uk, no instante de 
amostragem k, é dado pelo primeiro elemento Uk|k de S:{Uk|k, Uk+i|kv • •, Uk+Nm-i|k} na forma da 
Eq. 2.16. Esta equação pode ser manipulada algebrieamente resultando na Eq. 2.17, que 
transforma o horizonte de eontrole para horizonte finito. A seqüêneia de ações de eontrole é a 
solução do problema de otimização dado pelas Eq. 2.17-2.20, eom: Q = > 0, R = > 0, 

S = > 0, Q = C^QC -I- A^QA (Equação de Lyapunov), e Auk+j= Uk+j- Uk+j-i- No instante k + 

Nm,e posteriormente, o vetor de entradas u é um vetor de zeros e z representa o vetor de 
estados reais medidos ou estimados das saídas medidas. 

jk=minx;:oyLQyk+i +«k+.R«k+i +^»LsAUk+, 

u 


( 2 . 16 ) 
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Jj, = minS^HS + 2S^(Gz]j -Fu]j_l) 


Os termos do índice dado pela Eq. 2.17 são mostrados a seguir: 

B^QB + R + 2S B^A^QB-S B^A^ ” QB 

„ B^QAB-S B^QB + R + 2S B^A^"^” QB 

Jd = 

BTqan». B B^QA*"” "B B^QB + R + 2S 

B^QA 

G= 

B^QA""” 

rs] 


(2.17) 


(2.18) 


(2.19) 


( 2 . 20 ) 


2.1.2 Controlador MPC com Trajetória de Referência - MPC-RS (Model 
Predictive Control with Reference System) 


Essa estratégia é fundamentalmente uma transformação de um problema de 
otimização com restrições em um problema de programação quadrática (MONTANDON, 
2005). A experiência do grupo de pesquisa nessa classe de controladores e implementações 
anteriores na planta piloto de controle de nível e pH foram as grandes motivações para o 
estudo dessa técnica. O objetivo foi utilizar conhecimentos prévios, validação das rotinas de 
operação em tempo real e aquisição de sensibilidade no sistema experimental. Não há nessa 
formulação consideração de incerteza, tratando-se de um controle MPC convencional cujos 
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resultados formarão uma base eomparativa em termos de simulação para os eontroladores 
IHMPC e RMPC. Essa formulação será denominada nesse doeumento eomo MPC-RS (Model 
Predictive Control with Reference System). 

A téeniea de eontrole parte da linearização jaeobiana do modelo dinâmieo do 
sistema (Eq. 2.21-2.22). As matrizes A^_j, e são dadas pela Eq. 2.23-2.25. As 
variáveis x'(k), u'(k) e y^(k) passam a ser variáveis desvio em relação ao ponto 
imediatamente anterior ao ponto eorrente de operação. Para simplifiear a simbologia, fez-se, 
A = , B = e C = C]j_i ef = Além disso, adotou-se, também, para efeito de 

simplifieação, Xj, = x'j,, = u'j, = y'j,, eonforme se apresenta nas Eq. 2.26-2.27. 


dx X 

dt 

(2.21) 

y = h(i) 

(2.22) 

A 

^ ^ / C/X-^ i-x(k-l),u=u(k-l) 

(2.23) 

B ] 

^ ^ /ÔU'^ i=x(k-l),u=u(k-l) 

(2.24) 

c 

^ ^ k=x(k-l) 

(2.25) 

dx^ 

-= A, ,x' + B, ,u +f , 

dt 

(2.26) 

y =Ck-iX 

(2.27) 


A diseretização do sistema resulta na forma da Eq. 2.28. As matrizes diseretas de 
espaço de estado (O, T) são então definidas pela Eq. 2.29. Se, e somente se, A for não 


singular e subtraindo-se Xk+i e Xk, ehega-se à Eq. 2.30. 
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=a)x,+TBu,+^f (2.28) 

a) = e^^^ e T = A“^(e^^^-l) (2.29) 

^k+i =^k(2-30) 

Ax,^ = - Xj, j que é o mesmo que AXj^ = x^, uma vez que Xj, é variável desvio e 

Au,^ = Uj^ ““k-i ““k - ^ avaliação dessa formulação sujeita aos horizontes de eontrole e de 
predição, juntamente eom a trajetória de referêneia levará ao algoritmo do eontrolador, 
expresso pela função objetivo a seguir. Para maiores detalhes ver Montandon (2005). 

A trajetória de referêneia da variável eontrolada é definida eomo uma equação 
difereneial linear de P ordem, apresentada na Eq. 2.31. Kj^g serve eomo um parâmetro de 

Op 

sintonia do eontrolador, y é o valor desejado da variável eontrolada, y é a trajetória de 
referêneia da variável eontrolada e y é a predição da variável eontrolada. 


dy 


RS 


dt 




(2.31) 


A eseolha de uma trajetória de referêneia de P ordem se deve, primeiramente, à 
simplieidade matemátiea do problema e, segundo, pelo fato de reduzir a possibilidade de 
oeorrêneia de overshoot no desempenho do eontrolador. No entanto, uma vez que uma 
trajetória de referêneia proporeional não garante a eliminação de offset, um estudo foi 
realizado por Montandon (2005) a fim de resolver este problema. 

A função objetivo do problema de otimização, apresentada a seguir, já expressa a 
forma adequada para eliminação de offset e é definida na Eq. 2.32. 
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min 

Au(k),-",Au(k+Nn,))t(k), 


J = 

,k(k+Np) 


-(au^RAu + ).^S). 


(2.32) 


são parâmetros de sintonia do eontrolador. X é o vetor da 

variável de folga do eontrolador e pode ser eserita eomo uma função objetivo de um problema 
de programação quadrátiea padrão dadas pelas Eq. 2.33-2.36, de forma que o mesmo se 
eneontra sujeito a restrições de limites para as variáveis manipuladas, além da trajetória de 
referêneia: 



Em que: 



(2.34) 


(2.35) 






o^ 

(N 
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2.2 Controle Robusto 


O controle bem sucedido requer que se detenha algum conhecimento das 
características básicas do sistema em questão. Os parâmetros do controlador são ajustados 
com base nesse conhecimento e, quanto mais conhecimento se possui, melhor é o 
comportamento da malha fechada. Várias são as técnicas de controle moderno que exploram o 
referido conhecimento sobre o sistema, destacando-se as técnicas de controle preditivo, que se 
baseiam em um modelo do processo que se pretende controlar. 

Naturalmente, quanto melhor o modelo matemático utilizado para representar as 
características reais da planta, mais eficiente será o esquema de controle baseado nesse 
modelo. Infelizmente, a complexidade do modelo matemático é diretamente ligada à sua 
representação real dos processos. Dessa forma, um controlador que se baseia em um modelo 
incorreto de um dado processo pode levar o sistema a um comportamento completamente 
diferente do previsto. 

A caracterização de um conjunto de incertezas significa encontrar uma 
representação matemática da incerteza do modelo, baseada no conhecimento do sistema. Em 
sistemas MIMO a direcionalidade da planta se configura como um problema adicional. 

As incertezas podem ter muitas origens, como: parâmetros de modelo linear 
aproximados com erro; parâmetros do modelo linear variando devido à característica não 
linear dos sistemas reais associados com o ponto de operação; imperfeições na medida; 
modelos e estruturas desconhecidos em altas freqüências; escolha de modelos mais simples, 
negligenciando dinâmicas como incertezas; controlador implementado diferentemente do 


controlador sintetizado. 
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Como sugeriu Skogestad e Postlethwaite (1997), as principais formas de modelar 

incertezas são: 

• Incerteza paramétrica ou incerteza estruturada, em que a estrutura do modelo é conhecida, 
porém alguns parâmetros são incertos, como apresentada a Eq. 2.37. Gs constitui uma 
classe de modelos associados ao modelo nominal Gn, cujo ganho 13 é incerto, porém 
limitado a i^mín e i3máx- As desvantagens para esse tipo de incerteza são: a demanda de um 
profundo conhecimento sobre o processo, a necessidade de modelo exato, juntamente com 
a impossibilidade de lidar com incertezas de dinâmicas não modeladas, e a exigência de 
perturbações reais no sistema. 

Gs =1^Gn(s);i3mín (2-37) 

• Incerteza de dinâmicas negligenciadas ou não-modeladas, em que o modelo se encontra 
em erro, uma vez que perdeu representação de parte de sua dinâmica, usualmente em altas 
freqüências, que foi deliberadamente negligenciado devido à complexidade de obtenção de 
um modelo matemático ou à dificuldade de lidar com o modelo obtido. A forma usual de 
representação é mostrada na Eq. 2.38, que expressa a norma infinita da incerteza A 
limitada a um número finito que, nesse caso, é a unidade, considerando-se sistemas 
escalados. A Figura 2.4 ilustra o sistema com esse tipo de representação de incerteza. 

NI < 1 (2.38) 

II II oo 



Figura 2.4 - Processo com incerteza dinâmica negligenciada ou não-modelada. 
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• Incerteza distribuída ou ineerteza não estruturada, em que a deserição de ineerteza 
representa uma ou muitas fontes de ineerteza paramétriea e/ou não-modelada, eombinadas 
em uma simples perturbação distribuída de uma estrutura eseolhida. A Equação 2.39 
expressa a ineerteza multiplieativa na entrada e é ilustrada na Figura 2.5. Aqui, Wi é uma 
função peso que, juntamente eom a ineerteza Ai, gera regiões de ineerteza eentrada em Gn, s 
é a variável de Laplaee, j é a variável eomplexa e co representa o domínio da freqüêneia. 
Outros tipos possíveis para essa elasse de ineerteza são: ineerteza aditiva, ineerteza 
multiplieativa na saída, ineerteza aditiva inversa, ineerteza multiplieativa na entrada 
inversa, e ineerteza multiplieativa na saída inversa. Essas formas são semelhantes, porém a 
aplieação em sistemas de eontrole é distinta. Vale ressaltar que é uma eonsideração natural 
0 estudo da ineerteza no domínio de freqüêneia, uma vez que as ineertezas estão 
assoeiadas ás dinâmieas de altas freqüêneias ou transitórias. Esse tipo de ineerteza é uma 
suposição de que todas as ineertezas podem ser globalizadas em um únieo vetor 
(RAMÍREZ; CAMACHO, 2002). 


Gs(s) = Gn(s)(l + Wi(s)Ai(s));Ai(jco)|<l,Vco^||AiL<l 



Figura 2.5 - Proeesso eom ineerteza multiplieativa na entrada. 


(2.39) 


• Eíma alternativa para a aproximação na deserição de ineertezas é eonsiderar uma 
distribuição probabilístiea dos parâmetros eom uma resposta média, eomo apresentou 
Kailath (1980). Esta é a ineerteza estoeástiea. Desse modo, a saída no instante de tempo k 
não pode ser exatamente determinada a partir de dados passados. Considerar a ineerteza. 




33 


nesse easo, signifiea levar em eonta variáveis aleatórias presentes nos proeessos. Porém, 
essa abordagem é de difíeil análise e, até eerto grau, ineoerente eom o estudo de robustez, 
uma vez que, eom as ineertezas estoeástieas, earaeterizadas eomo ruidos na planta, o 
sistema se tornará estável e, na realidade, a planta poderá tomar-se instável. 

• Quando a ineerteza é parametrizada por uma lista finita de possíveis modelos na forma da 
Eq. 2.40, tem-se a forma de politopos, eomo mostrou Bemporad e Morari (1999). Quando 
é permitido que o sistema real varie dentro de uma família de sistemas definida por uma 
lista de modelos, obtém-se a ehamada ineerteza politópiea, eomo apresenta a Eq. 2.41. 
Esse tipo de ineerteza tem sido estudado já há muitos anos e é uma ferramenta flexível e 
versátil para deserever sistemas eom ineerteza e sistemas variantes no tempo (GRIEDER; 
PARRILO; MORARI, 2003). 


n e {n j • • • n } 


yk =cx, 

[a. bJeD 

n=colA,, By] ■■■ [Aj„ b^„]} 


(2.40) 


(2.41) 


o eonheeimento do proeesso é a ehave para a definição das ineertezas. Para sistemas 
eom atuadores e eomponentes meeânieos, pode obter-se um modelo de ineertezas na entrada. 
Sistemas variantes no tempo e não lineares podem ter as suas ineertezas melhor modeladas na 
forma de potitopo. A definição da ineerteza exige eonheeimento do proeesso e integração eom 
a estrutura de eontrole robusto que se pretende utilizar. Gustafsson e Mãkilã (2001) 
apresentam um roteiro para modelagem de ineertezas para aplieação em eontrole robusto. 
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Um controle é robusto se ele for insensível às diferenças entre o sistema atual e o 
modelo utilizado para o cálculo do controlador. Os problemas de robustez estão ligados a 
inúmeros fatores, tais como: defeitos diversos no sistema; sensores; atuadores; 
envelhecimento; restrições físicas; alterações nos objetivos de controle; abertura e fechamento 
de malha, e problemas matemático-computacionais ligados à descrição do problema real de 
controle, algoritmos e convergência. Pode-se, ainda, encontrar o ponto de ótimo desempenho 
que seja excessivamente conservativo e proibitivo para aplicação, em particular, se o pior caso 
raramente acontece. 

A teoria de controle i/oo, na década de 1980, reviveu fortemente a idéia de 
robustez, levando a um desenvolvimento bastante rigoroso e produzindo idéias e métodos 
associados (ÂSTRÕM, 2001). O paradigma chave do controle robusto Hoo é verificar se as 
especificações de projeto são satisfeitas mesmo para a pior condição. 

O estudo de autovalores da planta não oferece a melhor análise da direcionalidade 
da planta. A justificativa para tal fato é que os autovalores medem o ganho para o caso 
especial em que as entradas e as saídas estão na mesma direção. Uma ferramenta para análise 
de direcionalidade é a decomposição de valores singulares (SVD), expressa pela Eq. 2.42. 
Algumas vantagens dessa ferramenta podem ser descritas, como: melhor informação sobre os 
ganhos da planta, as direções obtidas pela decomposição de valores singulares são ortogonais 
e podem ainda ser obtidas de plantas não quadradas. O máximo valor singular também é 
muito útil para o desempenho robusto. Zhu e Jutan (1998) apresentam um controlador robusto 
baseado no conceito de SVD. |3 é a matriz diagonal decrescente de valores singulares não- 
negativos (Oi); ^ é a matriz de saída de vetores singulares, e Ç é a matriz de entrada de vetores 
singulares. O índice H indica a forma transposta conjugada (SÂGFORS; WALLER, 1998). 
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G = (2.42) 

O número de eondieionamento (NC) e a matriz de ganho relativo (RGA) são 
ferramentas matemátieas importantes utilizadas para quantifiear o grau de direeionalidade e de 
interação em sistemas MIMO. 

O número de eondieionamento é definido, na Eq. 2.43, eomo uma relação entre o 
máximo e mínimo valor singular. É interessante notar que o eondieionamento do sistema 
depende fortemente do esealamento das entradas e saídas. Se o número eondieional é grande, 
ele pode indiear alguns problemas de eontrole eomo, por exemplo, a presença indesejável de 
mínimo valor singular pequeno (SÂGFORS; WALLER, 1998). Entretanto, na prátiea, o 
sistema poder perder desempenho drastieamente devido a diferenças no tipo ou, mais 
geralmente, á direeionalidade da perturbaçãoAe/pomt. O número eondieional dessa forma 
provideneia uma medida da robustez do sistema em malha feehada eontra a direeionalidade de 
perturbação e setpoint (ZHU; JUTAN,1998). 

NC(G) = -^^ (2.43) 

o(G) 


A matriz de ganho relativo (RGA) é definida na Eq. 2.44 e ofereee informações 
importantes sobre o sistema em todo o espeetro de freqüêneia. Contudo, os pontos mais 
utilizados para análise são os valores de RGA na forma de estado estaeionário e também em 
freqüêneias próximas ás freqüêneias de eruzamento. 

RGA = G,,H,,, eom H = (g“^ f 


(2.44) 



36 


Algumas das propriedades úteis para eontrole, dadas pela RGA, são listadas a seguir: 

• A RGA é um bom indieador da sensibilidade às ineertezas nos eanais de entrada (ineerteza 
diagonal). Pode indiear, também, a impossibilidade de aplieação de desaeopladores; 
Wolovieh e Falb (1976) e Weller e Goodwin (1996) aplieam o eoneeito de RGA no 
desenvolvimento de desaeopladores. Os trabalhos de busea de pareamento adequado de 
variáveis são feitos desde antes do eoneeito de RGA, eomo apresentados por Niederlinski 
(1971); 

• Variações no sinal da RGA indieam a existêneia de zeros do lado direito do plano 
(eonsiderando-se sistema eontinuo); 

• Análise de plantas não quadradas podendo indiear a neeessidade de eliminação de entrada 
extra ou a impossibilidade de eontrole de saída; 

• A RGA pode ser utilizada eomo medida de dominâneia diagonal e pode indiear o 
pareamento adequado para o eontrole; 

• A RGA forneee informações sobre integridade e estabilidade para sistemas eom eontrole 
deseentralizado. 

Embora pieos de sensibilidade e elementos de RGA, dentre outros, sejam 
indieadores úteis de problemas de robustez, eles não indieam exatamente se uma dada fonte 
de ineerteza promoverá instabilidade ou um desempenho pobre do sistema. Isso motiva a 
neeessidade de melhores ferramentas para a análise dos efeitos de modelos de ineerteza. É 
neeessário que se evitem proeedimentos baseados em tentativa e erro ao verifiear a 
estabilidade e o desempenho para um grande número de plantas eandidatas (CHANG; YU, 
1990; GAGNEPALN; SEBORG, 1982; MCAVOY et ai, 2003; SKOGESTAD; MORARI, 
1987; WITCHER; MCAVOY, 1977). 

O eoneeito de estabilidade interna remete diretamente á idéia de eontroladores 


estabilizantes expressos por parametrização de Youla. A estrutura de eontrole preditivo IMC 
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pode ser eompreendida eomo uma forma dessa elasse de eontroladores (MORARI; 
ZAFIRIOU, 1989). 

Considerando-se o eoneeito de estabilidade interna, pode-se enuneiar um 
importante teorema, eujo resultado é a base para estudo de estabilidade robusta e derivação de 
demais formulações. 


Teorema 2,2 (Teorema do Pequeno Ganho): Supondo M e e seja ooO. Então, a 

intereonexão do sistema mostrado na Figura 2.6 é bem posta e intemamente estável para todo 
A(s) e 9í//oo,oomo mostram as Eq. 2.45-2.46: 

||A|| < 1 /a se e somente se ||m|| <l/a (2.45) 

||A|| < 1 /a se e somente se ||m|| <l/a (2.46) 

II Iloo II Iloo ^ ^ 

Prova: ver Zhou e Doyle (1998). 



Figura 2.6 - Intereonexão do sistema para o Teorema do Pequeno Ganho. 

As limitações dos sistemas podem ser impostas por: tempo morto, fase não- 
minima, pólos instáveis, perturbação do sistema, restrições na entrada, atraso de fase e 
ineerteza (HOLT; MORARI, 1985a; 1985b; MORARI, ZAFIRIOU; HOLT ,1987; MORARI, 


1983; HAVRE, 1998). 
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Dado um conjunto de modelos ineertos H e uma série de objetivos de desempenho, 

supondo que G e H é o modelo nominal e K é o eontrolador que estabiliza internamente G, 

então a malha feehada é dita possuir: 
f 

• Estabilidade nominal (NS): se K estabiliza internamente o modelo nominal G; 

• Estabilidade robusta (RS): se K estabiliza internamente todas as plantas em D; 

I • Desempenho nominal (NP): se os objetivos de desempenho são satisfeitos para G; 

• Desempenho robusto (RP): se os objetivos de desempenho são satisfeitos para todas as 
plantas em IT. 


2.3 RMPC {Robust Model Predictive Control) 


Em 1968, Witsenhausen propôs uma téeniea emergente de eontrole baseada na 
otimização da pior eondição do sistema para eomportamento futuro originando, assim, uma 
das prineipais idéias do eontrole robusto (WEINMANN, 1991). 

O eontrole robusto trata de projetar eontroladores que preservem a estabilidade e o 
desempenho apesar das ineertezas do modelo. A estabilidade robusta é um problema sério 
para a aplieação industrial do MPC (QIN e BADGWELL, 1996). Os diversos métodos de 
Lyapunov podem ser utilizados para prever o desempenho de sistema linear eom ineerteza de 
modelo e também podem ser utilizados para projetar eontroladores adequados a essa situação. 
Alguns métodos são revistos por Osburn e Bernstein (1997) e Ding et al. (2007). 

Os modelos no MPC são utilizados para a predição do que irá aeonteeer nas 
trajetórias futuras. A forma apropriada para se deserever ineerteza nesse eontexto é um 
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modelo ou uma família de modelos, que, ao invés de gerar uma trajetória futura, possa gerar 
uma faixa de trajetórias na qual a trajetória do proeesso esteja inelusa quando a entrada for 
aplieada, apesar das ineertezas. Deve-se esperar que, para um bom modelo, essa banda de 
trajetórias seja estreita, assim eomo o eontrário também pode aeonteeer (CAMACHO; 
BORDONS, 1999). 

A forma eomo o sistema irá eneontrar a trajetória desejada será indieada pela 
função objetivo e dependerá dos sinais presentes e futuros de eontrole e da ineerteza. Por 
eonseqüêneia, pode-se notar que a faixa dos eomportamentos possíveis será sempre maior que 
0 eomportamento real da dinâmiea do proeesso. Por essa razão, o eontrolador projetado 
baseado em modelo de ineerteza é sempre eonservador em suas ações (CHRISTEN et al, 
1997). 

Os modelos matemátieos, em espeeial os modelos para eontrole, podem deserever 
a dinâmiea de um proeesso apenas de forma aproximada (WEINMANN, 1991). A maioria das 
téenieas de projeto de MPC neeessita de um modelo de eontrole da planta eom estrutura e 
parâmetros fixos, que serão utilizados no proeedimento de projeto. O eontrole robusto 
eonsidera explieitamente as diserepâneias entre o modelo e o proeesso real. 

As ineertezas podem ser definidas sobre a eapaeidade de predição do modelo 
(CAMACHO; BORDONS, 1999). Como parte dos esforços de modelagem, é neeessário obter 
uma deserição apropriada para a deserição das ineertezas. Por um lado, a deserição de 
ineerteza deve ser justa, não eontemplando modelos extras ou situações irreais, devendo haver 
0 eompromisso entre o realismo e o resultado da eomplexidade eomputaeional da síntese e 
análise do eontrolador. 

Quando o easo usual de ineertezas limitadas é eonsiderado explieitamente, os 
limites das trajetórias preditas podem ser ealeulados e pareee que o eontrole robusto pode ser 
obtido se 0 eontrolador tentar minimizar a função objetivo para o pior easo, eomo mostrado na 
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Eq. 2.47. A função a ser minimizada é o maximum de uma norma, que mede quão bem a saída 
do proeesso segue a trajetória de referêneia. Diferentes tipos de normas podem ser utilizados 
para esse propósito. Ineertezas estruturadas eonduzem ao eaminho natural de análise do pior 
easo e, por eonseqüêneia, do problema mín-máx. O problema de robustez na presença de 
ineertezas não estruturadas é usualmente tratado por introdução de graus de liberdade 
baseados em parametrização de Youla. Essa téeniea não afeta a função de sensibilidade 
eomplementar e permite a minimização da função eusto de robustez na forma do bem 
eonheeido Teorema do Pequeno Ganho. 

minmaxj(u,0) (2.47) 

ueU ee0 


Mayne e Langson (2001) mostram uma forma alternativa para que a robustez 
possa ser adieionada em algoritmos eonveneionais de MPC, por meio da introdução de 
modifieações relativamente simples. Essa transformação eonsidera a solução nominal do 
sistema, o uso de um eontrolador estabilizante do modelo nominal e a solução ótima para o 
problema eonveneional nominal. 

O eomportamento do MPC pode ser muito eomplieado devido ao fato de a ação de 
eontrole ser determinada eomo um resultado de um problema de otimização em tempo real. 
Enquanto a intuição de engenharia pode ser freqüentemente usada na análise do 
eomportamento ou no projeto de sistemas de MPC, a teoria pode desenvolver valiosa ajuda, 
ampliando o julgamento e a intuição humana no desenvolvimento e implementação de um 
melhor sistema RMPC que possa explorar todo seu o poteneial (NIKOLAOU, 2001). 

A estabilidade é uma propriedade fundamental dos sistemas dinâmieos. Uma 
importante implieação dessa deelaração é que eontroladores RMPC devem resultar em malhas 
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fechadas que sejam estáveis. Como resultado, se um ótimo sistema RMPC é projetado, apenas 
candidatos que pertençam ao conjunto de controladores RMPC estabilizantes devem ser 
considerados. Para assegurar a estabilidade robusta na presença de incerteza do modelo dois 
caminhos podem ser seguidos (NIKOLAOU, 2001). O primeiro é escolher os parâmetros do 
MPC padrão (matrizes Q e R), de forma a atender a Eq. 2.48. 

0^Jopt(k)-Jep.(k + l) (2.48) 

Em outras palavras, a robustez é feita pela tentativa por meio do melhor ajuste do MPC com 
menor modificação em sua estrutura. A vantagem desta aproximação é que a forma do MPC 
padrão é mantida. A desvantagem, por sua vez, é que a otimização on-line pode tornar-se 
irreal e conservadora, dependendo da magnitude da incerteza de modelo no processo. 
Vuthandam et al. (1995) apud Mayne et al. (2000), apresentou um algoritmo no qual a 
estabilidade é garantida pelo atendimento de restrição que implica na redução gradual dos 
pesos da variação dos movimentos de controle (Au). Resultados similares podem ser 
encontrados para MPC com objetivo em tempo real baseado na norma-1. Variações para 
outras formulações de MPC também podem ser encontradas na literatura. Como exemplo, Lee 
e Yu (1997) apresentaram resultados de formulações de MPC empregando o problema de 
otimização on-line na forma da Eq. 2.49. Para esse caso, 0 refere-se a parâmetros de modelo 
do processo que são incertos ou á incerteza propriamente dita. 


min max j[k,u,0] 


(2.49) 
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A idéia dessa formulação é que a ótima entrada para o pior easo possível do 
proeesso seja eomputado a eada passo k. Outra forma é modifiear o problema MFC de 
otimização on-line adieionando a restrição que possa ajudar a estabilizar a malha feehada do 
sistema. Como exemplo de restrição, pode-se eitar a mesma Eq. 2.30. Dessa forma, 
adieionando a restrição no problema de otimização on-line, pode-se assegurar que o 
eontrolador MFC automatieamente perteneerá a um eonjunto de eontroladores estabilizante. 
Cheng e Krogh (1996) apud Nikolaou (2001), desenvolveram a ehamada ''stability 
constrained receding horizon control". Badgwell (1997a; 1997b) desenvolveu a formulação 
de MFC na qual a restrição utilizada no problema de otimização on-line garante a estabilidade 
robusta para a malha feehada do MFC para proeessos lineares estáveis e instáveis. Essa 
formulação é um problema global de otimização eonvexa e tem solução viável que, eomo 
resultado, garante ter uma solução ótima úniea. Tal aproximação tem eomo vantagem o fato 
de que o objetivo em tempo real pode ser formulado de forma a refletir o verdadeiro objetivo 
de eontrole, eonsiderando-se a estabilidade. For outro lado, a desvantagem reside no fato de 
que 0 problema de otimização pode tornar-se bastante eomplexo. Outra questão erítiea é o uso 
de uma família de modelos. O ideal seria a utilização de uma família de modelos bastante 
numerosa euja aproximação linear do proeesso em eada ponto eapturasse as ineertezas. 
Entretanto, a eomplexidade para a solução do problema aumenta na medida em que aumenta o 
número de modelos deseritivos. 

Existem eonsiderações para melhorar o desempenho de MFC, que podem se 
eitadas a seguir: 

• Modifieação na estrutura do MFC para o desempenho robusto; 

• Ajuste do MFC para desempenho robusto; 
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• Desenvolvimento de algoritmos efieientes na solução numériea do problema de 
otimização on-line do MPC, de modo a possibilitar formulações mais eomplexas e 
realistas nos problemas de otimização on-line levando a uma melhoria no desempenho. 

Como em um MPC nominal, eonsiderações espeeiais devem ser levadas em eonta 
a fim de garantir estabilidade robusta. O Teorema 2.3 resume os prineipais fatos. 

Teorema 2,3: Supondo que as seguintes eonsiderações se mantenham para o sistema eom 
ineerteza 0e 0, eontrolador nominal K(x), um dominio de estado terminal Xj, um eusto de 
estágio (p(xk,Uk) e um peso de estado terminal 'P(x), eomo apresentado pelas Eq. 2.50-2.56. 
Então, eonsiderando viabilidade de estado inieial, um eontrolador RMPC usando a seguinte 
formulação mín-máx, expresso pelas Eq. 2.57-2.59, irá garantir estabilidade assintótiea. 


Xk+l = f(Xk, Uk, 0k) 

0 ç XtÇ X 

f(Xk, K(Xk), 0k) e Xt, V XkG Xj, V 0e 0 
cp (0,0) = 0, cp (x,u) = 0, V (x,u) 0 

T' (0) = 0, 'P (x) > 0, V X 0 

'P(f(x,K(x), 0)) - ^{x) < cp ((x,K(x)), V XkG Xt, V 0e 0 

K(x) e U, V X e Xt 

Então: 

Züõ ' ^(^kUlk ’«kU|k ) + Jk ) 

Sujeito a: 

Uk+i|k e U, V 0e 0 

Xk+i|ke X, V 0e0, Xt, V 0e0 


(2.50) 

(2.51) 

(2.52) 

(2.53) 

(2.54) 

(2.55) 

(2.56) 

(2.57) 

(2.58) 

(2.59) 


Prova: (MAYNE et al, 2000). 
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A análise de robustez pode ser definida em duas situações distintas: 

• Análise de robustez de propriedade de MFC padrão projetados para o modelo nominal 
sem levar em eonsideração a ineerteza. Zheng (1999) desenvolve uma ferramenta para 
aplieação em algoritmos MFC que, sem alterar suas earaeterístieas atrativas, apresentam 
estabilidade assintótiea nos easos nominal e robusto, pelo uso de pesos variantes no tempo, 
formulados eomo um problema eonvexo envolvendo LMI; 

• Síntese de algoritmos MFC que são robustos por eonstrução. 

Existe a neeessidade de ferramentas para análise e, devido ao fato de os 
algoritmos padrões de MFC tipieamente requererem um menor esforço eomputaeional, 
tomam-se desejáveis para implementação em tempo real. Como ferramenta de análise, pode- 
se eitar (BEMFORAD; MORARI, 1999): mapeamento de eontração que deriva um eonjunto 
de eondições sufieientes para estabilidade nominal e robusta de MFC; índiee de desempenho 
baseado em norma 1, eondição de estado terminal eomo uma restrição de estabilidade e teste 
de robustez eomo simples satisfação de desigualdades; teste de análise de robustez off-line de 
eontrole de horizonte finito eom restrição que requer a solução de um eonjunto de LMI. 

Duas estratégias são possíveis para o easo da síntese de eontroladores RMFC: 

• Definir um modelo nominal eonsiderando uma perturbação nominal igual a zero e 
otimizar o desempenho nominal - téeniea eonheeida eomo robustez inerente; 

• Resolver o problema mín-máx para otimizar o desempenho robusto. Esta opção é 
eomputaeionalmente muito mais exigente que a primeira. Entretanto, sob eertas 
eonsiderações para as restrições e ineertezas, algoritmos muito efieientes podem ser 
aplieados. For outro lado, a ação de eontrole pode ser exeessivamente eonservadora. Na 
presença de ineertezas, as restrições das variáveis de estado podem ser impostas a todas as 
plantas possíveis ou somente à planta nominal (SCOKAERT; RAWLINGS, 1997). 
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O requerimento mínimo para a malha feehada é a estabilidade robusta, ou seja, a 
estabilidade na presença de ineertezas. Nos vários projetos de MPC a estabilidade robusta é 
eneontrada por dois eaminhos: 

• Indiretamente, pela espeeifieação de objetivos de desempenho e deserição de ineertezas de 
tal forma que a eomputação do eontrole ótimo leve à estabilidade robusta; 

• Diretamente, obrigando a um tipo de restrição de eontração robusta que garanta que o 
estado eneolherá para todas as plantas perteneentes ao eonjunto de ineertezas. 

A generalização do índiee de otimização nominal para o easo robusto expresso na 
formulação mín-máx da Eq. 2.49 pareee ser o eaminho natural. Contudo, para se assegurar a 
estabilidade robusta é neeessário que a ineerteza seja eonsiderada eomo variante no tempo 
(Eq. 2.60). Esse fato adieiona um eonservadorismo que torna a ação de eontrole proibitiva 
para aplieação (ZHENG, 1999). 

minmax J(k,u. ,0, ) (2.60) 

U€U 0G0 


A restrição de eontração apresentada na Eq. 2.61 é uma restrição de estabilidade 
que força a eontração do estado. Para se atingir a estabilidade robusta é neeessário que o 
estado se eontraia para todos os modelos perteneentes ao eonjunto de ineertezas (ZHENG, 
1999). 


“■k+llk 


< a X, ,a < 1 


(2.61) 


O eonjunto terminal elipsoidal invariante tem sido proposto no eontexto nominal 
eomo relaxação da restrição de igualdade terminal. Tal téeniea pode ser estendida para a 
formulação de MPC robusto pelo uso de téenieas eomo LMI (KOTHARE et al, 1996). 
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A predição em malha aberta avalia a direção do sistema para o pior easo de 
perturbação. Isso pode ser muito eonservador uma vez que os efeitos da perturbação são 
reduzidos na malha feehada. Logo, uma forma de se melhorar o eomportamento do sistema é 
realizar a predição eonsiderando a malha feehada (LEE; YU, 1997). Para uma larga faixa de 
ineertezas de modelo, Kothare et al (1996) projetaram, a eada instante de tempo, uma lei de 
eontrole em realimentação que minimiza o “pior easo” de uma função objetivo de horizonte 
infinito, sujeita às restrições de entrada e saída. O problema foi transformado em uma 
minimização do limite superior na função objetivo do pior easo em uma otimização eonvexa 
envolvendo EMI resultando em um algoritmo MPC robustamente estável. Por um lado, a 
formulação usando predição em malha feehada reduz o eonservadorismo. Por outro lado, a 
exigêneia de ineerteza variante no tempo pode ser eonservadora para algumas aplieações. Essa 
formulação será deserita na seção 2.4.1. 

Outra aproximação para a questão da estabilidade robusta é um problema 
aumentado eom restrição adieional no algoritmo MPC nominalmente estável de forma a 
garantir a estabilidade robusta. Essa formulação possui duas vantagens sobre a idéia mín-máx. 
Primeiramente, o desempenho pode ser melhorado quando o modelo nominal se toma tão 
próximo da planta atual, resultando num problema de otimização que pode ser 
signifieativamente menos eomplexo, e a outra vantagem é a poteneial estabilização de 
sistemas instáveis em malha aberta. 

O problema do desempenho robusto raramente é eonsiderado analitieamente, mas 
é avaliado tipieamente por simulações numérieas. Se vários eontroladores forem investigados, 
haverá uma grande exigêneia eomputaeional. Riehards e How (2004) ofereeem uma 
eontribuição para a questão pela formulação de um método analítieo para predição de 
desempenho de um sistema eontrolado por MPC usando a forma de MPC robusto, habilitando 
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uma rápida investigação dos efeitos nos parâmetros do eontrolador sem o uso intensivo da 
simulação. 

A estabilidade robusta garantida reduziria signifieativamente o ajuste e teste de 
algoritmos de MPC. Muitos dos algoritmos industriais disponiveis atualmente sofrem das 
limitações inerentes desde a origem dessa teenologia eom DMC e IDCOM. Esses problemas 
ineluem: eseolha de modelos limitados, realimentação sub-ótima, perda de estabilidade 
nominal e solução sub-ótima ou inefieiente da otimização dinâmiea (QfN; BADGWELL, 
2003). 


2,3.1 Formulação Clássica de MPC Robusto - O problema Mín-Máx 


Considerando que em eapítulos posteriores deste trabalho será feito o uso de 
formulação elássiea de MPC robusto, sua deserição será feita nesta seção. O Algoritmo 
adotado é deserito em Eee e Cooley (2000) e, por simplieidade, será apresentado e trabalhado 
apenas usando a versão para sistemas estáveis. Este algoritmo foi eseolhido por ser 
representativo do eontrole robusto na forma mín-máx. 

O MPC min-máx é definido de forma que o movimento de eontrole implementado 
Uk, no instante k de amostragem, é dado pelo primeiro elemento Uk|k de 
H:{Uk|k, Uk+i|kv • •, Uk+Nm-i|k}, que é a solução do problema de otimização dado pelas Eq. 2.62- 
2.66. O horizonte de eontrole neste easo é igual ao horizonte de predição. Essa formulação é 
definida eomo RMPC-MMTV (Robusí Model Predictive Control with Min-Max Optimization 


and Time Variante Uncertainty). 
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mm 


max 


,e„..A 


w=i: 




k+i 


X IN^- 

■ n 
J=0 


N„-l T Q 

u k+jRu 


k+j 


(2.62) 


Sujeito a: 


Xku =A‘z,+X‘.jA‘ B(0,^j^j)u,^j_;,i = l,...,oo 

(2.63) 

Uk+i e U, i =0,...,Nm-l 

(2.64) 

Uk+Nm+i = 0, i > 0 

(2.65) 

0k+, e 0, i=0,...,Nm-l 

(2.66) 


A definição dos termos neeessários para a eompleta formulação do problema eneontram-se a 
seguir: 


Zk+i =AZk +B(0j,)Uk 


(2.67) 


Definição: Diz-se que uma função A: é affine se existe uma função linear 

L: e um vetor b em tal que A(x) = L(x) -i-b V xe 

B são funções affine de 0, e 0 é um vetor desoonheeido que pertenee ao subespaço eompaeto 
0. As entradas são restritas e perteneentes ao espaço eompaeto e eonvexo U, que ineluem a 
origem eomo ponto interior. Usando uma aproximação similar eomo a utilizada no problema 
IHMPC, deserito na seção 2.1.1, é possível eserever na forma apresentada nas Eq. 2.68-2.81. 
Resultados dessa téeniea podem ser vistos em Silva et al (2006) e Silva et al. (2004). 


mm max 

Heâ” ©ee” 






( 2 . 68 ) 


Sujeito a: 
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+B(0JS^ 


Em que: 

^k+l ~ [^k+l ■■■ ^k+Nm- 


'^k+l~[*^k ■■■ *^k+Nm-l] 

0w=K - 8™-,]' 

6k=[ôi ■■■ ôdim(u)F 


A = 


H = diag 


Q Q Q 


Q = Q +A^QA 


A = diag 


R ...R 

Nm 


u"*” =u. .u 

x"*- 




B(0,) 0 

AB(0,) B(0,,,) 


0 

0 


A"”'-*B(0,) A""--'B(0,,O - B(0,,, ) 


B(0J = B 


l + ô‘ 0 
0 l + ô" 


-1- g<>m(“) 


(2.69) 

(2.70) 

(2.71) 

(2.72) 

(2.73) 

(2.74) 

(2.75) 

(2.76) 

(2.77) 

(2.78) 

(2.79) 

(2.80) 

(2.81) 


0 
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Para se ter o algoritmo de eontrole robustamente estável é neeessário formular o 
projeto do eontrolador usando a deserição de ineerteza variante no tempo. Os maiores 
problemas para essa abordagem residem nos seguintes fatos (LEE e COOLEY, 2000): (a) o 
grande esforço eomputaeional; (b) a natureza não eonvexa do problema de otimização. 

A fim de minimizar essas difieuldades, estes autores propuseram uma formulação eonvexa 
usando a estrutura mín-máx eom ineerteza eonstante no tempo - RMPC-MM (Robust Model 
Predicíive Control with Min-Max Optimization and Time Invariante Uncertainty). 
Basieamente, a diferença deste algoritmo para a forma RMPC-MMTV é que a Eq. 2.80 é 
reeserita na forma da Eq. 2.82 - a ineerteza, nesse easo, é eonstante durante o horizonte de 
predição. Entretanto, o eontrolador obtido ainda demanda um grande esforço eomputaeional e 
a formulação de ineerteza invariante no tempo não é robustamente estável. 


B(0J = 


B(0k) 

AB(0J 


0 

B(0k) 




0 

0 

B(0J 


(2.82) 


Alamo et al. (2008) apresentam um algoritmo de RMPC em formulação min-máx 
baseado no problema de rede. A implementação eomputaeional eresee de forma polinomial 
eomo 0 horizonte de predição, de forma análoga ao que oeorre no proeessamento dos 
problemas em LMI. O problema de otimização mín-máx é transformado em um únieo 
problema QP. A estabilidade do sistema eom esse eontrolador se mantém preservada. Alamo 
et al. (2008) mostram também o desenvolvimento de RMPC assoeiado á transformação do 
problema mín-máx em um problema de otimização quadrátieo e tratável do ponto de vista 
matemátieo-eomputaeional. O eontrolador também mostra a manutenção de earaeterístieas de 


estabilidade e de desempenho. 
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2.4 Desigualdades Matriciais Lineares e o Controle Robusto 

A idéia básica das inequações matriciais lineares (LMI - Linear Matrix Inequalities) 
em controle é interpretar um dado problema de controle como um problema de programação 
semidefinida (SDP), ou seja, um problema de otimização com objetivo linear e restrições 
positivas semidefinidas envolvendo matrizes simétricas que são affme nas variáveis de 
decisão (para maiores detalhes ver Ghaoui e Niculescu (2000)). As técnicas de LMI surgiram 
como uma ferramenta de análise em várias áreas da engenharia de eontrole, identifieação de 
sistemas e projeto estrutural (GAHINET et al, 1995). 

O iníeio da aplieação de LMI na análise de sistemas dinâmieos aeonteeeu há mais 
de eem anos atrás, quando Lyapunov publieou seu trabalho que passou a ser eonheeido, desde 
então, eomo a Teoria de Lyapunov (BOYD et al, 1994). O referido autor mostrou que a 
equação difereneial apresentada na Eq. 2.83 é estável, isto é, todas as trajetórias eonvergem 
para zero se, e apenas se, existir uma matriz positiva definida P tal eomo a Eq. 2.84. As 
exigêneias P > 0 e Eq. 2.85 são ehamadas de desigualdades de Eyapunov em P, sendo esta 
uma forma espeeial de uma EMI. Eyapunov mostrou que a EMI dada pela Eq. 2.84 pode ser 
resolvida explieitamente, eseolhendo-se qualquer matriz Q = > 0 e se resolvendo o sistema 

de equações lineares na forma da Eq. 2.84 para P, que é garantida ser positiva definida se o 
sistema da Eq.2.85 for estável. 

^x(t) = Ax(t) (2.83) 

a’^P + PA<0 (2.84) 


A^P + PA = -Q 


(2.85) 
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Em 1984 N. Karmarkar introduziu um novo algoritmo de programação linear que 
resolve programas lineares em tempo polinomial, eomo o método elipsoidal, porém de forma 
muito mais efieiente na prátiea. Em 1988 Nesterov e Nemirovskii desenvolveram o método 
dos pontos interiores aplieado diretamente a problemas eonvexos envolvendo EML 
Uma desigualdade de matriz linear tem a forma da Eq. 2.86: 

A 

L(x) = Lo+Z[^lXiLi>0 (2.86) 

X é a variável e Li = Li^ perteneente ao eonjunto de números reais de ordem apropriada, 
i=0,...,M são dadas. A desigualdade signifiea que L(x) é positiva definida, isto é, u^L(x)u > 0 
para todo u pertenee ao eonjunto de números reais de ordem apropriada diferente de zero. 
Pode-se observar que a desigualdade é equivalente a um eonjunto M de desigualdades 
polinomiais em x. As EMl também podem ter a forma não estrita, eomo apresentada na Eq. 
2.87. 

L(x) > 0 (2.87) 

As EMl estritas e não estritas são proximamente relaeionadas, porém algumas 
eonsiderações são neeessárias. A EMl da Eq. 2.87 é uma restrição eonvexa em x, isto é, o 
eonjunto {x | L(x) >0} é eonvexo. Essa forma representa uma larga variedade de restrições 
eonvexas em x. Múltiplas EMl, L*'\x) > 0,..., L''“\x) > 0, podem ser expressas eomo uma 
úniea e simples EMl na forma de diag(L*^'\x) > 0,..., L^^^^^x)) > 0. Quando as matrizes Li(x) 
são diagonais, a EMl L(x) > 0 é apenas um eonjunto de inequações lineares. Inequações não 
lineares são eonvertidas na forma de EMl usando o eomplemento Sehur, eonforme será 
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descrito a seguir: a idéia básica é mostrada nas Eq. 2.88-2.89, com: Q(x) = Q(x)^, 
R(x) = R(x)^ e S(x), funções lineares em x. A forma matricial equivalente é mostrada na Eq. 
2.90. Resolver um problema em LMI significa determinar se o problema é viável ou não, e, 
sendo viável, calcular o ponto de viabilidade com um valor objetivo minimo global dentro de 
uma exatidão pré-especificada. 


R(x) > 0 

Q(x)-S(x)R(x)“^S(x)^ >0 
Q(x)>0 

R(x)-S(x)'^Q(x)“lS(x)>0 


( 2 . 88 ) 

(2.89) 


Q(x) S(x) 
S(x)T R(x) 


(2.90) 


Os problemas padrões apresentados a seguir são tratáveis tanto do ponto de vista 
teórico, quanto pelo lado prático. Eles podem ser resolvidos em tempo polinomial, o que 
significa que eles possuem baixa complexidade computacional. Eles também podem ser 
resolvidos na prática de forma bastante eficiente e atingem o ótimo global com critérios de 
parada não heurísticos. E possível provar que o ótimo global pode ser obtido dentro de uma 


precisão pré-especificada. 
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2,4,1 Formulação Clássica - LMI-RMPC 


A técnica de controle MPC tem recebido muita atenção devido às suas muitas 
vantagens. MPC pode facilmente lidar com sistemas com restrição e variantes no tempo com 
bom desempenho. Entretanto, uma das dificuldades da técnica é a incorporação de modelos de 
incerteza da planta. De forma a resolver a questão, Kothare et al. (1996) apresentaram uma 
nova abordagem de MPC robusto que permite incorporar explicitamente a descrição das 
incertezas da planta na formulação do problema na forma de LMI (LEE e PARK, 2008). O 
objetivo de tal formulação é projetar, em cada ciclo, a lei de controle em realimentação de 
estado que minimize o pior caso na função objetivo, considerando o sistema linear variante no 
tempo na forma da Eq. 2.91. A Equação 2.92 mostra a função objetivo considerando-se o 
horizonte infinito. O problema de otimização é apresentado na Eq. 2.93. Definindo uma 
função quadrática na forma V(x) = x^Px com P > 0, tal que V(x) é um limite superior em J(k), 
como mostra a Eq. 2.94, o problema é resolvido encontrando uma lei de controle linear de 
realimentação de estado dada pela Eq. 2.95, tal que V(xk|k) é minimizado. Supondo que não 
exista restrição de estado e de entrada, e que o modelo de incertezas seja definido pelo 
conjunto convexo (co) definido pelos vértices [Ai, Bi], mostrados na Eq. 2.96, qualquer planta 
[A,B]en pode ser expressa como na forma da Eq. 2.97. Nessas condições, o MPC robusto 
pode ser transformado no problema LMI das Eq. 2.98-2.100. Alguns resultados da aplicação 
dessa técnica podem ser encontrados em Silva et al. (2004), Silva et al (2005), Silva et al. 


(2007) e Hu e Yuan (2008). 
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Yk =CXk 

[a, bJeD 

(2.91) 

J.(k) = Lr.«(*wfQ*wn +“k.rR“w) 

(2.92) 

min max J(k) 

“k+iik-iíoK+i 

(2.93) 

[Ak+i Bk+ifn,i>0 

(2.94) 

Com os termos definidos eomo a seguir: 


'^k+ilk ~ ^^k+i|k 

(2.95) 

n = eo{[Aj,BJ,[A2,B2],...,[A^,B^]} 

(2.96) 

[A,B] = ^“ aJA_,BJ eomai>0e =1 

(2.97) 


O problema de eontrole robusto formulado em LMI é expresso a seguir: 


min T 

x,W,T 


Sujeito a: 


1 


z 


k|k 



>0 


W 

WA^ +T^B^ 

1 1 

WQ‘ " 

TTr1/2 

A W + B T 

1 1 

Q 

0 

0 

q' 'w 

0 

xl 

0 

r1/2t 

0 

0 

xl 


(2.98) 


(2.99) 


( 2 . 100 ) 


Uma vez que o problema LMI seja resolvido, a lei de entrada pode ser obtida pela 
Eq. 2.101. Problemas que envolvam restrições de estado, restrições de entrada e outros tipos 
de ineerteza também podem ser formulados e resolvidos eomo problemas LMI, eomo mostrou 
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Kothare et al. (1996). Vale ressaltar aqui que para eada eielo de otimização é obtida uma 
matriz estaeionária de ganho. A prova de estabilidade pode ser derivada dos resultados 
apresentados por Bemussou et al (1989) apud Kothare et al (1996), para estabilização 
quadrátiea de ineertezas politópieas de sistemas eontínuos e suas extensões para o easo 
disereto. 


K = WT'^ (2.101) 

A utilização de algoritmos de eontrole eom aplieação de LMl tem se popularizado bastante. 
Kanev e Verhaegen (2006) mostram um eontrolador robustamente estável baseado em MFC 
de horizonte finito usando LMl para sistemas euja ineerteza não é restrita a alguma elasse de 
representação de ineerteza. Feng, Wang e Poh (2007) desenvolveram um eontrolador 
utilizando LMl aplieado a sistemas eom ineerteza formuladas na forma de LFT {Linear 
Fractional Transformation). Esse eontrolador parte das idéias apresentadas por Kothare em 
1996, direeionando para sistemas eom ineerteza estruturada variante no tempo. Foram 
introduzidos no algoritmo variáveis de folga e parâmetros dependentes de funções de 
Lyaponov. Essa formulação é eapaz de estabilizar sistemas para uma região de eondições 
inieiais maiores e sob fortes restrições nas entradas. Entretanto, a formulação, eonsiderando as 
variáveis de folga no problema de otimização em tempo real, resulta em um aumento 


eonsiderável da demanda eomputaeional. 
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2.5 RMPC com Incerteza na Forma de Perturbação Aditiva 


Lõfberg (2003) apresentou outra abordagem de MPC robusto utilizando LMI e, 
nesse easo, utilizando ineerteza na forma de perturbação externa (RMPC-AD - Robust Model 
Predictive Control with Additive Disturbances). Para esse easo, o algoritmo utiliza o sistema 
na forma linear, disereto e invariante no tempo dado pela Eq. 2.102. Os vetores x, u, d e y são 
0 vetor de estados, o vetor de entradas, o vetor de perturbação externa e o vetor de saídas, 
respeetivamente perteneentes ao eonjunto de número reais de ordem apropriada. Além da 
dinâmiea, o sistema possui restrições de estado e de entrada na forma das Eq. 2.103-2.104. A 
ineerteza, nesse easo, é eonsiderada na forma de uma perturbação externa limitada. 

Xku =Ax,+Bu,+Bgd, 

yk =Cx, 

u eU, U = {u: EuU < fu}, eom 0 e 0 
X e X, X = {x: ExX < fx}, eom 0 e 0 


( 2 . 102 ) 

(2.103) 

(2.104) 


O índiee de desempenho de horizonte infinito é aproximado pelo índiee de 
desempenho finito da Eq. 2.105, eonsiderando-se um horizonte de eontrole N^, e eonforme 
efetuado anteriormente para os eontroladores da elasse de horizonte infinito. O horizonte de 
predição é igual ao horizonte de eontrole. As matrizes Q=Q^>0eR = R^>0 são matrizes 
de ponderação positivas definidas. A fim de se eoloear o problema de otimização na forma 
adequada para programação quadrátiea, são definidos os vetores de saídas futuras, entradas 
futuras e estados futuros nas Eq. 2.106-2.109. Os estados preditos e as saídas podem ser 
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escritos convenientemente na forma da Eq. 2.110. As definições necessárias se encontram nas 
Eq. 2.111-2.114. 


mmmax = (yk^.i/QYk^.ik + «k..|k^R«k+.|k) (2-105) 

^Imi “ (y k|k yk+l|k •" y k+N„-l|k ) (2.106) 

XM=(xk|k ... (2.107) 

UM=(uk|k »L|k - «U^-iikT (2-108) 

D,„-(d;k d;:,,,, ... (2.109) 
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AB 

B 0 


0 0 0 0 

Bg 0 0 0 

Bg = ABg Bg 0 0 

a""” "Bg ABg Bg 0 


( 2 . 110 ) 


( 2 . 111 ) 


( 2 . 112 ) 


(2.113) 
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C 0 0^ 

0 0 
0 0 


(2.114) 


As matrizes Qimi e Rimi são mostradas nas Eq. 2.115-2.116 e possuem ordem 
apropriada no eonjunto dos números reais. Como as restrições são eseritas na forma de 
desigualdades lineares, elas podem ser reeseritas na forma das Eq. 2.117-2.118. 
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Imi 


^R 0 0^ 

0 ■•. 0 
^0 0 Ry 
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0 0^ 
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cr 


cr 


(2.115) 


(2.116) 


(2.117) 


(2.118) 


Nesse ponto, o problema de otimização pode ser re-eserito na forma da Eq. 2.119. 
Esse problema é então eserito na forma de LMI dado pela Eq. 2.120. 


min 


Sujeito a: 
E U<F 

“iml “ 

E U<F 


(2.119) 
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min r 

r,u,„, 


Sujeito a: 
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E 

u 

E B 

V ^Imi 


Imi 




F -E A, .x... , 

V *lmi *lmi 1™ k|k ) 


( 2 . 120 ) 


Considerando a transformação da perturbação externa na forma de uma ineerteza 
em bloeo diagonal A, após manipulações algébrieas, obtém-se a forma final dessa formulação 
em LMl dada pela Eq. 2.121. O termo E é o esealar a ser minimizado, Uimi é o vetor de ações 
de eontrole, y é a máxima ineerteza e T é a matriz diagonal de Xi, parâmetro relaeionado eom 
as ineertezas. Torna-se importante ressaltar que, para este easo, a restrição de estado é 
robustamente satisfeita se E^(Ax^ii^ -t-BUj^^^j-l-y < for satisfeita. A Eq. 2.122 expressa a 

LMl referente à restrição de estabilidade que, adieionada ao problema, garante a estabilidade 
robusta. Resultados da aplieação dessa téeniea podem ser vistos em Silva et al. (2005). 


min r 


Sujeito a: 
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2.6 Considerações Finais 


A evolução dos sistemas de eontrole e de toda a teoria assoeiada sempre esteve 
vineulada à resolução de problemas prátieos e relevantes para a humanidade. A utilização de 
ferramentas matemáticas levou o controle de sistemas a utilizar-se de toda uma teoria 
previamente elaborada, inclusive contribuindo para o desenvolvimento de novos métodos e 
linhas de pesquisa. Se, por um lado, os problemas chegaram então a uma solução calculada, 
por outro, perdeu-se, em algum ponto dessa evolução, parte da sensibilidade de engenharia. 

Com informações cada vez mais disponíveis e, máquinas e equipamentos cada vez 
mais acessíveis, permitiu-se a propagação das possibilidades de automação. 

Diversos ramos das ciências convergem para a multidisciplinaridade e a 
transdisciplinaridade. Um exemplo é a utilização cada vez maior de teorias ligadas aos 
sistemas físicos e biológicos juntamente com ferramentas matemáticas e computacionais em 
aplicações de controle de processos. A visão integral dos sistemas é uma forma de recuperar a 
sensibilidade sobre os processos. Essa mesma convergência é responsável pelas diversas 
possibilidades de construção de algoritmos de controle preditivos. 

Algumas técnicas de controle ainda são penalizadas por exigirem mão-de-obra 
especializada, e, dentre elas, encontra-se o controle robusto. No entanto, trata-se de uma linha 
de pesquisas de grande profundidade que realmente impressiona aqueles que resolvem estudá- 
la com afínco. 

O estudo resultante deste capítulo é muito maior do que se apresenta. Constrói-se 
todo um universo de idéias que, muito além de descrever sistemas, trouxe entendimento para 
questões práticas e formação para resolução de problemas industriais. 



62 


O objetivo deste eapítulo foi apresentar, de forma sueinta, eoneeitos esseneiais ao 
estudo de eontroladores preditivos robustos. Ao longo do eapítulo foi possível situar o 
eontrole preditivo robusto no eontexto dos sistemas de eontrole e dos eontroladores 
preditivos. As definições de ineerteza e robustez são impreseindíveis para a análise e projeto 
de eontroladores robustos. Levando-se em eonsideração as informações obtidas eom a revisão 
bibliográfiea deste eapítulo, serão propostos algoritmos de RMPC no Capítulo 3. 



CAPITULO 3 
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3 DESENVOLVIMENTOS EM RMPC 


3.1 Introdução 

As idéias a serem desenvolvidas neste eapítulo podem ser resumidas eonforme as 
eonsiderações dispostas a seguir: 

• RMPC-MMR (Robust Model Predictive Control with Min-Max optimization, time 
invariant uncertainty and constRained stability): MPC robusto baseado em otimização 
min-máx eonsiderando ineerteza na entrada (matriz B) invariante no tempo não é 
robustamente estável. A modifieação, nesse easo, será a adição da restrição de robustez. 
Esse algoritmo foi baseado nas formulações de MPC robusto baseado em otimização min- 
máx - Lee e Cooley (2000) - RMPC-MMTV e RMPC-MM e nas idéias para forçar a 
robustez segundo Badgwell (1997a); 

• LMI-RMPCRli^oò^/xt Model Predictive Control with stability constRained lype Lyaponov 
and uncertainty as yolytoye in LMI form) : MPC robusto utilizando LMI baseado na 
formulação apresentada por Lõfberg (2003), que usa ineertezas na forma paramétriea. A 
ineerteza adotada foi a forma de multimodelos. Esse algoritmo foi baseado utilizando-se 
MPC robusto baseado em LMI - Lõfberg (2003) - RMPC-AD. 
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3,1,1 Restrição de Robustez 


Existem muitas estratégias para tornar um algoritmo MFC robusto. De forma a 
pesquisar algoritmos de eontrole robustamente estável, as idéias apresentadas por Bagdwell 
(1997a) irão forçar a robustez na forma de uma restrição. Forçar a robustez por uma restrição 
se assemelha a requerer que a função custo seja menor ou igual à função custo de malha aberta 
(isto é, 0 custo obtido implementando a restrição de trajetórias ótimas previamente 
computadas) para todos os valores de parâmetros. Esse requerimento pode ser expresso pela 
Eq. 3.1. S* é uma versão deslocada do valor ótimo de j para o problema de otimização 
min-máx. A vantagem de tal formulação reside no fato de que a restrição de contração dada 
pela Eq. 3.1 se baseia diretamente no índice de desempenho e, por conseqüência, é mais 
intuitivo que a maior parte das restrições de contração. Os resultados de estabilidade robusta 
para os controladores nos quais a restrição for adicionada são confiados inteiramente às 
restrições de contração. 


J.(S„ej<E(S:,0J,V0,e0 


(3.1) 


O MFC robustamente estável é definido de forma que o movimento de controle 
Uk, no instante k de amostragem, é dado pelo primeiro elemento Uk|k de 
S:{Uk|k, Uk+i|kv, Uk+Nm-i|k}, que é a solução do problema de otimização dado pelas Eq. 3.2- 


3.3. 
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Jk =minmaxXo'í'(xk+i.«k+i>^«k+i) (3-2) 

S 0 
Sujeito a: 

«ku|k^U,V0,e0 

Xku|k^X,V0,e0 (3.3) 

Jk(Sk,0k)^Jk(3:,0J,V0,e0 

Será apresentado aqui apenas um esboço da prova, que pode ser feita de forma semelhante às 
provas apresentadas por Lee e Cooley (2000). 

Prova: . Inicialmente é definido o conceito de deslocamento da ação de controle, dado pelo 
operador de deslocamento (rj) dado pela Eq. 3.4 (BAGDWELL ,1997). Note que a restrição 
de contração introduzida é viável e a inequação sempre se mantém, 

Jk(Hk+i) < Jk(Sk+i) = Jk(“D“'í^(xk+L“k+L0k+l) • Dadas as considerações necessárias 
para a função custo, então Jk(Hk,0k)^0 é uma seqüência convergente e os estados e 
entradas irão se anular para k ^ oo , o que faz o projeto robustamente estável.# 

0 I 0 0 

001:0 

0 I 

0 0 . 0 



(3.4) 
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3.2 Formulação de MFC Robusto Mín-Máx com Restrição de Robustez 
Adicionada - RMPC-MMR 


Para o estudo de RMPC a eseolha do problema de otimização mín-máx eom 
ineerteza nas entradas foi feita levando-se em eonsideração que em sistemas reais este tipo de 
problema é relativamente comum, devido à facilidade de simulação e pela natureza das 
mainipulações mecânicas dos sistemas. 

O MPC min-máx é definido de forma que o movimento de controle Uk, no instante 
k de amostragem, é dado pelo primeiro elemento Uk|k de H:{Uk|k, Uk+i|kv, Uk+Nm-i|k}, que é a 
solução do problema de otimização dado pelas Eq. 3.5-3.9 (SILVA et ai, 2004; SILVA et al, 
2005). O horizonte de controle é igual ao horizonte de predição para a formulação descrita. 
Essa formulação é denominada RMPC-MMR (Robust Model Predictive Control with min- 
max optimization, time invariant uncertainty and constrained stability). 


min max 

Ur 0R v^Qr-Nj^-i 

Sujeito a: 

Xku =A‘z,+X‘^jA‘ B(e,^j_j)u,^j_;,i = l,...,oo 

Uk+i e U, i =0,...,Nm-l 
Uk+Nm+i = 0, i > 0 

Sk+iG 0, Í=0,...,Nm-l 


w=i: 




k+i 


j=0 


N„-l T Q 

u k+jRu 


k+j 


(3.5) 


(3.6) 

(3.7) 

(3.8) 

(3.9) 


Por simples manipulação algébrica, o problema de horizonte de controle infinito é 
transformado no problema de horizonte de controle finito expresso pelas Eq. 3.10-3.28. E 
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importante ressaltar que a incerteza neste tipo de algoritmo é um valor a ser obtido pelo 
programa de otimização. Note que z é o vetor de estados reais e x é o vetor de estados 
estimados. 


A 

minmax +S/AH, 

Sujeito a: 

0<l|0| 

Jk(Sk,ej<j,(H:,0j,V0,G0 
e U, i = 0,...,Nin-l 
= 0’ quando i > 0 

^k+i = AZk +B(0 j^)H^ 


(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 




(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 
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B(0J 


B(0,) 

0 

0 


AB(0J 

B(0k) 

0 

(3.22) 

A"'”-*B(0J A''”'^'B(0J •• 

B(0J' 


[l + ô* 0 

0 

1 



B(0J = B 


0 l + ô" 


0 


\ ^ g<lim(u) 


(3.23) 


H = diag 


Q Q Q 


Q = Q + A^QA 
A = diag R...R 

U"*” =u. .u 

N„ 


(3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 


Teorema 3.1: Quando as ações de eontrole do problema de otimização min-máx (3.10), 
sujeito às restrições (3.11-3.15) e às definições (3.16-3.28), são aplieadas ao sistema de 
eontrole na forma Xk+i = Axk + Buk, o resultado será em laço feehado robustamente estável, 
em que: Zk ^ 0, Uk ^ 0 quando k ^ oo para toda a eondição inieial zo e para todos os 
possíveis parâmetros da seqüêneia juj, e U, k>0}. 


Prova: A estabilidade robusta pode ser provada eonsiderando o eusto de uma planta real. O 
valor da máxima ineerteza é determinado pelo problema de maximização e eneontra-se 
obviamente limitado a uma faixa. Sendo assim, o problema de minimização eneontrará a ação 
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V ís* 0)<v Íh* 0 )= 

, ^ , . _£■ j k+l\“k+l’'^/— '^k+l\“k+l’^/ 

de controle para 0 pior caso, na íormade: , . 

vdH:.0)-z;^,Qz,.,-uL,Ru,., 


Como 


Vk 



> 0 é uma seqüência convergente e Zk ^ 0, Uk ^ 0# 


3.3 Formulação de MFC Robusto Baseado em LMI - LMI-RMPCR 


Considerando 0 sistema na forma linear, discreto e invariante no tempo, dado pela 
Eq. 3.29, os vetores x, u e y são 0 vetor de estados, 0 vetor de entradas e 0 vetor de saídas, 
respectivamente pertencentes ao conjunto de número reais de ordem apropriada. Além da 
dinâmica, 0 sistema possui restrições de estado e de entrada na forma das Eq. 3.30-3.31. A 
incerteza, nesse caso, é considerada na forma politópica, como mostra a Eq. 3.32. O horizonte 
de controle é igual ao horizonte de predição para a formulação descrita. 

Xku =Ax,+Bu, 

yk =Cx, 

u eU, U = {u: EuU < fu}, com 0 e 0 
X e X, X = {x: ExX < fx}, com 0 e 0 
[a B]g n 

n=colA, B,] [a„ b„]} 


(3.29) 

(3.30) 

(3.31) 

(3.32) 


O índice de desempenho de horizonte de controle infinito da Eq. 3.33 é 
aproximado pelo índice de desempenho finito da Eq. 3.34, considerando-se um horizonte de 
controle N^. A dedução, neste caso, considera que os horizontes de controle e de predição são 
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iguais. As matrizes Q = > 0 e R = > 0 são matrizes de ponderação positivas definidas. 

A fim de se eoloear o problema de otimização na forma adequada para programação 
quadrátiea, são definidos os vetores de saídas futuras, entradas futuras e estados futuros na Eq. 
3.35. Os estados preditos e as saídas podem ser eseritos eonvenientemente na forma da Eq. 
3.36-3.38. As definições neeessárias se eneontram nas Eq. 3.39-3.41. 


Jk =Zr=o(yk..^Qyk.. +»ku^R»ku) 

J k = Zür* (yku^Qyku+«kun^RUku.i) 

V — r x 

Imi Imi ^ Imi 

^Imi '^Imi^klk 


(3.33) 

(3.34) 

(3.35) 


Com os vetores e matrizes definidos eonforme a seguir: 


Y, . = 

Imi 

(yik 

T 

yk+l|k 

T \T 

y^+N^-llk/ 

(3.36) 


kik 

T 

^k+l|k 

T \T 

•• ^k+N„-l|k/ 

(3.37) 

U,..= 

kik 

T 

'^k+llk 

T \T 

•• “k+N„-l|k/ 

(3.38) 


í ^ I 

A 

(3.39) 



0 

0 

0 
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... 0 


AB 
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AB 

B 0 


(3.40) 
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Imi 


0 0^ 
0 •. 0 
^0 0 Cj 


(3.41) 


As matrizes Qimi, Rimi pertencem ao conjunto de números reais de ordem 
apropriada e são mostradas nas Eq. 3.42-3.43. Como as restrições são escritas na forma de 
desigualdades lineares, elas podem ser reescritas na forma das Eq. 3.44-3.45. 
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Imi 


^R 0 0^ 

0 •. 0 
^0 0 Ry 
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Imi 


0 0^ 
0 •. 0 


0 0 Q 
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0 ^ 
0 •. 0 

0 0 E 


E.. 0 


F = íf^ 

’ ui™ Vu 


E. 


E. 0 
0 • 
0 


uy 

0 ^ 
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0 E 




xy 


f.ú’ 


f/)’ 


(3.42) 


(3.43) 


(3.44) 


(3.45) 


Nesse ponto, o problema de otimização pode ser re-escrito na forma da Eq. 3.46. 
Esse problema é apresentado na forma de EMI dado pela Eq. 3.47. 


minY, .'^Q, .Y, .+U, .'^R, .U 

y Imi ^ Imi Imi Imi Imi 

Sujeito a: 


E 
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“iml 


*lmi 


Um. 
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Imi 


<F 

<F 


“iml 


xlmi 


(3.46) 
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min X 


Sujeito a: 




r E,. ^ 

f 

ulmi 

u< 

E B, , 

V X Imi Imi J 

V 


F , . ^ 

ulini 

-E A, .x,„ , 

*lml 1™ k|k J 


(3.47) 


As Equações 3.48-3.50 expressam a LMI para o problema eom restrição resultante 
que garantem a estabilidade robusta. 


mm y 


(3.48) 


Sujeito a: 

U, 0 R 


-1 
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Imi 


Imi 


>0 


, E -B , 

V xlmi i,lmi J 




, F -E , .A., .x.|. , 

V *lml *1™ ‘>™‘ J 


(3.49) 


1 

v(A.Zk|k+B.u,|J 


(A.Zk|k+B.»k|k) 




.-1 


>0 


(3.50) 


Torna-se importante ressaltar que, para este easo, a restrição de estado é robustamente 
satisfeita se E^(^Xj^ij^ -l-BU) < F^ for satisfeito. A Eq. 3.50 expressa a EMI referente à 

restrição de estabilidade que, adieionada ao problema, garante a estabilidade robusta. 
Resultados da aplieação dessa téeniea podem ser vistos em Silva et al. (2007). Outro aspeeto a 
ser ressaltado é que a ineerteza está eontida na forma de um politopo. Quando a função 
objetivo é eonvexa e o modelo de ineertezas gera um eonjunto feehado de predições, o 
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atendimento dos requisitos dos vértiees do politopo assegura que todos os modelos internos da 
região sejam atendidos, eonforme indiea o Teorema 3.2 (LÔFBERG, 2003). 

Teorema 3.2: Seja C um eonjunto feehado e seja f: C ^ R uma função eomplexa. 
Então, se f atinge o máximo sobre C, ele ehega ao máximo em algum ponto extremo de C 
(Lõfberg, 2003). 

Teorema 3.3: Quando as ações de eontrole do problema de otimização (3.48), sujeito ás 
restrições (3.49-3.450) e ás definições (3.36-3.45) são aplieadas ao sistema de eontrole na 
forma Xk+i = Axk + Buk, resultará em laço feehado robustamente estável. 

Prova: por indução, eonsiderando-se que o problema é viável para X]^_i|i,_l. Como o 
problema foi viável para X]j_i|]j_l, eonheee-se xj^jj^eX. A viabilidade da solução no 
instante k é V = 0. = (A -I- BK)xj^^ji^. Para a eonsideração de que K e P estão eontidos 

em {x : x^Px < l]p > 0, garantindo a viabilidade das restrições de estado e de eontrole. A nova 
restrição Xj,^j g {x : x^Px < l]p > 0, se mantém para U]j|]j = Kx]j|j.por eonstrução. 


Conseqüentemente, o problema é viável no instante k. # 
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3.4 Considerações Finais 


Foram escolhidas formas de incerteza representativas, de fácil implementação para 
aplicações de controle em tempo real, que buscam reduzir o esforço computacional e que 
incorporem novas técnicas matemáticas e computacionais. Basicamente, o objetivo deste 
capítulo foi apresentar algoritmos desenvolvidos e que são robustamente estáveis. Um deles 
foi obtido pela formulação mín-máx, em que estabilidade robusta foi adicionada através da 
restrição de robustez. O outro foi obtido utilizando incertezas na forma de LMI aplicado a 
sistemas com incerteza politópica. Espera-se que a implementação desses controladores possa 
ajudar a reduzir a distância entre a teoria e a prática, contribuindo ainda para a divulgação e 
aplicação de controladores RMPC. No Capitulo 4 serão apresentados os resultados de 
simulação para estas técnicas e no Capítulos 7 os resultados para comparativos em termos de 
simulação com controladores clássicos implementados experimentalmente na unidade de 


controle de nível e pH. 



CAPITULO 4 
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4 SIMULAÇÃO 


4.1 Introdução 


De forma a demonstrar a aplieabilidade dos algoritmos de eontrole propostos no 
Capítulo 3 serão avaliados, neste eapítulo, as respostas e o eomportamento dos mesmos 
aplieados a proeessos eomuns na indústria que, por suas caraeterísticas, exigem de alguma 
forma a robustez, o desempenho ou a estabilidade do eontrolador. Serão estudados proeessos 
que possuem desvio planta/modelo, diserepâneia entre tempo morto da planta e do modelo, 
fase mínima/fase não-mínima, sistema mal eondieionado, modelo linear/planta não linear, 
diserepâncias devido à diseretização do sistema e modelo estável/planta instável. Essas 
respostas serão comparadas aos algoritmos clássicos da literatura apresentados no Capítulo 2: 
IHMPC (MUSKE; RAWEINGS, 1993), RMPC-MMTV (LEE; COOLEY, 2000) e LMI- 
RMPC (KOTHARE et al, 1996). Durante a pesquisa foram simulados diversos ajustes e são 
mostrados a seguir apenas os resultados comparativos e representativos dos estudos, sem que 
isto implique em indicar a melhor condição. Os ajustes dos controladores, como pode se ver 
nas seções a seguir, foram bastante similares, senão idênticos, de forma a tornar a comparação 


mais visível. 
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4.2 CSTR Isotérmico 


Nesse tópieo será avaliado o eomportamento dos eontroladores frente a um 
sistema SISO {Single Input Single Output), não linear, eom desvio planta/modelo inerentes ao 
proeesso de diseretização e presença de fase não mínima em determinadas regiões. Considere 
um reator eontínuo de tanque agitado (CSTR) usando a taxa de diluição eomo variável 
manipulada. O volume do reator é eonstante. A dinâmiea e o eontrole desse reator são 
bastante estudados (HENSON; SEBORG, 1997). O modelo do reator é dado pelo sistema 
apresentado na Eq. 4.1-4.2, que representa as reações da Eq. 4.3-4.5. Esse sistema é eonheeido 
eom reação de Van der Vusse. 


Xj = -kjXj -k3X^ +(Xf. -Xj)u 

(4.1) 

^2 =kiXi-k2X2-X2U 

(4.2) 

k, 

A^B 

(4.3) 

k2 

B^C 

(4.4) 

ks 

2A^D 

(4.5) 


As variáveis de estado xi e X 2 representam as eoneentrações das espéeies A e B, 
respeetivamente, dadas em mol/L, xp expressa a eoneentração de A na alimentação e u 
representa a taxa de diluição (da alimentação dado em h"'). Os parâmetros {ki, k 2 , ks, xp} 
possuem os seguintes valores nominais, respeetivamente: {50 h'\ 100 h'^ , 10 L/(mol.h), 10 
mol/L}. O objetivo de eontrole primário é manter a saída y = X 2 no valor de setpoint desejado. 
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Foi resolvido o estado estaeionário do sistema, apresentado na Figura 4.1. É importante 
ressaltar que este sistema possui multiplieidade de estado estaeionário na ação de eontrole nos 
valores de u = 25 h'' e u = 233,3 e que para uma máxima eonversão de xi, onde u é 
aproximadamente igual a 82 h‘\ a derivada do sinal de eontrole é nula representando um 
ponto de difieil aplicação de estratégias de controle. O programa foi feito de forma a admitir 
simulações de vetores de valores desejados para as saídas e perturbações simultâneas ou não. 
Os casos simulados foram os mesmos para todos os algoritmos a fim de tomar as respostas 
comparativas na mesma base. 



Figura 4.1 - Comportamento do sistema em estado estacionário. 

O sistema não linear foi linearizado segundo a expansão em série de Taylor em 
tomo do estado estacionário de u = 25 h'*, resultando nas matrizes A e B apresentadas nas Eq. 
4.6-4.7. O sistema também foi discretizado com tempo de amostragem de Td = 0,0012 h, 
obtido conforme regra apresentada por Seborg et al. (2004) e expressa pela Eq. 4.8, onde Td é 
0 tempo de amostragem e Ts é o tempo de assentamento do processo. Com este sistema foram 
calculados os controladores MFC e RMPC nos diversos valores de referência. O sinal de 


controle calculado pelo modelo linearizado discreto foi aplicado na planta não linear e 
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contínua sendo retido este valor durante a amostragem. Os estados neste easo foram 
eonsiderados eomo medidos, não sendo projetados observadores de estado. 


A = 


u-kj -2k3Xj 


0 

-u-k. 


B = 


Xp-Xi 


-X, 


^<Td<^ 

15 6 


(4.6) 

(4.7) 

(4.8) 


As restrições consideradas para este tópico são as seguintes: limite máximo de u igual a 
325 h"'; limite mínimo de u igual a zero; limite para variação de u igual a 10 h'*; limite de 
saída igual a 10 mol/L, e limite mínimo da saída igual a zero. Os parâmetros de ajuste para 
todos os eontroladores eom restrição são horizonte de controle de 5 passos, e todas as matrizes 
de peso iguais à matriz identidade de ordem adequada, eonforme indica a Tabela 4.1. A 
incerteza nesse easo é limitada a 10% da entrada e o horizonte de predição é igual ao 
horizonte de controle. A Figura 4.2 mostra as respostas para o sistema linear com desvio 
planta/modelo submetido ao controlador IHMPC com restrição aplicada ao sistema não linear 
eom ineerteza aleatória (MUSKE e RAWLINGS, 1993a). O mesmo vetor de referêneia será 
aplieado a todos os eontroladores a seguir. Por questão de simplieidade, será mostrado o 
eomportamento da variável Xi apenas no primeiro eontrolador mostrado a seguir. Para os 
demais, será mostrada apenas a saída (X 2 ) e a variável manipulada (u), variáveis estas que 
realmente possuem relevâneia para diseussão deste tópieo. 


TABELA 4.1: AJUSTES PARA CONTROLADORES. 


Controlador 

Q 

R 

S 

Incerteza 

IHMPC 

I 

I 

I 

Não se aplica 

RMPC-MMTV 

I 

I 

I 

0,11 

RMPC-MMR 

I 

I 

I 

0,11 

LMI-RMPCR 

I 

I 

I 

0,11 
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TemDOÍh) (a) 



Temooíh) (b) 



Tempo(h) 


Figura 4.2 - Resposta do sistema submetido ao eontrole HíMPC eom restrição e proeesso não 
linear eom incerteza aleatória. 
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A Figura 4.3 mostra as respostas para o sistema com controlador RMPC-MMTV 
(LEE e COOEEY, 2000) calculado sobre a maximização de uma incerteza de ±10% na 
entrada com restrição aplicada ao sistema não linear com incerteza aleatória. 



Temooíh) (a) 




Figura 4.3 - Resposta do sistema submetido ao controle RMPC-MMTV com restrição 
aplicada ao sistema não linear com incerteza aleatória. 


A Figura 4.4 mostra as respostas para o sistema com controlador RMPC-MMR 
calculado sobre a maximização de uma incerteza de ±10% na entrada com restrição aplicada 


ao sistema não linear com incerteza aleatória. 
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(a) 



(b) 


A Figura 4.5 mostra o comportamento do sistema para o eontrolador LMI- 
RMPCR. Neste easo, para a situação de desvio planta/modelo dada por diferenças entre 
modelo/planta e sistema continuo/disereto, o sistema responde adequadamente. 

As Figuras 4.2-4.5 apresentam dinâmica razoável. Pode-se pereeber que o sistema 
não apresenta simetria quando da alteração do valor de referêneia, característica esperada 
eonsiderando-se que neste easo o sistema onde o controlador foi aplicado é não linear e o 
modelo foi linearizado no ponto inieial sem atualização. 
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(a) 


(b) 


A ação de controle apresentada atende aos requisitos de desempenho do sistema. 
Não é notada a presença de offset. Espera-se que o caminho do sistema de u= 233,3 h'' para 
u = 25 h"' possa ser melhorado pela utilização de modelo linear local em u=233,3 h"'. 

A Figura 4.6 mostra as respostas do sistema de forma comparativa, submetido aos 
seguintes controladores: IHMPC, RMPC-MMTV, RMPC-MMR e LMI-RMPCR. 
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Figura 4.6 - Respostas dos diversos controladores aplicados ao processo da seção 4.2. 

As Figuras 4.7-4.8 mostram o controlador projetado usando Ricker (1990), 
algoritmo implementado no Matlab®, com o modelo linearizado aplicado à planta não linear. 
As respostas apresentam grande oscilação, indicando a perda de desempenho. Para o caso de 
desvio planta/modelo, o sistema mostra-se instável. 

As respostas do MPC clásssico (IHMPC), do RMPC clássico (RMPC-MMTV, 
LEE e COOLEY, 2000) e dos RMPCs: RMPC-MMR e LMI-RMPCR, mostraram-se 
similares. Pode-se notar diferentes comportamentos dinâmicos para a elevação de u de 25 h"' 
para 233 h'\ comparando-se com o retorno de 233 h"' para 25 h"', característico da natureza 
não liner do problema de controle. Durante esse estudo não houve a preocupação de se buscar 
0 melhor controlador, e sim de obter um controlador que atendesse de maneira satisfatória as 
exigências do processo em estabilidade e desempenho. 
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(a) 



Figura 4.7 - Sistema submetido ao eontrole MFC (Matlab®): problema regulador. 



(a) 
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O controlador MPC possui características que o tomam mais adequado para 
determinados sistemas, eomo por exemplo, para a multiplieidade de estado estaeionário. A 
adição de restrições é de fácil desenvolvimento e implementação, reforçando a adequação 
dessa téeniea. Os eontroladores propostos no Capítulo 3 apresentaram respostas similares às 
respostas dos controladores clássicos como IHMPC e RMPC-MMTV. 


4.3 Modelo Estável/Planta Instável em Malha Aberta 


Este tópico visa avaliar o comportamento do sistema eonsiderando a incerteza no 
modelo para o caso extremo de planta estável modelo de controle instável e/ou vice-versa 
(WANG e RAWLING, 2004 a e b). Para este easo, apenas as formulações LMI são aplicáveis, 
uma vez que a formulação apresentada para as demais é válida apenas para sistemas estáveis. 
Os modelos utilizados são expressos pela Eq. 4.9-4.10. O tempo de amostragem neste easo é 
de 1 s, 0 horizonte de eontrole (Nm) é igual a 5 com Np é igual a Nm e as matrizes de pesos são 
matrizes de dimensões adequadas, conforme indica a Tabela 4.2. 




1 


z-0,5 


(4.9) 


G 


2 


1 


z-1,6 


(4.10) 


TABELA 4.2: AJUSTES PARA CONTROLADORES. 


Controlador 

Q 

R 

s 

Incerteza 

LMI-MPC 

I 

0,0011 

I 

Não se aplica 

LMI-RMPCR 

I 

0,0011 

I 

0,11 







o objetivo de eontrole neste sistema SISO é que a saída passe de 0 para -1. Para 
aplieação eom sueesso da téeniea LMI RMPCR, faz-se neeessário o uso do sistema na forma 
de Au, como indica Muske (1995). As restrições foram feitas de forma a não influenciar no 
comportamento do controlador. Para este dois casos, a planta foi adotada como o modelo G1 
(estável) e o controlador foi calculado considerando-se o modelo G2 (instável). Os resultados 
pelo algoritmo do toolbox do MPC/Matlab® (RICKER,1990) apresentam-se instáveis para 
qualquer condição de incerteza ou de desvio planta/modelo. Na Figura 4.9 é mostrado o 
comportamento do sistema utilizando o controlador LMI RMPC. Pode-se perceber que a 
resposta atende a condição requerida sem offset e com dinâmica bastante rápida e com 
oscilações tanto na saída quanto na entrada. 


0 

- 0.5 

> -1 

- 1.5 

-2 

0 20 40 60 80 100 120 140 

Tempo(s) 



(a) 


0.5 
0 

- 0.5 

3 

-1 
- 1.5 
-2 

0 20 40 60 80 100 120 140 

Tempo{s) (b) 

Figura 4.9 - Comportamento do sistema para modelo instável e planta estável: LMI RMPC. 
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A Figura 4.10 mostra o comportamento do sistema eom eontrolador LMI-RMPCR 
que apresenta as seguintes earaeterístieas: saída de aeordo eomo o valor de referêneia sem 
offset, variação da entrada oseilante e deereseente, e dinâmiea da saída razoável. 

0 

- 0.2 
- 0.4 
- 0.6 
- 0.8 
-1 

0 20 40 60 80 100 120 140 

Tempo{s) 

0.2 

0.1 

^ 0 

- 0.1 

- 0.2 

0 20 40 60 80 100 120 140 

Tempo(s) (b) 

(b) 

Figura 4.10 - Comportamento do sistema para modelo instável e planta estável LMl-RMPCR. 




























O Controlador LMl-RMPCR se mostrou adequado no eontrole deste proeesso que 
explora uma eondição limite de estabilidade e de forma menos oseilante que o LMI-RMPC 
elássieo. WANG e RAWLING (2004 b) mostraram que nem todos os eontroladores 
conseguem eontrolar proeessos nesta forma. A Figura 4.11 mostra de forma eomparativa os 
eomportamentos dos controles LMl-RMPC e LMl-RMPCR. 
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Figura 4.11 - Comparação dos controladores LMl-RMPC x LMl-RMPCR. 


4.4 Coluna de Destilação 


Nesse tópieo será eonsiderado o efeito da ineerteza e do eontrole robusto em um 
sistema MIMO mal condieionado. O exemplo esooUiido foi o da coluna de destilação, 
bastante estudada por Skogestad e eo-autores (WANG; RAWLINGS, 2004a; 2004b). Os 
modelos são apresentados nas Eq. 4.11-4.12. Pela variação do elemento (1,1) da matriz B de 
0,1591 para 0,2393, o número de eondieionamento deeresce de 132 para 9,33. 


Xk+l = 
y = x 

Xk+1 = 


0,8187 

0 

Xk + 

"0,2393 

-0,1563 

0 

0,8187_ 

0,1958 

-0,1985 

0,8187 

0 

Xk + 

"0,1591 

-0,1563 

0 

0,8187 

0,1958 

-0,1985 


y = x 


(4.11) 


(4.12) 
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A Figura 4.12 mostra o comportamento em malha aberta da saída y(l). 



Figura 4.12 - Malha aberta: sistema eoluna de destilação mal eondieionado 


O problema servo foi exeeutado para o vetor de setpoint: [1 0,1]. As matrizes de 
pesos adotadas foram iguais às matrizes de ordem apropriada eonforme indieado na Tabela 
4.4. O horizonte de eontrole adotado foi o de 5 passos, a máxima variação das entradas foi 
eonsiderando o valor 10, e o horizonte de eontrole e o horizonte de predição são iguais. Os 
easos a serem analisados são: easo 1- eontrolador ealeulado em G1 e aplieado em G2, e easo 
2- eontrolador ealeulado em G2 e aplieado em Gl. O tempo de amostragem é de 1 s. 


TABELA 4.3: AJUSTES PARA CONTROLADORES. 


Controlador 

0 

R 

s 

Incerteza 

RMPC-MMR 

I 

I 

I 

0,11 

LMl-RMPCR 

I 

0,011 

I 

0,11 


As Figuras 4.13-4.14 mostram o eomportamento nos easos 1 e 2, respeetivamente, 
para o eontrolador RMPC-MMR. 

As Figuras 4.15-4.16 mostram o eomportamento nos easos 1 e 2, respeetivamente. 


para o eontrolador LMl-RMPCR. 
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Figura 4.13 - Comportamento do sistema eom RMPC-MMR: easo 1. 
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Figura 4.14 - Comportamento do sistema eom RMPC-MMR: easo 2. 
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Figura 4.15 - Comportamento do sistema com LMI-RMPCR: caso 1. 


(b) 




(b) 


Figura 4.16 - Comportamento do sistema com LMl-RMPCR: caso 2. 
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As respostas mostradas nas Figuras 4.13-4.16 mostram-se plenamente 
satisfatórias. O controlador LMl-RMPC só funcionou na forma do sistema em Au. As saídas 
não exibiram offset. Considerando-se o algoritmo apresentado no Apêndice A, o mesmo 
apresenta respostas instáveis para qualquer condição de desvio/planta modelo neste caso. As 
Figuras 4.17-4.18 apresentam as respostas comparativas para os controladores RMPC-MMR e 
LMI-RMPCR para os casos 1 e 2, respectivamente. 



Figura 4.17 - Comparação RMPC-MMR x LMI-RMPCR caso 1. 
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Figura 4.18 - Comparação RMPC-MMR x LMI-RMPCR caso 2. 
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4.5 Coluna de Flotação 


Os processos de concentração de minerais são baseados nas diferenças das 
propriedades dos minerais eonstituintes. O proeesso de flotação utiliza eomo prineipal 
propriedade difereneiadora as earacterístieas físico-químicas da superfície das partíeulas dos 
diferentes minerais (SILVA et al., 2004). Através da adição de agentes espeeífícos, aumenta- 
se ou diminui-se o grau de hidrofibicidade de uma partícula mineral. Fazendo-se passar um 
fluxo de bolhas de ar em eondições hidrodinâmieas de baixa turbulência, através da polpa do 
minério, as partículas hidrofóbicas são coletadas ao longo da seção de coneentração e levadas 
à seção de limpeza. As hidrofílieas não coletadas permanecem em meio aquoso. Os agregados 
bolhas-partíeulas transportados para a região superior da eélula de flotação eneontram o fluxo 
de água de lavagem em contracorrente, cuja fínalidade é retornar para a seção de eoncentração 
as partíeulas hidrofílieas arrastadas hidraulie amente, como apresentaram DEL VILLAR, 
GRÉGOIRE e POMERLEAU (1999). A Figura 4.19 mostra um esquema básieo da coluna de 
flotação. 



Figura 4.19 - Esquema de uma coluna de flotação (BERGH e YIANATOS, 2003). 
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Controle Estabílízante: Quando o objetivo da estratégia de eontrole se restringe a estabilizar a 
eoluna de flotação dentro de seus limites operaeionais, pode-se atribuir ao sistema de eontrole 
somente a tarefa de eontrolar a altura da eamada de espuma (SILVA et al, 2004) A prineipal 
vantagem desta estratégia é a simplicidade do esquema de controle. Contudo, para se alcançar 
resultados metalúrgicos razoáveis, é conveniente utilizar a estratégia de controle chamada 
controle otimizante. 

Controle Otimizante: Esta estratégia se baseia em modelos matemáticos ou em regras 
heurísticas que relacionam objetivos secundários do esquema de controle (altura da camada de 
espuma, holdup de ar e hias) com objetivos primários (teor e recuperação) (CARVALHO, 
DURÃO e FERNANDES, 1999; PERSECHINI, JOTA e PERES, 2004). Contudo, a grande 
interação entre as variáveis de processo e a presença de não linearidades agudas na relação 
entre elas, são alguns dos principais fatores responsáveis pela dificuldade no planejamento de 
uma estratégia de controle otimizante. 

A aplicação do movimento de controle foi feita utilizando a forma contínua, e o 
cálculo do controlador foi feito utilizando-se da forma discreta. Um pequeno erro que possa 
ocorrer na discretização, juntamente com uma pequena incerteza adicionada ao sistema em 
que foi aplicado o controle, constituem uma forma de se avaliar o comportamento do sistema 
sob 0 desvio planta/modelo. O aspecto variante no tempo da matriz de entrada utilizada no 
modelo para cálculo do controlador, B(6), tenta capturar os efeitos da incerteza multiplicativa 
composta por ôi, Ô2 e Ô3, quando for 0 caso, uma vez que os sistemas avaliados possuem 
ordens diferentes. A resposta do sistema com 0 controle RMPC-MMR foi satisfatória. O 
horizonte de controle (N^) foi fixado em 5, as matrizes Q e R são matrizes identidade de 
ordem adequada e a perturbação foi da ordem de 10% para cada estado. 

Para 0 caso de aplicação de controle robusto foi utilizado 0 modelo, cuja descrição 


é mostrada a seguir: 
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• Modelo apresentado por Perseehini, Jota e Peres (2000) de uma eoluna de flotação piloto. 
Neste easo, as variáveis de saída são: altura da eamada de espuma (y(l)), holdup de ar na 
seção de limpeza (y(2)), e vazão de hias (y(3)). As variáveis manipuladas são: vazão de 
água de lavagem (u(l)), vazão de ar (u(2)) e vazão de flotado (u(3)). A vazão de flotado é 
eonstituída pela relação entre vazão de alimentação e a vazão de não flotado. A estrutura 
do sistema é apresentada na Eq. 4.13. As funções de transferêneia no domínio de Laplaee 
são apresentadas nas Equações 4.14-4.22. As eonstantes y e (3 valem, respeetivamente, 
-4,8.10"^ e -3,84.10^. O tempo de amostragem é de 10 s. 





G,, 

G." 

■u(l) 

y(2) 

= 


G,. 

G33 

u(2) 

_y( 3 )' 


G3, 

G32 

G 33 _ 

_u(3) 


-(l,029.10'^s +2,3.10"^) 

(s + 4,02.10 ")(s + l,92.10"') 

-l,59.10""s + 4,33.10'' 

(s + 4,02.10"')(s +7,981.10"') 

(1,029.10"'s +2,3.10"') 

(s + 4,02.10"')(s +1,92.10"') 

7,6.10"' 

(s +1,92.10"') 

7,78.10"' 

(s + 7,981.10"') 


-7,6.10"' 

(s +1,92.10"') 


(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 


G3j=s(yG,^+|3G2,) (4.20) 

G33=s(yG,3+|3G33) (4.21) 


G33=l + s(yGj3+pG33) 


( 4 . 22 ) 
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A Figura 4.20 mostra o comportamento das saídas no domínio do tempo do 
sistema para uma perturbação inieial de 10 % nos estados em problema regulador. A Figura 
4.21 mostra a máxima ineerteza ealeulada, as ações de eontrole mínimas ealeuladas, o tempo 
de eáleulo do eontrolador e o eomportamento da função de otimização. As respostas mostram- 
se satisfatórias. 



Figura 4.20 - Comportamento do sistema de flotação eom o eontrole RMPC-MMR. 
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Figura 4.21 - Comportamento do sistema de flotação com o controle RMPC-MMR. 
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4.6 Implementação Computacional 


Os algoritmos de controle MPC, IHMPC, RMPC-MMTV e RMPC-MMR tiveram 
as suas formulações manipuladas de maneira a se obter a forma de otimização quadrática 
possibilitando o uso de pacotes computacionais disponíveis. O método usado para resolver o 
problema de programação quadrática é uma estratégia de conjunto ativo como método de 
projeção (COLEMAN; ZHANG, 2000). Para isso, foi utilizado o algoritmo ^wa^i-Newton, 
composto de duas fases que constituem a definição da direção de busca (usando o método 
BFGS - Broyden, Fletcher, Goldfarb, and Shanno) e a pesquisa em linha. Os algoritmos de 
controle robusto usando FMI foram resolvidos com a utilização do método dos pontos 
interiores de Nesterov e Nemirovski (GAHINET et al, 1995). 

Os programas foram implementados em computador pessoal com as seguintes 
configurações: 256 MB de memória RAM e processador AMD Atlhon 1.1 GHz (máquina 1) e 
processador Core 2 Duo 2.2 GHz com 1 GB de memória RAM (máquina 2). As 
implementações de maior dificuldade, como esperado, foram as que consideram 
explicitamente a incerteza. Contudo, todos os resultados apresentados tiveram suas ações de 
controle calculadas dentro do tempo de amostragem. 
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4.7 Considerações Finais 


Nesse eapítulo foram aplieados os algoritmos de eontrole: IHMPC, RMPC- 
MMTV, RMPC-MMR, LMI-RMPC e LMI-RMPCR. Os diversos sistemas avaliados foram: 
CSTR (sistema de Van der Vusse), sistemas lineares diseretos, aqueeedor, eoluna de 
destilação eom modelo mal eondieionado e sistema de flotação. Esses modelos exploram as 
diversas difieuldades para os sistemas de eontrole reais, eomo por exemplo: aplieação de 
modelo linearizado em proeesso não linear, diserepâneias entre tempo morto, fases, ganhos, 
eonstantes de tempo, tempo de diseretização e, ainda, difieuldades assoeiadas a sistemas 
MIMO, eomo o mal eondieionamento, por exemplo. 

O objetivo desse eapítulo foi avaliar os eontroladores RMPC-MMR e LMI-RMPCR, 
desenvolvidos no Capítulo 3, em diversas eondições que, de alguma forma, exigiram alguma 
earaeterístiea de robustez, de estabilidade ou de desempenho. Foram feitas eomparações eom 
MPC elássieo (IHMPC, Toolbox MPC do Matlab®) e RMPC elássieos (RMPC-MMTV e 
LMI-RMPC). No Capítulo 6 os eontroladores serão aplieados ao sistema experimental 
multivariável de eontrole de nível e pH por manipulação das vazões de entrada, que será 


apresentado no Capitulo 5. 



CAPITULO 5 
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5 O SISTEMA EXPERIMENTAL - DESCRIÇÃO 


5.1 Introdução 


Nesse Capítulo desereve-se em detalhes a montagem experimental, a 
instrumentação, o hardware e o software empregados na eonstrução experimental. 

O sistema experimental eonstruído é um proeesso de neutralização que oeorre em 
um reator tanque agitado. As vazões de áeido e de base são as variáveis manipuladas, a vazão 
de solução tampão é a perturbação não medida ou medida, a depender de eomo a mesma é 
tratada. Para o easo desse estudo, a vazão de solução tampão foi utilizada apenas nos estudos 
de simulação, de forma a se avaliar o eomportamento do sistema em determinadas eondições. 
Na parte experimental não foi utilizada a vazão de solução tampão em nenhum easo. A altura 
do nível no reator (h) e o pH são as variáveis eontroladas ou saídas. O sistema foi 
espeeialmente eseolhido para fins de eontrole por apresentar um forte eomportamento 
dinâmieo não linear e por resultar em um baixo eusto de operação do ponto de vista de 
eonsumo de insumos eomo energia elétriea e matéria-prima (MONTANDON, 2005). O 
sistema se toma interessante do ponto de vista de eontrole por se tratar de um proeesso, 
multivariável e não linear, em que o ganho estátieo experimenta fortes variações dentro da 
faixa operaeional de vazões. A Figura 5.1 mostra a eurva de titulação do sistema. As eurvas 
eonsiderando a adição de solução tampão são todas resultantes de simulação. Nessa Figura, qi 
é a vazão de áeido e q 2 é a vazão de base. Pode-se observar que o pH do proeesso sofre uma 
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mudança de 4 para 9 no intervalo [0,87; 1,22] de qi/qi. É justamente essa mudança no ganho 
estátieo do proeesso que torna o problema de neutralização difíeil do ponto de vista de 
eontrole, pois pequenas modifieações na relação q 2 /qi podem eausar grandes ou pequenas 
variações no pH, dependendo da loealização do ponto eorrente de operação dentro da faixa 
operaeional. Ressalta-se que num proeesso analítieo de titulação apenas uma gota é sufieiente 
para deslocar o pH de 4 para 9. Desta forma, pode-se inferir a sensibilidade do processo 
dinâmico associado e a dificuldade para o seu controle em tempo real. Pode-se notar que a 
adição de solução tampão suaviza o comportamento não linear do sistema, facilitando a ação 
do controlador, uma vez que se o sistema fosse linear, a relação dada pelo gráfico da relação 
de vazões pelo pH seria uma reta, cuja aproximação pode ser percebida pela curva simulada 
de maior quantidade de solução tampão. 



Figura 5.1 - Curva de titulação: medidas experimentais x modelo simulado variando a vazão 
de tampão. 
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5.2 O Sistema Experimental 


5.2,1 O Reator 


O sistema consiste de um reator tanque continuamente alimentado com uma 
corrente de ácido (HNO3) e outra de base (NaOH), impulsionadas por bombas de 
engrenagens. O conteúdo do reator é misturado por um agitador mecânico que, juntamente 
com as chicanas, promovem 0 contato entre as moléculas do ácido e da base. Após a reação, 
os produtos deixam 0 reator por gravidade, sendo coletados em uma tubulação e enviados ao 
sistema de tratamento de resíduos e descarte de líquidos do laboratório. Dessa forma, 
possibilita-se uma interação entre as saídas do processo. O pH da solução do reator é medido 
na tubulação de saída. Montandon (2005) executou a medição do pH dentro do reator, sendo 
percebido 0 aparecimento de um intenso ruído durante a operação da unidade. A adoção do 
novo ponto de medição limitou significativamente esse ruído. Os valores nominais de cada 
variável chave do processo e de cada parâmetro são dados na Tabela 5.1. A Figura 5.2 
apresenta uma foto do sistema experimental. 


TABELA 5.1: VARIÁVEIS CHAVES E PARÂMETROS DO REATOR. 


Variáveis 

Símbolo 

Valores nominais 

Altura 

h 

25 cm 

Área 

Ar 

168,38 cm^ 

Volume 

Vr 

4209,67 cm^ 

Vazão de ácido 

qi 

11,9130 mL/s 

Vazão de base 

q 2 

11,8235 mL/s 

Vazão de Solução Tampão 

q3 

0,01 mL/s 

pH 

pH 

7,0 

Cone. ácido em qi 

[HNO 3 ] 

3,510.10'“ M 

Cone. base em q 2 

[NaOH] 

3,510.10'“ M 
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Figura 5.2 - Foto do sistema de neutralização. 


5,2,2 Equipamentos 


Dois tanques de 200 litros eonfeoeionados em polietileno são usados para 
armazenamento das soluções de áeido e de base para o abasteeimento do proeesso. Esses 
tanques são inieialmente preenehidos eom água de abasteeimento públieo e, posteriormente, 
são adieionadas quantidades eonheeidas de áeido nitrieo P.A. e hidróxido de sódio P.A. de 
forma a se obter as soluções áeida e alealina nas eoneentrações dadas pela Tabela 5.1. Com 
esses volumes de soluções, o sistema experimental pode operar aproximadamente por 4,5 
horas sob as eondições nominais. Durante as operações inieiais foi deteetada a neeessidade de 
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se adaptar um filtro na entrada de água para o sistema devido a uma grande quantidade de 
partíeulas em suspensão, eomo mostra a Figura 5.3. Estas partíeulas, além de sujarem o 
sistema, acarretando erro de medição, danificam os transmissores de vazão. Foi então 
adaptado um fdtro constituído de uma parte de tubo de PVC cujas extremidades continham 
uma tela de aço com malha fina e preenchida com algodão. Este algodão foi trocado 
diariamente para assegurar a eficiência na retenção de partículas. A Figura 5.4 mostra os 
elementos filtrantes utilizados para a operação em duas semanas durante o mês de Maio/2008, 
período em que foram realizados os testes experimentais finais. 



Figura 5.3 - Filtro adaptado no sistema de enchimento de água nos tanques. 



Figura 5.4 - Elementos filtrantes utilizado no período de testes experimentais em Maio/2008. 
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A adição desse aeessório ao sistema, entretanto, aearreta uma lentidão no 
enehimento dos tanques. Uma alternativa foi a adaptação de uma bombona de 1000 L e o 
abasteeimento da mesma em regime eontínuo. A filtragem foi feita em seu abasteeimento e 
durante a transferêneia para os tanques de trabalho. A Figura 5.5 mostra o tanque adaptado, 
loealizado na área externa do laboratório onde se eneontrava a unidade experimental. 



Figura 5.5 - Tanque alternativo para armazenamento de água filtrada. 

Além dos tanques de 200 litros, existem mais dois tanques de polietileno de 100 
litros eada um, para preparação de volumes adieionais de soluções de áeido e base que são 
transferidos por bombas eentrífugas para os tanques operaeionais de 200 litros. Dessa forma, 
pode-se operar o sistema sem o riseo de esvaziamento dos tanques de armazenamento de 
reagentes durante os ensaios. De forma a completar o sistema e possibilitar em outros 
experimentos a avaliação experimental com a solução tampão (NaHCOs), existe um tanque de 
120 L apropriado para o uso dessa solução integrando o sistema, mas que nesse trabalho não 
foi utilizado. Cada um dos tanques de 200 litros e da solução tampão são equipados com uma 
micro-bomba de engrenagem, modelo Cole-Parmer U-73004-02 dotada de um Gear Pump 
Drive modelo Cole-parmer U-07144-02, que é usada para fornecer solução ao processo 
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(Figura 5.6). As vazões de ácido, de base e de solução tampão podem ser reguladas através da 
corrente elétrica (4-20 mA) enviada às bombas de engrenagem pelo sistema de aquisição e 
controle de dados. Esta foi uma opção de controle do sistema, feito diretamente pelas 
correntes das bombas quando houve, durante os testes experimentais, a queima de um dos 
tramsissores das bombas. No entanto, como houve a aquisição de um novo equipamento para 
a conclusão dos testes experimentais, optou-se pelo controle das vazões, em que a vazão entra 
como um valor de referência para um sistema escravo de controle de corrente para as bombas, 
conforme será mais detalhado no Capitulo 6. 



Figura 5.6 - Bomba de engrenagem dotada de Gear Pump Drive. 

Nas tubulações que conectam os tanques ao reator e antes das bombas de 
engrenagem, existem filtros para reter o material particulado e executar a purificação das 
correntes de ácido e base. Esses filtros removem sólidos que podem danificar os medidores de 
vazão e são trocados periodicamente, a cada seis meses. O conteúdo do reator é misturado 
usando um agitador mecânico que permite a manipulação manual da taxa de agitação. O 
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agitador é fixo à parede e foi ajustado de forma a tentar evitar a propagação de ruídos no 
sistema (Figura 5.7), conforme sugeriu Montandon (2005). 



Figura 5.7 - Agitador mecânico do reator. 


5.2,3 Instrumentação 


A montagem experimental é equipada com instrumentos de precisão e resulta em 
um sistema de controle por computador dependente de medidas. O pH do efluente é medido 
com eletrodo de vidro da Leeds & Northrup modelo 7082, conforme Figura 5.8. O eletrodo é 
colocado na tubulação de saída do reator. Este eletrodo é conectado a um amplificador- 
transmissor de sinal (Figura 5.9), que tem uma faixa de 4-20 mA para o sinal transmitido. 

As vazões de ácido, de base e de solução tampão são medidas usando sensores de 
vazão de infravermelho tipo turbina modelo Cole-Parmer U-32250-00 (Figura 5.10). Cada 
sensor de vazão é conectado a um transmissor-conversor pulso-corrente modelo Cole-Parmer 
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U-94788-11 (Figura 5.11). Esses conversores transformam um sinal de freqüência em um 
sinal de corrente na faixa de 4-20 mA. 



Figura 5.8 - Sensor de pH. 



‘ ^ 1 




Figura 5.9 - Foto do transmissor de pH. 
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Figura 5.10 - Sensor de vazão de infravermelho tipo turbina. 



Figura 5.11 - Transmissor-eonversor pulso-eorrente. 
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Figura 5.12 - Sensor-transmissor de nivel. 

O nível do reator é medido por um sensor-transmissor de pressão difereneial, eujo sinal 
apresenta uma relação linear eom a altura de líquido, mostrado na Figura 5.12. O equipamento 
apresenta dois orifíeios de saída, sendo que um deles está ligado à eoluna de líquido dentro do 
reator de neutralização e o outro está aberto para a atmosfera. A diferença de pressão que 
surge entre estes dois pontos provoea o desloeamento de uma membrana loealizada dentro do 
equipamento que, por sua vez, origina um sinal em voltagem que eorresponde á altura de 


líquido dentro do reator. 
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5,2,4 O Sistema de Aquisição de Dados 


O sistema de aquisição de dados é o SCXI, fabrieado pela National Instruments. 


Este sistema possui uma arquitetura flexível e eonsiste em um eondieionador de sinais de 
entrada e saída eoneetado a uma placa de aquisição (DAQ). Esta, por sua vez, está conectada a 
um computador PC/AT Pentium IV 2.4 MHz através de um slot ISA padrão. A Figura 5.13 
mostra esquematicamente o chassi SCXI-IOOO usado neste trabalho. A Figura 5.14 mostra 
uma representação do sistema de aquisição e da estação de trabalho utilizada. Esta 
configuração ilustra a entrada e saída de sinais condicionados para ambientes DAQ (aquisição 
de dados). O chassi abriga uma variedade de condicionadores de sinais e módulo de 
instrumentação. Sinais analógicos são condicionados e transmitidos para a placa de aquisição 
diretamente ligada ao computador. A placa ainda controla e monitora a entrada/saída de sinais 
digitais e sinais analógicos. 


Placa de 



Módulo 
condicionador 
de sinais 


Figura 5.13 - Representação esquemática do chassi SCXI-1000 e a placa de aquisição de 
dados. 
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Figura 5.14 - Representação esquemátiea do sistema experimental e da aquisição de dados. 


5.3 Softwares Utilizados 


O Labview® é baseado na linguagem de programação G e é responsável pelo 
eontrole em tempo real do sistema. É uma programação tal eomo C ou BASIC que possui uma 
biblioteea de funções para qualquer tarefa de programação. O software inelui ainda biblioteeas 
para aquisição, análise e armazenamento de dados, além de ferramentas eonveneionais eom 
uma exeeução animada por onde se vêem os dados sendo repassados através do programa. 
Difereneia-se das outras linguagens por ser uma linguagem gráfiea e não baseada em textos. 
Possui oaraeterístioas gráfieas iterativas, sendo uma poderosa linguagem de programação 
gráfiea. A versão do Labview® utilizada foi a 8.2. 
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O Matlab® versão 7.5 foi utilizado na etapa de solução on-line do problema de 
otimização assoeiado ao problema de eontrole preditivo. A rotina “quadprog” do Matlab® foi 
utilizada para resolver um problema de programação quadrátiea padrão eom restrições na 
variável manipulada. O método usado para resolver o problema de programação quadrátiea é 
uma estratégia de conjunto ativo conhecido como método da projeção (Gill et al, 1984 e Gill 
et al, 1991). O modelo dinâmico do processo utilizado nessa otimização foi o modelo 
fenomenológico, descrito no Capítulo 6. Foi utilizado também um pacote chamado 
“YALMIP” (LÔFBERG, 2004) para resolução de problemas na forma de LMFs. Esse pacote 
implementa uma interface amigável e utiliza o método dos pontos interiores para resolver o 
problema de otimização. O método dos pontos interiores e os algoritmos para solução de LMI 
foram os do Toolbox LMITool®. Foi também utilizado o Toolbox MPC® para programação de 
controladores comparativos e cujo algoritmo encontra-se descrito no Apêndice. A Figura 5.15 
mostra o sistema de aquisição e a estação de trabalho usada para executar o controle em tempo 


real do sistema experimental. 
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Figura 5.15 - Sistema de aquisição e do microcomputador utilizado no trabalho. 


5.4 Considerações Finais 


Nesse Capítulo o sistema experimental de controle de nível e pH utilizando as 
vazões de ácido e base foi detalhadamente descrito. Os componentes de hardware e software 
do sistema foram apresentados. No Capítulo 6 será mostrada a operação em malha aberta do 
sistema e no Capitulo 7 os resultados de simulação e experimentais. 
















CAPITULO 6 
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6 MODELAGEM DO SISTEMA EXPERIMENTAL 


6.1 Modelo Fenomenológico 


O modelo físico do processo, apresentado na Figura 6.1, foi desenvolvido por Flall 
(1987) e é baseado nas hipóteses de mistma perfeita, densidade constante e solubilidade 
completa dos íons envolvidos. As reações químicas envolvidas no processo de neutralização 
ácido-base (HN 03 -Na 0 H) são mostradas nas Eq. 6.1-6.5. De forma a se ter um modelo 
completo, será considerada a presença de solução tampão (NaHC 03 ). Contudo, para os 
resultados experimentais, não foi utilizada essa solução. 



Figura 6.1 - Diagrama esquemático do reator de neutralização em malha aberta. 
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HNO 3 ^ H+ + NO 3 

( 6 . 1 ) 

NaOH ^ Na+ + OH“ 

( 6 . 2 ) 

kw 


(6.3) 

kal 

NaHC03+H20^Na0H + H2C03 , H 2 CO 3 + HC03‘ 

(6.4) 

ka2 

HCO 3 +C03^' 

(6.5) 


As quantidades Wa e Wb são definidas para o equaeionamento deste problema e são 
ehamadas de invariantes porque suas eoneentrações não são afetadas ao longo da reação. As 
reações são rápidas o bastante para permitir que o sistema seja eonsiderado em equilibrio. 
Logo, as reações de equilíbrio podem ser usadas para determinar a eoneentração de íons de 
hidrogênio através da eoneentração dos invariantes reaeionais. As eoneentrações de equilíbrio 
são dadas pelas Eq. 6 . 6 - 6 . 8 . 


kal = 


HCOj|H\ 

[H 2 CO 3 


jco^d 

HCOj 


kw = [H^J[0H 


( 6 . 6 ) 


(6.7) 

( 6 . 8 ) 


Segundo Gustaffson e Waller (1983), dois invariantes da reação são definidos para a 
i-ésima eorrente: 


w. = [h • 1 - [oh - 1 - [hCO d - 2 [cO f 1 (6.9) 

W„ =[HjC0d + [HC0 d-dcO^d, 


( 6 . 10 ) 
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Combinando as Eq. 6.6-6.10, a seguinte relação implíeita entre 
eneontrada na Eq. 6.11: 



, Wa e Wb pode ser 


W 


a 



K. 


K 


w. 




-Wb 


ai 

+ 2 

H+ 

H+ 


1 + 


K 

H 


al . ^al^a2 

^ .12 

H 


( 6 . 11 ) 


A Equação 6.12 é o balanço de massa (m) no reator. Uma vez que a densidade é eonstante por 
hipótese, ehega-se ao sistema de equações difereneiais expresso nas Eq. 6.13-6.15. 


dm .... 

-= m^ -I- m2 -I- m3 - 


dh 1 

- = -(q.+q,+q,-q.) 

^ ^ (q, (w., - w. )+ q, (w., - W. )+ q, (W., - W.)) 


^ (q, (W„ - W J + q, (W„ - W J + q, (W„ - W J) 

dt Vr 


( 6 . 12 ) 

(6.13) 

(6.14) 

(6.15) 


2 3 

Ar é a área do reator em em ; Vr é o volume reaeional em em ; h é a altura da solução no 
reator dado em em; qi é a vazão de áeido; q 2 é a vazão de base; qs é a vazão de solução 
tampão, e q 4 é a vazão de saída do reator. Todas a vazões são expressas em mL/s. Wa é o 
invariante reaeional do áeido e Wb é o invariante reaeional da base, que são expressos eomo 
funções da eoneentração de íons de hidrogênio. As Equações 6.12-6.15 permitem ealeular a 
altura e os invariantes reaeionais eomo uma função do tempo, admitindo-se que sejam 
eonheeidas as eondições inieiais do sistema. O pH da solução pode então ser determinado a 
partir da Eq. 6.11 e dos valores das eoneentrações dos invariantes reaeionais. As eondições 
nominais de operação foram apresentadas na Tabela 5.1. 
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A vazão de saída do sistema é por gravidade. Conforme mostra a Figura 6.2, uma 
das válvulas, que nesse easo é do tipo válvula globo, é ajustada de forma a estabeleeer a vazão 
do sistema e eujos parâmetros (Ov e p7) devem ser estimados para a definição eompleta do 
modelo do sistema eonforme expresso na Eq. 6.16. Essa válvula é indieada em destaque na 
Figura 6.2. Após a eorrida/estimação de parâmetros foi retirada a alavanea de movimentação 
da válvula, eomo pode ser observado na Figura 6.2, a fim de assegurar a manutenção dos 
parâmetros estimados durante o eonjunto de eorridas. Caso oeorressem alterações nas 
respostas, ou modifieações indevidas da posição da válvula, seria neeessária a exeeução de 
nova estimação de parâmetros. 

q4=eX^ (6.16) 



Figura 6.2 - Válvulas de saída do sistema de neutralização para estimação de parâmetros. 

Para estimação desses parâmetros e validação do modelo, faz-se neeessária a 
operação da planta através de entradas que exeitem todos os modos dinâmieos e em ampla 
faixa de freqüêneia, de forma que as saídas eontenham informações sufieientes para tal 
proeedimento. As exeitações mais usuais são as do tipo: rampa, degrau e impulso. Contudo, 
para se obter informações em uma ampla faixa de freqüêneia, aplieam-se entradas aleatórias 
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formadas por seqüências de degraus eom amplitude e duração variáveis (ML-PRBS), 
eonforme indiea Montandon (2005). O tempo de amostragem para esse proeesso foi de 10 s. 
As vazões de áeido e de base variaram eonforme deserito nas Eq. 6.17-6.18, eom uma 
probabilidade de passo igual a 0,8, em que a duração de eada degrau foi de 120 s. As vazões 
sobre as quais ineidiram as variações foram de 12 mL/s em eada uma delas e de forma a se 
manter a soma das vazões dentro do limite indieado na Eq. 6.19. As vazões dadas pelas Eq. 
6.17-6.18 foram determinadas por simulação e avaliadas experimentalmente em testes de 
operação eom a unidade experimental. Os valores apresentados abaixo foram definidos de 
forma que durante a operação o reator não esvazie e nem transborde, resultando nos limites 
impostos pela Eq. 6.19. Considerando a parte eonstrutiva da unidade, o transbordamento de 
solução aearreta em danifieação da instrumentação, sujeita nesse easo ao eontato eom o 
líquido. 


q^e [9,856 13,310] 

(6.17) 

q^ e [9,977 13,215] 

(6.18) 

qj-Fq2e[21,5 24,5] 

(6.19) 


A partir do eomportamento dinâmieo do sistema frente às alterações aleatórias nas 
vazões de entrada, é possível avaliar se o modelo fenomenológieo é eapaz de predizer o 
eomportamento real da planta. Para o modelo utilizado na predição das saídas reais, deseja-se 
que 0 erro de modelagem seja aproximado por uma distribuição normal eom média zero e 
variâneia finita, indieando que erro é aleatório e não eorrelaeionado. A Figura 6.3 mostra a 
exeitação para estimação de parâmetros. A Figura 6.4 mostra os resultados experimentais para 
a vazão de áeido. São mostrados os valores de referêneia, ealeulados pelo algoritmo de 
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geração de números aleatórios, e o resultados experimentais. A Figura 6.5 mostra uma 
ampliação nos instantes inieiais da eorrida. 



Figura 6.3 - Exeitação das vazões de áeido e de base para estimação de parâmetros da válvula. 



Figura 6.4 - Exeitação da vazão de áeido: valor de referêneia x valor experimental. 

Pode-se notar uma pequena diferença entre os valores de referêneia e experimentais 
que são resultantes de ruídos na entrada, oriundos do proeesso, da rede elétriea e da 
instrumentação. Existe também uma dinâmiea na vazão experimental dada pela ação de 
eontrole do tipo PI, pois quando a exeitação é feita no proeesso, na realidade o que oeorre é a 
alteração de setpoint para a malha de eontrole de vazão. 
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Figura 6.5 - Excitação da vazão de ácido: valor de referência x valor experimental - 
Ampliação de escala para os instantes inieiais da eorrida. 

A finalidade dessa implementação é manter a vazão eontrolada em parâmetros desejados, 
prineipalmente eonsiderando-se que na medida em que os tanques de armazenamento têm o 
seu nivel alterado, eonforme a operação, podem oeorrer variações na vazão, oeasionando 
alguma interferêneia no proeesso em eontrole. Cabe ressaltar ainda que todas as medidas 
experimentais foram feitas usando filtros digitais, eonforme indieam as Figuras A.5-A.6 e 
Gonçalves (2007). 

A Figura 6.6 mostra os resultados experimentais para a vazão de base, 
eonsiderando-se as entradas aleatórias previamente ealeuladas. 



Figura 6.6 - Exeitação da vazão de base: valor de referêneia x valor experimental. 
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Para a estimação dos parâmetros da Eq. 6.16 utilizou-se um eonjunto de dados experimental 
eomposto por 6000 pontos de aquisição. Para a estimação foi utilizada metade dos dados e, 
para a validação, os dados restantes. Esse proeedimento é feito de forma a assegurar que o 
modelo tenha eapaeidade de predição de dados deseonheeidos. O proeedimento é uma solução 
não linear de aproximação a uma eurva, utilizando o método dos mínimos quadrados. A 
estimativa inieial foi dada por um valor aleatório. Os valores estimados resultantes foram: 
e^ =20,4477 e p7 = 0,0523 . A Figura 6.7 mostra o resultado da etapa de estimação, que foi 

feita eonsiderando-se uma predição eom horizonte de passo igual a 15, isto é, a eada 15 passos 
é feita a atualização eom o dado experimental. Este valor de passo foi adotado já que se 
esperava trabalhar aproximadamente nessa faixa de predição do modelo no horizonte do 
eontrolador. A Figura 6.8 mostra os resultados para a etapa de validação. A média do erro 
eometido na etapa de validação foi de - 0,0084 e eom um desvio padrão de 0,117. A Figura 
6.9 mostra o histograma do erro eometido na predição da altura. Pode-se notar que o gráfieo 
se aproxima razoavelmente de uma distribuição gaussiana de média zero, indieando que a 
estimação não é tendeneiosa e por isso válida para os propósitos de eontrole pretendidos. 

Estimação 



Figura 6.7 - Respostas experimental e teóriea- Etapa de estimação. 
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Validação 



Figura 6.8 - Respostas experimental e teóriea - Etapa de validação. 



Figura 6.9 - Erro na predição da altura. 


Foi então realizada outra eorrida em malha aberta. Dessa vez, foram feitas 
alterações manuais e aleatórias, tanto em valor eomo em duração, variando o nivel do reator 
entre 5 e 23 em. As Figuras 6.10-6.11 mostram os valores experimentais para as vazões de 
áeido e de base. A Figura 6.12 apresenta os dados experimentais eomparados eom a predição 
da altura eom passo igual a 1, isto é, a eada passo o valor do modelo é atualizado eom o valor 
experimental, e a Figura 6.13 é o resultado para predição da altura eom passo infinito. 





























Vazão de Base (mL/s) o Vazão de Ácido (mL/s) 
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- Vazões de áeido para eorrida em malha aberta. 



Figura 6.11 - Vazões de base para eorrida em malha aberta. 



Figura 6.12 - Dados experimentais de eorrida em malha aberta x Predição de altura usando 
modelo fenomenológieo eom passo 1. 
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Figura 6.13 - Dados experimentais de eorrida em malha aberta x Predição de altura usando 
modelo fenomenológieo eom passo infinito. 

Para esse easo apenas o valor experimental inieial é utilizado fazendo-se a integração do 
modelo eom as entradas. O distaneiamento observado entre os valores teórieos e 
experimentais deve-se ao fato de que o sistema é muito próximo de um sistema integrante. Tal 
fato pôde ser eomprovado na operação da unidade, uma vez que os valores de vazão 
simulados (12 mL/s) indieam a estabilização do nivel e, na realidade, aeonteee um enchimento 
ou esvaziamento continuo do reator. Desta forma, diferenças entre os valores usados para a 
simulação e o valore real se configuram em desvios que são adicionados e integrados no 
modelo resultando na diferença entre as curvas da Figura 6.13. O modelo fornece predições 
qualitativas da dinâmica do sistema, mesmo para o caso de predição com passo infinito. As 
informações que se afastam do valor real podem ser deslocadas para a incerteza, levando a 
exigir dos controladores sua característica de robustez. 

A Figura 6.14 mostra uma corrida em malha aberta e os resultados para a predição 
do modelo. Como pode ser observado, o modelo atende aos propósitos de controle predizendo 
de forma muito próxima o nível, e de forma razoável o pH. Como esperado, é demonstrada 
uma certa dificuldade do modelo diferencial não-linear predizer o pH na faixa entre os valores 
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de 6 e de 8 que, eonforme indiea a Figura 5.1 (eurva de titulação), trata-se da região de alta 
não linearidade. 



Figura 6.14 - Malha aberta: resultados teórieos e experimentais. 
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6.2 Modelo na Forma de Espaço de Estados 


A formulação em espaço de estados resultante da manipulação algébriea é 
mostrada nas Eq.6.20-6.27. Os termos dWadz e dWbdz são resultantes da linearização da Eq. 
5.11 e são apresentados nas Eq. 6.28-6.29. A Eq. 6.30 mostra a relação de dydz. 


x = [h W, Wb] 

(6.20) 

e 

II 

(6.21) 

y = [h pH] 

(6.22) 

d = [93 ] 

(6.23) 


A = 
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dydz = -logio — 

Irl 


(6.30) 


6.3 Considerações Finais 


Nesse eapítulo foi apresentado o modelo difereneial não linear do sistema 
experimental. Foi também obtido o modelo na forma de espaço de estados, utilizado nos 
cáleulos dos eontroladores preditivos. Foram desenvolvidas as rotinas no Labview®, 
apresentadas no Apêndiee A, para permitir a operação da planta experimental em tempo real. 
Realizou-se a operação em malha aberta do sistema para estimação/validação de parâmetros 
neeessários e do modelo. Os resultados mostraram que as estimativas realizadas levaram a um 
modelo eonsistente que atende aos requisitos para a aplieação de eontroladores. 



CAPITULO 7 
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7 MPC E RMPC APLICADOS AO SISTEMA EXPERIMENTAL: 
SIMULAÇÕES, RESULTADOS E DISCUSSÕES 


7.1 Introdução 


Nesse Capítulo serão apresentados os resultados de simulação eom os seguintes 
eontroladores: IHMPC, MPC-RS, RMPC-MMR, LMI-RMPCR (forma nominal) e LMI- 
RMPCR (ineerteza na forma de politopo). Os eontroladores robustos (RMPC-MMR e LMI- 
RMPCR) terão suas respostas eomparadas eom os eontroladores elássieos. A eomparação de 
maior relevâneia nesse easo é a que foi realizada eom o IHMPC. Os eontroladores IHMPC e 
os eontroladores RMPC desenvolvidos nesse trabalho são eontroladores da elasse de eontrole 
de horizonte infinito. 

Para a implementação em tempo real foi neeessária a manipulação dos algoritmos 
para atender a forma da linearização loeal, eomo apresentado na Seção 7.2. O trabalho eom 
algoritmos sueessivamente linearizado no ponto de operação leva à forma de sistemas LTV 
{Linear Time Variant), eujas formas eneontram-se indieadas na Seção 7.3. O Capítulo 
também trás um breve estudo sobre observadores de estado, os resultados de operação em 
malha feehada e alguns dos prineipais problemas que oeorreram durante a operação da planta 
experimental. 

Cabe ressaltar que os resultados experimentais referem-se apenas aos eontroladores 


elássieos: IHMPC e MPC-RS. 
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7.2 Modificações nos Controladores em Função da Linearização Sucessiva 

Considerando que a linearização do sistema experimental apresentado nas Eq. 7.1-7.2 
resulta nas Eq. 7.3-7.4, euja forma difere da formulação olássiea para espaço de estados, torna- 
se neeessário modifiear os algoritmos para atender a eondição de linearização. Isto devido aos 
termos fk i e Ok i. Caso a linearização fosse feita na eondição de estado estaeionário, chegar- 
se-ia à forma clássica, uma vez que estes termos seriam iguais a zero. Contudo, a linearização 
no estado estacionário não atenderia á exigência do processo, pois devido á sensibilidade, um 
controlador projetado baseado neste ponto não conseguiria controlar o sistema, já que para 
algumas regiões o estado estacionário poderia encontrar-se muito distante do ponto desejado, 
mesmo com linearização sucessiva, tendo assim um modelo não representativo para controle 
do sistema, principalmente nas regiões de alto ganho. As modificações para esses 
controladores IHMPC, RMPC-MMR e LMI-RMPC seguem as idéias mostradas nessa seção. 
O controle MPC-RS já considera os termos resultantes da linearização em sua formulação, 
conforme pode ser observado no Capítulo 2. 

^ = f(x,u) (7.1) 

dt 

y = h(x) (7.2) 

dx' 

^ = + + (7.3) 




(7.4) 
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Existem diversas formas de tratar o problema. Uma delas é fazer a expansão dos 
estados eonforme a Eq. 7.5, eriando uma nova matriz de estados e assim pode-se manter a 
estrutura original do problema de otimização, eonforme Capítulo 3. O termo relaeionado eom 
a linearização das saídas não aearreta problemas para o eontrolador, uma vez que apenas faz o 
mapeamento estado-saída e não interfere na dinâmiea do sistema eomo o termo fk-i- Feita essa 
modifieação, é possível trabalhar eom o eontrolador definido anteriormente. Desta forma, 
eonsiderando a faeilidade que esta formulação ofereee, a mesma foi adotada para a operação 
em tempo real. Outra forma pode ser obtida em Reis (2008). 


dx' 

A,, 

0 “ 

/V / 

X 

dt 

0 


1 


+ B, jU 


7.5) 


7.3 Controladores Preditivos na Forma LTV 


No Capítulo 2 foram apresentados os eontroladores preditivos e robustos olássieos 
utilizados eomo base para esse trabalho. No Capítulo 3 foram desenvolvidos dois algoritmos 
de eontrole robusto. Em todos os easos eitados, eonsiderou-se que o sistema era invariante no 
tempo {Linear Time Invariant - LTI). Para a aplieação real no sistema experimental, que 
possui oaraeterístioas não lineares, faz-se neeessária a linerização sueessiva do modelo do 
proeesso de forma a se apliear os algoritmos lineares. A evolução temporal do proeesso 
resulta, então, em um sistema variante no tempo eujos algoritmos são deseritos a seguir. 

O IHMPC eom matrizes variantes no tempo é definido de forma que o movimento 
de eontrole Uk, no instante de amostragem k, é dado pelo primeiro elemento Uk|k de 
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S:{Uk|k, Uk+i|kv • •, Uk+Nm-i|k} na forma da Eq. 7.6. Esta equação pode ser manipulada 
algebrieamente resultando na Eq. 7.7, que transforma o horizonte de eontrole inifinito para 
horizonte de controle finito. A seqüência de ações de controle é a solução do problema de 
otimização dado pelas Eq.7.8-7.10, em que: Q = > 0, R = > 0, S = > 0, 

Qk=Ck^QC,+A/Q,A, (Eq. de Lyapunov), e Auk+j Wk+j" Wk+j-i. No instante k 3“ Nj^^, e 

posteriormente, o vetor de entradas u é um vetor de zeros. O vetor z representa o vetor de 
estados reais medidos ou estimados das saídas medidas. 

Jk = 

J, =mmS^H,S + 2S^(G,z,-F,u,_j) 

Com as seguintes matrizes: 

b/Q^b,+r + 2S b/a/q,b,-s b/a/" q,b, 

B/QkA,B,-S B/Q,B,+R + 2S BJa/" Q,B, 

B/QkAk""'B, B/Q,A/"^B, B/Q,B,+R + 2S 

' B/q,A, ■ 

Bk^QkA,^ 

s~ 

0 

0 



(7.8) 


(7.9) 


(7.10) 


(7.6) 

(7.7) 
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Para o caso do controlador MPC-RS, a função objetivo do problema de 
otimização eonsiderando o sistema LTV {Linear Time Variant) é mostrada nas Eq. 7.11-7.15. 





b 
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K,s(yf-yk)-c,A,Ax,-^^j^ 

K,sAyf..-(C,A,+K,,C,)0,Ax, 

K,sAyf.2-(CkA,+K,3C,)0,^Ax, 

K,sAyL„ -(C^A, +K,3C,)0,""-Ax, 
K,sAyL„.i -(C,A, +K,3C,)0,^"-^'Ax, 
K.sAyr.N„.2 -(C^A, +K,3C,)0,""-^^Ax, 

K,,Ayf,,^ -(C,A, +K,3 C,)^,'^^Ax, 
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O RMPC-MMR na forma de sistema variante no tempo é definido de forma que o 
movimento de eontrole u(k), no instante k de amostragem, é dado pelo primeiro elemento Uk|k 
de S:{uk|k, Uk+i|k,- • •, Uk+Nm-i|k}, que é a solução do problema de otimização dado pelas Eq. 
7.16-7.34. 


mm max 

Hea” ©ee” 


JdHi,eJ=x„'nx„+3jAs. 


(7.16) 


Sujeito a: 


0k^i0k 


Jk(Sk,0k)^Jk(s:,0kXV0,^0 


(7.17) 

(7.18) 


«k+q^ U,j=0 ,... 

,N^-1 



(7.19) 

Ujj+q+j = 0, quando 1>0 



(7.20) 

^k+i = +®k(®k)“k 



(7.21) 

Cujas definições neeessárias eneontram-se a seguir: 
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0 0 

• OJ 



^k+1 ~ Í*k+1 




(7.23) 

'^k+1 ~ i^k 

"í.N.-.l 



(7.24) 

®k+l = [0k ■■■ 

8i+n„-,F 



(7.25) 

0 k=[ôi ••• Ô 
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dim(u) J 



(7.26) 



140 


A 


k 


A 


k 


B,(0J = 


B,(0J 


0 

BkOJ 


A A/"X(0k) - B,(0J 


B,(0J = B, 


l + ô‘ 0 
0 l + ô' 


l+ô 


dim{u) 


fi = diag 


Q...Q Q, 


Qk -Q +Ak QkAk 

A = diag 


R...R 

N„ 


u"*” =^..u 


X‘ 


X...X 


(7.27) 

(7.28) 

(7.29) 

(7.30) 

(7.31) 

(7.32) 

(7.33) 

(7.34) 


O LMI-RMPC na forma de sistema variante no tempo é definido de forma que o 
movimento de eontrole u(k), no instante k de amostragem, é dado pelo primeiro elemento Uk|k 
de S:{uk|k, Uk+i|k,. •Uk+Nm-i|k}, dado pelas Eq. 7.35-7.45. 
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(7.35) 


Sujeito à: 
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7.4 Observador de Estado 


Uma vez que o estado pode ser definido apenas eomo uma variável matemátiea, 
neeessária à resolução do sistema, ela pode não ter representação fisiea. Desta forma, é preeiso 
estimar os estados de modo a informar ao sistema um valor de razoável representação, sob 
pena de levar o eontrolador a ealeular uma ação de eontrole para valores inieiais irreais, 
podendo assim apresentar eomportamentos diferentes das respostas desejadas. 

No sistema experimental, definiu-se eomo estados h, Wa e Wb. Para esse easo, os 
dois últimos estados não são medidos, eomo foi apresentado na seção de modelagem do 
sistema, sendo então neeessário estimá-los. Considerando o sistema sem o tampão, o valor de 
Wa pode ser ealeulado pela Eq. 6.9 e Wb é igual a zero, eonforme Eq. 6.10. Apesar de não 
trabalhar na parte experimental eom o tampão, para fins de simulação, faz-se neeessária a 
eonsideração de uma fração mínima de solução tampão. Isto oeorre devido ao fato de a água 
de abasteeimento não ser exelusivamente H 2 O, eontendo íons dissolvidos que levam 0 sistema 
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real a apresentar eomportamentos de uma solução levemente tamponada (MONTANDON, 
2005). Para tal, foi utilizado o filtro de Kalman (KAILATH, 1980), eomo apresenta a Eq. 
7.46: 

Xkuik = + Bu, + ) (7.46) 

Para esse easo, a matriz de transmissão direta é igual a zero e Lk é a matriz de ganho de 
Kalman, eneontrada pela solução da Equação Disereta de Rieeati. ymedk ^ ® vetor das saídas 

medidas. Os estados são estimados a partir das saídas medidas e das entradas enviadas para a 
planta. Os valores para Qobs e Robs foram iguais 0,1.1 de ordem apropriada. 

Os resultados apresentados nessa Seção foram obtidos exelusivamente por simulação 
eomputaeional. Para a obtenção das respostas a seguir foram feitas as eonsiderações: 

• Plant: modelo difereneial não linear; 

• Ruído: modelo difereneial não linear eom adição de ruído braneo, gerado por seqüêneia de 
números aleatórios, na ordem de 20% e y(l) e 10% em y(2). As respostas desse sistema 
devem aeompanhar as respostas de Plant, exeeto pela diferença imposta pelo ruído 
adieionado; 

• Model: modelo linearizado loealmente; 

• Estim: modelo linearizado loealmente eom estimação de estados. As respostas desse 
sistema aeompanham a resposta de Model, euja linearização é utilizada para o eáleulo das 
respostas, para o easo da saída y(2), já que y(l) = x(l) é uma variável medida. 

A Figura 6.1 apresenta uma simulação em que se eonsidera que a planta e o modelo 
possuem um poueo de solução tampão(q 3 = 0,05 mE/s), uma vez que eomo a água de 
abasteeimento naturalmente já possui alguns íons gerando um efeito levemente tamponante 
está é uma eonsideração bastante razoável. As entradas neste easo foram dadas por uma 
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seqüência aleatória de distribuição uniforme. Percebe-se que os resultados para o nivel (y(l)) 
são bastante próximos. Já para o caso do pH (y(2)), percebe-se que o ruído mostra um efeito 
mais visível que para o nível. Pode-se observar uma degradação na resposta do modelo com o 
estimador de estados na região próxima ao pH igual a 7. 




(b) 



Tempo (min) 


(c) 


Figura 7.1 - Simulação para observador de estado: caso modelo perfeito. 
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A planta se tornaria o mais não linear em q 3 = 0 mL, eonforme indiea a Figura 5.1 
(eurva de titulação). A Figura 7.2 mostra o eomportamento do sistema em uma situação de 
desvio planta-modelo. As entradas para esta simulação são também dadas por seqüêneias 
aleatórias de distribuição uniforme (Figura 7.2 e). Nesse easo, a planta inieialmente possui qs 
= 0,05 mL/s e após 160 minutos de operação este valor é zerado, levando o modelo ehamado 
de Plant a apresentar o eomportamento de maior não linearidade. Para Model, Ruído e Estim 
q 3 eontinua a eonsideração de que a vazão de solução tampão se mantenha no valor de 0,05 
mL/s. A partir desse instante, pode-se observar que o modelo e a Plant, na integração dos 
sistemas, apresentam valores diferentes. Observa-se também um aumento na degradação da 
resposta do sistema eom estimador de estado, mesmo em regiões distantes do ponto erítieo do 
sistema representado pelo pH igual a 7. Essa situação de desvio planta modelo pode oeorrer, 
por exemplo, ao repor solução no tanque, pela agitação da solução pronta, uma vez que na 
operação em tempo real foi observada variação eonsiderável da medida de pH nessas 
oeasiões. A situação eontrária é mostrada na Figura 7.3. Para esse easo, o modelo toma-se o 
mais não linear (Model eom q3 = 0 mL/s após 160 minutos de simulação) e a planta se 
manteria tamponada (Plant eom q3 = 0,05 mL/s durante toda a simulação) mostra a situação 
eontrária, e a planta, nesse easo, tem o eomportamento não linear suavizado pela presença de 
tampão. Para a resposta de y(l), ou seja, o nível (Figura 7.3 a), a divergêneia refere-se à 
integração dos modelos, e reflete a diferença na eonsideração de q 3 nos modelos. Pode-se 
notar que o estimador de estado funeiona bem, pois ele eorrige essa diferença. Já para a Figura 
7.3 b, pereebe-se que os valores de pH pelo Model e pela Plant diferem, e, além disso, houve 
degradações nas respostas para o Estim mesmo após a retirada de q 3 em Model, o sistema não 
eonseguiu ehegar á estimação sem offset. 



Entradas (mL/s) pH(adm) Nível (cm) 
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Figura 12 - Simulação para observador de estado - qa = 0 apóslóO min no modelo Plant. 
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(b) 



Figura 7.3 - Simulação para observador de estado - qa = 0 após 160 min em Model 
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A Figura 7.4 apresenta uma simulação eom o eontrolador MPC-RS, eonsiderando 
a situação de desvio planta modelo deserita para a Figura 7.2. É interessante ressaltar que easo 
a planta esteja próxima da máxima não linearidade teóriea (q3 = 0 mL/s), existem problemas 
na região próxima ao pH igual a 7. A, nesse easo, são as variáveis de folga (eonforme Eq. 
7.12) para viabilidade do problema de otimização assoeiado á formulação do eontrolador 
MPC-RS. Isto impliea que o valor desejável para essa variável é zero quando o sistema 
eneontra-se próximo ao estado estaeionário e é diferente de zero quando o sistema é 
submetido á pertubação ou eomando. A simulação indiea que a eorrida não teria problemas na 
estimação de estados. Na operação real da unidade, eontudo, oeorreram problemas na 
estimação de estado, finalizando uma eorrida abruptamente, uma vez que não houve 
eonvergêneia da solução do observador e assim a ação de eontrole não foi ealeulada. A Figura 
7.5 mostra a eorrida eom o eontrolador MPC RS em que oeorreu esse fato. Após análise, 
através de simulação e testes eom operação em tempo real, a situação foi eontornada para 
outras eorridas, utilizando-se os estados do modelo sem o termo de inovação, sendo 
satisfatório para a aplieação dos eontroladores. Nesse easo, o valor para o estado equivale á 
integração loeal do modelo não linear, valor este que deve eoineidir eom o valor para a 
linearização loeal. Nesse easo a suposição é de que o sistema, sob algumas eondições, passe 
por regiões onde a planta é não observável, resultando assim em problemas na estimação de 
estados. A utilização de modelos entrada-saída, eomo por exemplo, as redes neuronais, 
dispensa a utilização de observadores de estado, a depender da formulação do eontrolador 


(MONTANDON, 2005). 
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Tempo (min) 


(a) 




Figura 7.4 - Simulação: MPC-RS com observador de estado. 
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(b) 



Tempo (min) 


Figura 7.5 - Corrida em malha fechada - MPC-RS: problemas com observador de estado. 
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7.5 Operação em Malha Fechada 

Esse tópico trata da operação em malha fechada do reator de neutralização. O nível 
do reator e o pH são as variáveis controladas ou saídas do sistema. As vazões de ácido e de 
base são as variáveis manipuladas ou entradas. A vazão de solução tampão foi modelada 
como uma perturbação, que pode ser tratada como perturbação medida ou não medida, a 
depender de como ela é utilizada no sistema/controlador e nos modelos de simulação. Para a 
operação em malha fechada experimental não foi utilizada a adição de solução tampão. 

Para a operação em tempo real foram desenvolvidas VI em Labview® para cada 
controlador. Esse programa realiza a aquisição de dados, a interface gráfica do sistema e a 
interface matemático-computacional. Os controladores foram implementados em linguagem 
Matlab®. O diagrama de blocos genérico é mostrado na Figura 7.6. 

Utilizou-se 0 controle cascata, em que a saída do controlador MPC ou RMPC é o 
valor de setpoint para um controle PID. A adoção deste tipo de controle se deu como uma 
forma de corrigir a vazão devido ao efeito que a variação do nível no tanque de 
armazenamento exerce sobre a saída da bomba. Durante a fase de estudos foi também 
avaliada a aplicação do controle na forma de corrente para a bomba e não na forma de vazão, 
considerando a possibilidade de problemas nos transmissores de vazão, uma vez que se trata 
de um equipamento bastante sensível. Para este caso, não se utilizou o controle do tipo cascata 
e a ação de controle seria aplicada diretamente na bomba, mas este caso foi apenas estudado e 
não será tratado nesse trabalho. A implementação de fato utilizou-se da opção de controle 
cascata. 

As implementações em Eabview® e em Matlab® para os controladores MPC-RS, 
IHMPC, RMPC-MMR e EMI-RMPCR encontram-se no Apêndice A. 
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Figura 7.6 - Diagrama de blocos da malha fechada para controladores MFC e RMPC. 
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7.5.1 Resultados Experimentais 


Os controladores MPC-RS e IHMPC foram implementados na operação em tempo 
real na planta experimental. Os demais eontroladores (RMPC-MMR e LMI-RMPCR), por sua 
natureza, eomo eontroladores de horizonte infinito, tiveram os resultados de simulação 
eomparados eom os resultados de simulação e experimentais dos eontroladores de fato 
implementados na unidade, prineipalmente objetivando-se avaliá-los em uma mesma base 
para fins de eomparação. 

Na literatura, eneontram-se diversas téenieas de sintonização de eontroladores. 
Adotou-se nesse trabalho o roteiro apresentado por Henson e Seborg (1997) e as indieações de 
Montandon (2005), eujos desenvolvimentos foram feitos no eontrolador MPC-RS aplieado ao 
eontrole SISO da unidade experimental (eontrole de pH eom manipulação da vazão de a) e 
que, nesse easo, será adotado eomo um balizador de ações, de forma a se adquirir experiêneia 
na operação da unidade e na solução dos problemas que porventura aeonteeessem. Foi ainda 
realizado um ajuste fino de eampo durante a eorrida inieial em malha feehada eom eada 
eontrolador, realizando pequenos ajustes nos parâmetros de forma que o sistema apresentasse 
melhor desempenho experimental na malha feehada. A Tabela 7.1 mostra os parâmetros de 
eada eontrolador resultante da simulação. Não houve preoeupação em se obter o melhor 
eontrolador, e sim respostas que atendessem, de forma satisfatória, ao desempenho e á 
estabilidade. Desta forma os eontroladores foram inieialmente sintonizados eom matrizes 
identidades e horizonte de eontrole igual a 10 e foi feita a sintonia até a obtenção de uma 
resposta satisfatória. Os eontroladores de horizonte infinito mantiveram a mesma sintonia. 
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TABELA 7.1: PARÂMETROS DE SINTONIA DOS CONTROLADORES. 


Controlador 

Parâmetros 

IHMPC 

N^=10 

R=10I 

S=10I 

Aumax=0,25 mL/s 

RMPC-MMR 

N^=10 

R=10I 

S=10I 

S=0,01 

Aumax=0,25 mL/s 

LMI-RMPCR 

N„,=10 

R=10I 

S=10I 

5=0,01 

Aumax=0,25 mL/s 

MPC-RS 

K = 0,0051 

R = I ^ R = 1001 

S= 100001 

N„, = 2 

Np =10 

AUmax=0,5 mL/s 

Todos os 
eontroladores 

At=10s 

Umax=38 mL/s 

Umin=4 mL/s 


Foram realizadas variações de setpoint para se avaliar a eapaeidade dos eontroladores 
em lidar eom transições. Partindo-se do sistema eom pH igual a 7 e altura do reator em 18 em, 
foram realizadas variações na forma do vetor dado pela Eq. 7.47, sendo que eada degrau teve 
a duração de 2000 segundos, eonforme apresenta a Figura 7.7. 
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Figura 7.7 - Variação dos setpoints para as saídas em malha fechada - MPC-RS. 


Considerando o tempo das corridas e os problemas apresentados durante a operação 
da planta experimental em malha aberta e nas corridas para testes da malha fechada, optou-se 
pela divisão da corrida em duas partes. Em uma primeira etapa fixou-se o nível e foram feitos 
os degraus no pH. Em uma etapa posterior, fixou-se o pH e variou-se o nível. As corridas 
foram realizadas na seqüência, fazendo apenas a divisão para fins de gravação de dados e para 
evitar alterações significativas nas concentrações, devido, por exemplo, à absorção de dióxido 
de carbono (CO 2 ) pela solução de NaOH. 

Para todas as simulações de controladores foi adicionado um ruído aleatório de 
média zero na ordem de 10% nas saídas (h e pH). A Figura 7.8 mostra, como exemplo, a 
realização de uma sintonia fina durante a operação inicial e em tempo real executada com 
todos os controladores. Para essa situação 0 nível foi mantido constante. Foram alterados os 
valores de ajustes, que podem ser percebidos nos gráficos como variação da amplitude das 


ações de controle. 
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(a) 



Figura 7.8 - Corrida experimental para sintonia fina no MPC-RS - Destaque para variação da 
amplitude da ação de controle nos pontos modificados. 


A Figura 7.9-7.10 mostra a simulação para a malha fechada com o controlador 
MPC-RS. A resposta apresentada é satisfatória. As Figuras 7.11-7.12 mostram os resultados 
para a corrida experimental em malha fechada para o controlador MPC-RS. A Figura 7.11 b 
indica que durante a operação foi necessário fazer alterações na sintonia do MPC RS para que 
0 sistema se mantivesse no valor de referência desejado. Pela Figura 7.12 é possível notar o 
aumento no valor de X no ponto onde o sistema desvia do valor de referência. 
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Figura 7.9 - Simulação em malha fechada - MPC-RS: entradas e saídas. 
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Figura 7.10 - Simulação em malha fechada - MPC-RS: variável de folga. 
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Figura 7.11 - Corrida experimental com MPC-RS: saídas para transições no pH. 




















































































159 



Tempo(min) 


0.05 


0 



-0.05 ■ 


50 100 

Tempo (min) 


150 


Figura 7.12 - Corrida experimental do pH eom MPC-RS: entradas e variável de folga. 


A eorrida experimental para o pH foi eonsiderada satisfatória. O nível se manteve 
pratieamente durante toda a eorrida no entorno valor determinado, prineipalmente 
eonsiderando-se que a preeisão de leitura neste easo é da ordem de 0,1 em. A exeeção foi 
quando o pH sofreu um desvio do valor desejado no fim da eorrida. Contudo, é possível 
observar a atuação do sistema de eontrole no sentido de levar o sistema para o setpoint 
desejado. 

A Figura 7.13-7.14 apresenta a eorrida experimental do sistema para transições no 


nível eom o eontrolador MPC-RS. 
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(c) 


Figura 7.13 - Corrida experimental eom MPC-RS: saídas para transições no nível. 
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Figura 7.14 - Corrida experimental do nível eom MPC-RS: variável de folga. 

Pela Figura 7.13, pode-se observar que o controle de nível do sistema MIMO foi bem 
executado. Nota-se, nos instante iniciais da corrida, indicado pelas setas, que durante a parada 
da unidade na corrida anterior de transição o pH, ocorreu esvaziamento do sistema que foi 
corrigido para os valores de referência desejados. O controle do pH na região de 7, 
sabidamente ponto de grande dificuldade, durante a operação em tempo real se mostrou 
bastante satisfatório e eficiente. Contudo, depois de finalizada a corrida e gerados os gráficos, 
pode-se notar que o sistema apresentou um desvio do comportamento esperado, oscilando 
abaixo do valor desejado para o pH, podendo-se concluir que apesar de prova matemática a 
respeito, o sistema apresentou offset. O fato gerado para essa situação pode ser o observador 
de estado, uma vez que os estados utilizados não foram dados por filtro de Kalman e sim pela 
integração do modelo linearizado localmente. Outro fato que reforça essa idéia, é que o 
sistema foi para essa região na Figura 7.4, que trata dos problemas como o observador de 
estado. Para os controladores de horizonte infinito, adotou-se o uso de um observador de 
estado em malha aberta, ou seja, informando ao controlador o erro (e) entra a planta e o 
modelo, apresentado na Eq. 7.48, de forma a se minimizar a interferência da estimação de 
estados no comportamento do sistema (MACIEJOWSKI, 2002). 
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Nota-se também que no final da eorrida o valor do pH foi para uma região indesejada e não 
foi eorrigido, isto devido à não alteração da sintonia do eontrolador. A variável de folga 
refletiu esse desvio, eonforme indieado na Figura 7.14. 

Durante a operação, apesar de haver bastante euidado na elaboração da solução, de 
esgotamento/limpeza dos tanques e também no aeompanhamento das eorridas, foi observado 
que em eada eorrida havia variações exigindo adaptações ou resultando em respostas 
ligeiramente diferentes. Entretanto, é bom ressaltar que em um ambiente industrial, 
possivelmente, estes aspeetos abordados tenham seus efeitos poteneializados. 

O uso dessa estratégia (MPC-RS), que se difere dos eontroladores desenvolvidos, 
fez-se prineipalmente de forma a ganhar experiêneia na planta, eonsiderando-se que já houve 
outras pesquisas utilizando-se deste eontrolador aplieado ao eontrole de pH (este eontrole 
SISO já estava implementado na unidade) nessa unidade experimental e, assim, os problemas 
deseritos para esse easo podem ser mais bem tratados em um outro estudo. De qualquer 
forma, essa aplieação ofereeeu um eonheeimento sobre o sistema e uma visão erítiea da 
problemátiea que envolve uma eorrida em malha feehada, bem eomo a operação em tempo 
real, eontribuindo para a melhoria de desempenho do eontrolador de horizonte infinito. 

A Figura 7.15 apresenta a variação para as saídas eonsiderando as simulações e a 
operação em tempo real eom malha feehada para os eontroladores de horizonte infinito. Este 
tipo de entrada para o sistema não é estranho á prátiea industrial, prineipalmente eonsiderando 
as faeilidades da informátiea, que podem faeilmente tornar uma transição na forma de uma 
resposta de primeira ordem. Essa entrada de setpoint foi adotada devido ao fato de o sistema 
apresentar resposta oseilante e insatisfatória quando um degrau de grande amplitude era 
solieitado á planta. A outra opção, mais eomumente utilizada, é a implementação de um 
modelo ou de uma trajetória de referêneia ineorporada ao eontrolador. 
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Figura 7.15 - Variação dos setpoints para corrida em malha fechada - IHMPC. 


A Figura 7.16 mostra a simulação do sistema com o controlador IHMPC. Esta 
resposta é satisfatória. Para essa simulação foi considerado que o sistema encontrava-se em 
desvio planta/modelo, que foi caracterizado por um efeito tamponante equivalente a 0,01 
mL/s da vazão q 3 na planta, enquanto o modelo possuía a consideração de q 3 igual a 0,05 
mL/s. Estas considerações fazem com que o controlador não seja tão sensível como ele seria 
sem a consideração da solução tampão. Na corrida experimental não houve adição de tampão. 

A Figura 7.17 apresenta o tempo de cálculo do controlador IHMPC. Nota-se que a 
montagem do problema e resolução da otimização é realizada em um tempo bastante rápido, 
considerando-se que o tempo de amostragem é de 10 s. 
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Figura 7.16 - Simulação para malha fechada - IHMPC. 
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A Figura 7.18 apresenta a eorrida experimental para a transição do nível utilizando 
0 eontrolador IHMPC. Pereebe-se que o eontrolador manteve o nível em um patamar de 
melhor desempenho eomparado eom o eontrolador MPC-RS. O nível aparentemente apresenta 
um offset, eontudo esse valor eneontra-se dentro da ineerteza da medida. Ao final da eorrida 
foi neeessário ajustar a sintonia do eontrolador, eomo já havia sido feito eom o eontrolador 
MPC-RS. A Figura 7.19 apresenta a eorrida experimental para a transição do pH utilizando o 
eontrolador IHMPC. O nível foi mantido dentro do valor da ineerteza da medida. A eorrida é 
eonsiderada satisfatória. 

Considerando que os eontroladores robustos são téenieas eonstruídas utilizando-se 
do eontrolador IHMPC é natural fazer a eomparação em simulação e esperar que os resultados 
experimentais sejam eonsistentes eom as eonelusões desse estudo. 
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Figura 7.18 - Corrida experimental IHMPC: transição do nivel. 
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Figura 7.19 - Corrida experimental IHMPC: tradição do pH. 
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A Figura 7.20 mostra o erro entre os valores implementados pelo controlador e os 
valores obtidos com resultados experimentais para as vazões durante a corrida de variação de 
nível e de pH. Esse erro é dado pela dinâmica da estrutura mestre/escravo do sistema de 
controle e pelas diferenças ocasionadas pela medição. Esses valores são de interesse para a 
definição da incerteza na entrada, necessária para os controladores robustos. Nesse caso, pode- 
se observar que os valores encontram-se entre 2,32-4,83% de incerteza. De forma a se avaliar 
as respostas dos controladores com algum nível de incerteza, e não buscar casos que 
raramente aconteçam, define-se que o valor da incerteza a ser considerada pelos controladores 
robustos é de 2%. 
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Figura 7.20 - Erro nas entradas: (a) transições de nível; (b) trasições de pH. 
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A Figura 7.21 apresenta a simulação do sistema com o controlador RMPC-MMR. 
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Figura 7.21 - Simulação do sistema com o controlador RMPC-MMR. 
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A simulação indica uma ligeira perda de desempenho do eontrolador, mas, ainda assim, a 
resposta apresenta-se satisfatória. A Figura 7.22 apresenta o tempo de eáleulo do eontrolador. 
Conforme esperado, existe um maior esforço eomputaeional, eonsiderando-se o problema de 
otimização mín-máx. O tempo de cálculo é maior, mas dentro de valores que permitem a 
implementação no tempo de amostragem, que é de 10 segundos. As ineertezas ealeuladas pela 
otimização são mostradas na Figura 7.23. 
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Figura 7.22 - Tempo de eáleulo do controlador RMPC-MMR em simulação do sistema. 
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Figura 7.23 - Máxima Ineerteza ealeulada - RMPC-MMR: (a) Ineerteza em qi; (b) Ineerteza 
em q 2 ; (e) Ampliação. 
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A Figura 7.24 mostra as respostas do sistema eom o controlador LMI-RMPCR 
aplicado com sistema nominal. 
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Figura 7.24 - Resultados de simulação - LMI-RMPCR: caso nominal. 
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Os resultados dessa figura se diferem de um eontrolador RMPC pelo fato de que não é 
utilizado um politopo para a deserição da ineerteza e sim tão somente o modelo nominal. Os 
resultados apresentam uma degradação de eomportamento na região de pH igual a 7, mas, 
ainda assim, obteve-se um resultado de simulação satisfatório. A Figura 7.25 apresenta o 
tempo de eáleulo do eontrolador. Existe um aeréseimo do tempo de eáleulo ao longo da 
eorrida e pereebe-se que em um horizonte de tempo finito oeorrerá a ultrapassagem do tempo 
de amostragem inviabilizando a utilização do controlador. Considerando-se a implementação 
e todo 0 cuidado para se trabalhar com variáveis de forma otimizada, supõe-se que o 
acréscimo de tempo para esse controlador deve-se a limitações de memória, em função da 
dimensão das matrizes geradas. Para o cálculo desse controlador fez-se o uso do toolbox para 
LMI (Yalmip, LÔFBERG, 2004) e do LMlTools®. O uso das rotinas do toolbox YAEMIP 
toma a programação bastante intuitiva e direta, mas, em contrapartida, acarreta um gasto de 
tempo para estmturar o sistema na forma solucionável pelo LMlTools®. A Figura 126 mostra 
as respostas para o controlador considerando a incerteza na forma de politopo. Do ponto de 
vista de simulação a resposta é satisfatória. A Figura 7.27 mostra o tempo de cálculo do 
controlador. Percebe-se que o tempo ultrapassa o tempo de amostragem inviabilizando a 
aplicação nessa condição. 
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Figura 7.25 - Resultados de simulação - EMI-RMPCR: caso nominal. 
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Figura 7.26 - Resultados de simulação - LMI-RMPCR: easo eom politopo. 






























































174 



Figura 7.27 - Resultados de simulação - LMI-RMPCR: easo eom politopo. 


7.6 Problemas Apresentados Durante a Operação em Tempo Real 


Durante a operação em tempo real, diversos problemas oeorreram, invalidando a 
eorrida em eurso. Foram situações que podem ser eonsideradas eomuns em se tratando de 
atividades industriais, e dentre as quais podem ser eitadas: ruídos de natureza elétriea, ruídos 
de proeesso e vazamento, dentre outros. A exeeução do experimento desse trabalho exigiu 
mais do que eonheeimento teórico. Em alguns casos é necessária uma criteriosa análise do 
fato para se obter uma solução satisfatória. Em outros casos, a solução é obtida por tentativa e 
erro. De qualquer forma, pode-se considerar que a experiência resultou em uma oportunidade 
para o uso de ferramentas teóricas bem como de senso de engenharia, tão necessários para a 
solução dos problemas práticos e para a formação do pesquisador. 

A Figura 7.28 mostra a região de transição para a agitação do reator. Não é 
recomendável trabalhar na faixa de 4 a 6 centímetros de nível, principalmente considerando-se 
a oscilação, conforme mostra a Figura 7.28 b. Devido ao acoplamento do sistema, é esperado 
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que haja interferência no pH, e dependendo da região de trabalho poderá ocorrer instabilidade 
no sistema. Outra fonte de interferência no processo derivada da agitação é a vibração causada 
pelo motor do agitador, mostrado na Figura 7.29. 



Figura 7.29 - Motor do agitador: fonte de ruídos no sistema. 

Durante a operação experimental, tanto a vazão de ácido, quanto a vazão de base 
apresentaram vazamentos de solução em diversos pontos. Alguns desses pontos encontram-se 
indicados nas Figuras 7.30-7.32. Contudo, os maiores problemas foram na vazão de base. A 
necessidade de utilização de vários filtros de partículas devido à baixa qualidade da água, um 
deles em especial, mostrado na Figura 7.32, tornou-se fonte de entrada de bolhas no sistema. 

























176 


cujos efeitos são apresentados na Figura 7.33. A Figura 7.34 mostra o ponto onde se pode 
aeompanhar a passagem de bolhas. 




Figura 7.30 - Pontos de vazamentos de solução no sistema: bombas. 



Figura 7.32 - Filtro de partíeulas: bolhas no sistema. 
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Figura 7.33 - Ruído causado pela passagem de bolhas no sistema. 



Figura 7.34 - Visualização da passagem de bolhas no sistema. 


A Figura 7.35 apresenta os ruídos de ordem elétriea na vazão de base. Foram usados 
filtros digitais de sinais, eontudo, quando o ruído é de grande amplitude, aeonteee uma 
retenção na medida em que leva a efeitos indedesejados no sistema. Na Figura 7.36 são 


apresentados os ruídos de proeesso. 
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Figura 7.35 - Ruídos de ordem elétriea na vazão de base. 


i 12.1 

0) 

(fí 

iS 12 

0) 

■D 

11.9 

N 
(O 
> 

11.8 

130 135 140 145 150 155 160 

Tempo (min) 


Figura 7.36 - Ruídos de proeesso na vazão de base. 
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A Figura 7.37 apresenta um ruído transitório, aleatório e indesejável que oeorreu 
durante uma eorrida. A Figura 7.38 mostra uma eorrida em que a vazão de base sofreu 
diversos problemas. Para esse easo, baseado em experiêneia própria, infere-se que o problema 
oeorra devido a algum aqueeimento ou anormalidade na instrumentação, eausando essa 
interferêneia, uma vez que o problema se manifestava após eonsiderável tempo de eorrida. 

A Figura 7.39 mostra os ruídos na altiua do reator. Cabe ressaltar que a variável real 
manteve seu valor inalterado durante as medições deseritas. 
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Figura 7.38 - Ruídos na vazão de base. 
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Figura 7.39 - Ruído na medição da altura do reator. 
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A Figura 7.40 mostra uma corrida em malha aberta eompleta, mostrando entradas e 
saídas onde são indieados efeitos de ruídos. 
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Figura 7.40 - Malha aberta: resultados teóricos e experimentais - ruídos. 
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As Figuras 7.41-7.42 apresentam os efeitos dos ruídos para eorridas eom o MPC-RS. 
A Figura 7.43 mostra o efeito de inieialização ineorreta para o eontrolador IHMPC. 



Figura 7.41 - MPC-RS: resultados experimentais e ruídos. 




Figura 7.42 - Malha fechada MPC-RS: ruídos. 
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Figura 7.43 - Malha fechada IHMPC: inicialização incorreta. 

As Figuras 7.44-7.45 mostram a montagem elétrica do sistema cujos pontos de 
conexão podem causar a captação e introdução de ruídos no sistema. 



Figura 7.44 - Sistema de aquisição de dados. 
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(a) (b) 

Figura 7.45 - Circuitos elétricos da instrumentação da planta experimental. 


A rede elétriea que alimenta o laboratório apresenta visíveis oscilações, indicando possíveis 
sobreeargas no cireuito de alimentação. A instalação elétriea do laboratório também neeessita 
de modifieações para eliminação de interferêneias, prineipalmente quando se trabalha eom 
equipamentos sensíveis, eomo os de instrumentação e aquisição de dados. Nos dias em que as 
corridas foram realizadas, foi necessário que outras atividades que pudessem provoear ruídos 
fossem paralisadas ou adiadas. A unidade experimental neeessita ser desmontada e remontada 
utilizando-se de fontes estabilizadas e conexões soldadas e com rotas otimizadas e isoladas de 
forma a se eliminar os pontos de ruído, inelusive fazendo-se uso de blindagem de cabos. A 
alocação da instrumentação deverá ser feita de forma que vazamento ou transbordamento no 
reator não eheguem aos instrumentos. Existe a neeessidade também de uma modifieação nos 
eircuitos hidráulieos para eliminação de pontos de vazamentos e bolhas e também a 
substituição ou modificação dos filtros. A questão fmaneeira diante das neeessidades de 
aquisição e modificação do sistema mostra-se crucial. Uma simples unidade como essa 
demanda um razoável volume de reeursos, o que inviabiliza, em alguns casos, compras 
necessárias, forçando a adaptações que permitam a eontinuidade da pesquisa sem que a 
melhor solução seja implementada. Para algumas implementações de rotinas, eomo filtros, 
deteeção de mudanças bruseas e outras nesse sentido, podem ajudar a minimizar os 
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problemas. Em outras situações seria neeessária a troea de fontes ou a instalação de baterias 
eomo fonte de eorrente eontínua. Outras situações poderiam ser resolvidas eom a duplieidade 
de sensores e atuadores. As sugestões de melhoria são para que se eoneentre o esforço em 
trabalhos futuros na rotina de eontrole. Contudo, as difieuldades se mostraram salutares para o 
exereíeio pleno da engenharia e da pesquisa. 



7.7 Seção Comparativa 
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Situação: 0 - respostas atendem aos requisitos; - Não avaliado; ^ - Necessidade de maior estudo para solução dos problemas apresentados. 






186 


7.8 Considerações Finais 


Os controladores RMPC-MMR e LMI-RMPCR foram aplicados a um problema real 
representativo dos sistemas de controle. O modelo fenomenológico foi obtido e avaliado 
experimentalmente, sendo considerado adequado para os propósitos desejados nesse trabalho. 
Os controladores MPC-RS e IHMPC mostraram respostas satisfatórias na aplicação em malha 
fechada, com modelo perfeito. As respostas experimentais são avaliadas pelas simulações de 
problema com desvio planta modelo. Os controladores MPC-RS e IHMPC foram 
implementados na operação em tempo real e os resultados das respostas com esses 
controladores apresentaram-se adequados e satisfatórios. O controlador RMPC-MMR foi 
avaliado por simulação, apresentando uma degradação em tomo do ponto de pH igual a 7, 
mas ainda sim aceitável. A resposta de simulação para o controlador LMI-RMPCR é aceitável 
do ponto de vista de simulação, mas é inviável para aplicação real, considerando-se o elevado 
esforço computacional, que supera inclusive o tempo de amostragem. A sua aplicação requer 
novo estudo de forma a se trabalhar com tempos de amostragem maiores. Foram também 
mostrados e discutidos diversos problemas ocorridos na operação da unidade experimental. A 
seguir, no Capítulo 8, são apresentadas as conclusões desse trabalho bem como as sugestões 


para futuros trabalhos. 



CAPITULO 8 
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8 CONCLUSÕES 


Nesse trabalho foram propostos dois algoritmos de eontrole robusto: 

RMPC-MMR: RMPC formulado eom otimização min-máx, eom ineerteza invariante no 

tempo na entrada do sistema e eom restrição de robustez; 

LMI-RMPCR: RMPC formulado eom otimização na forma de LMI utilizando ineerteza na 

forma politópiea e restrição de robustez. 

Os resultados apresentados nos Capítulo 4, 6 e 7 permitem listar as seguintes 

eonelusões: 

• As respostas de simulação das saídas dos sistemas submetidos aos eontroladores se 
mostraram adequadas e satisfatórias; 

• As respostas dos eontroladores desenvolvidos demonstraram desempenhos semelhantes 
aos algoritmos elássieos de MPC e RMPC. 

• Os algoritmos de eontrole RMPC desenvolvidos eonseguiram lidar bem eom os problemas 
freqüentemente eneontrados nos sistemas de eontrole apresentados pelos sistemas 
estudados, tais eomo: desvio/planta modelo, ineertezas aleatórias, diserepâneias entre fases, 
tempo morto, ganho, eonstante de tempo, modelo linear/planta não linear, modelo estável / 
planta instável, mal eondieionamento; 

• A planta experimental foi modelada e teve os parâmetros neeessários identifieados através 
de eorrida em malha aberta. Os resultados de simulação foram satisfatórios e indiearam que 
0 modelo é representativo do proeesso real e adequado para os propósitos de eontrole 


pretendidos; 
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• As respostas experimentais do eontrolador MPC-RS foram satisfatórias. 

• As respostas teórieas e experimentais do eontrolador IHMPC foram satisfatórias. O ajuste 
do eontrolador para operação em tempo real se mostrou adequado. 

• Baseado nas respostas do sistema experimental submetido ao eontrole IHMPC, foi 
definido que uma ineerteza na entrada representativa desse sistema é de 2% do valor de 
entrada. 

• As respostas da aplieação dos eontroladores LMI-RMPC e RMPC-MMR ao sistema 
experimental mostraram-se satisfatórias do ponto de vista de simulação. 

• O RMPC-MMR eom ineerteza na entrada limitada a 2% foi ealeulado dentro do tempo de 
amostragem e eom desempenho satisfatório; 

• O LMI-RMPCR eom modelo nominal (sem ineerteza) já indieou que haveria problemas 
para o eáleulo da ação de eontrole dentro do tempo de amostragem; 

• O LMI-RMPCR eom ineerteza na forma de politopo eobrindo uma região de ineerteza de 
2% na matriz de entrada não eonseguiu ser ealeulado dentro do tempo de amostragem, o 
que torna a aplieação real para esse easo inadequada. Faz-se neeessário o trabalho em outro 
tempo de amostragem, que permita o eáleulo do eontrolador e que não eomprometa o 
desempenho do sistema. É possível esperar, diante da evolução do proeessamento 
eomputaeional, que em um futuro próximo esta limitação seja eliminada ou bastante 
reduzida; 

• A interferêneia de questões relaeionadas eom vazamentos, ruídos de proeesso, ruídos de 
natureza elétriea e outros problemas assoeiados eom a aplieação em tempo real, trouxeram 
difieuldades adieionais que exigiram habilidade manual e eonheeimento do proeesso. 
Sugere-se a reeonstrução da unidade experimental de forma a se eliminar fontes de ruídos e 
vazamentos, para que futuras pesquisas eoneentrem esforços no desenvolvimento dos 


eontroladores. 
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8.1 Sugestões para Trabalhos Futuros 


• Utilização dos controladores robustos desenvolvidos em modelos empíricos, como por 
exemplo, modelos neuronais; 

• Aplieação de controle robusto no modelo experimental eonsiderando o eontrole da 
corrente das bombas e não o controle de vazão; 

• Reeonstrução da unidade buseando eliminar os problemas deteetados nesse trabalho; 

• Estudo experimental de controle robusto assoeiado a controle tolerante a falhas; 

• Análise e melhoria do esforço eomputaeional no sentido de reduzir o tempo de 
proeessamento dos controladores robustos. 
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APÊNDICE A 


A.l. Programas para Operação em Tempo Real 


Para a operação da planta foram implementadas algumas rotinas no software 
Labview®. Através destas rotinas é possível: 

• Iniciar/parar/eneerrar a operação da planta experimental; 

• Fazer a exeitação das vazões de entrada; 

• Adieionar/remover filtro digital de sinais; 

• Seleeionar modo de Operação Manual/Automátieo, ajustando setpoints desejados. 

Nesse módulo é ainda possível visualizar: a eorrente das bombas (modo manual ou modo 
automátieo); a tensão dos transmissores envolvidos no eontrole; as entradas e as saídas 
(valores adquiridos e valores desejados); tempo de aquisição; tempo de iteração e valores de 
ajustes dos eontroladores PI para as vazões de áeido e de base. A Figura A.l mostra o painel 
frontal da VI para a exeeução da eorrida em malha aberta e a Figura A.2 apresenta o diagrama 
de bloeos referente a essa implementação. 

As Figuras A.3-A.4 mostram o painel frontal e o diagrama de blocos, 

respectivamente, para o controlador PID das bombas. 

As Figuras A.5-A.6 mostram o painel frontal e o diagrama de blocos, 

respectivamente, para filtro digital de sinais implementado. Para maiores informações sobre 
os filtros implementados nesse trabalho ver Gonçalves (2007). 
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Figura A. 1 - Painel Frontal do VI Malha Aberta. 


























































































































Figura A.2 - Diagrama de blocos da VI Malha Aberta. 
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Input PID parametersl 


proc. var. 

P 

Sídíõõ I 

Slõiõõ 1 

setpoint 

I 

SPõ 1 

Slõõõ 1 

upper limit 

Slõiõõ 1 

D 

lower limit 

3IÕÍÕÕ 1 

Slíõõ 1 



previous controller outputi 

glO.OO I 


controller outputl 

|D,00 I 


Figura A.3 - Painel frontal da sub-VI do Controlador PID. 



Figura A.4 - Diagrama de bloeos da sub-VI Controlador PID. 



Figura A.5 - Painel Frontal da sub-VI Filtro Digital. 
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Figura A.6 - Diagrama de blocos da sub-VI Filtro Digital. 


As Equações A1-A5 mostram as curvas de calibração da instrumentação. Para 
maiores detalhes sobre o procedimento ver Gonçalves (2007). Essas curvas, juntamente com 
os coeficientes de regressão, são apresentadas nas Figuras A7.-A.11. 


qj =5,2788Vo-9,4776 

(A.1) 

q^ =5,0997V,-9,5541 

(A.2) 

q3 =0,2090V2 -0,3845 

(A.3) 

pH = l,8982V3 -3,5583 

(A.4) 

h =5,6347V4-9,0838 

(A.5) 
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Figura A.7 - Curva de calibração e regressão linear dos dados - Vazão de Ácido. 



Tensão Transmissor - Bomba de Base(V) 


Figura A.8 - Curva de calibração e regressão linear dos dados - Vazão de Base. 
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Tensão do Transmissor - Bomba de Solução Tampão (V) 


a A.9 - Curva de calibração e regressão linear dos dados - Vazão de solução tampão. 



Figura A. 10 - Curva de calibração e regressão linear dos dados - Altura (h). 
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Tensão do Transmissor de pH (V) 


Figura A. 11 - Curva de calibração e regressão linear dos dados - pFI. 

As Figuras A12-A13 apresentam o painel de eontrole e o diagrama de bloeos, 
respeetivamente, implementados para o controle em malha fechada com MPC-RS. 

As Figuras A14-A15 apresentam o painel de controle e o diagrama de blocos, 
respectivamente, implementados para o controle em malha fechada com IHMPC. Na Figura 
A.16 é mostrada a condição “7>we” para aplicação do controlador, que nesse caso entra em 
funcionamento controlando o sistema. A Figura A. 17 mostra a condição ‘‘"False” e, nesse caso, 
apenas as medidas de interesse são mantidas para a inicialização do controlador quando posto 
em funcionamento. A Figura A. 18 indica a condição “True” para a excitação dos valores de 
referência seguindo a transição desejada, conforme mostrada nas Figuras 7.7 e 7.15. 

Para os controladores RMPC-MMR e LMI-RMPCR a estrutura permanece a mesma 
realizada para o IHMPC, sendo acrescentado apenas o termo referente à entrada de dados 
relacionados com a incerteza a ser aplicada, considerando-se os controladores robustos, como 


pode ser observado na Figiua A. 19. 



211 


«HJ 




m 


m 




niTO rm Ita^ic 

liwmiTkl^l Pni^rp) ^nhi^nii'! 

Exõ |to 


(fr Snnrwfl^ 


UU] ■AKUnl 


rcRvtp 



í ÍMÍ ■ " 
ÍU ' 


(? 


# 


u-iniii 


KRh'* 


■^1 




w 


-kttpMc' ifMhl >MíM: 


íjM» •;»» iis,ii 

P L D 

i!!£É;Çà^Ê 5 


Í3 


loíin 'tt* KbbIH 

}üa~lBji I 


Sk-w lir™ igun 





«»i:i>í 




ur 

í 

*1 

í 


■r 


0 » 


v^w" H 


a,N 


PittinriÉ m 

■ 

0 » 

Srl 

[^ |iir¥ 





—iwraççrPBTTraííB- 

iewc 


^ id»* 

P S “ 

CsrP -» Im (P1 

|! iílB 


k| c 



fWrfi 

vwt 

Pll et 


OÍ»Kh 


Figura A. 12 - Painel Frontal da VI implementado para MPC-RS. 
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Figura A. 13 - Diagrama de blocos da VI implementado para MPC-RS. 
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Figura A. 14 - Painel Frontal da VI implementado para IHMPC. 
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Figura A.15 - Diagrama de blocos da VI implementado para IHMPC. 
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iMATLAe Scriptf 

global Ki^ Kal Ka2 Q2 At Wa Wb cv p u 
options = optimset('DisplayVoff); 

...System data 
At =168,3870; 
cv =20,4477; 
p = 0,0523; 

Q2 =0,001; 

HN03cl =0,0035011014778; 

NâHC03c3 = 0,0; 

NdOHc3 = 0,0035011014778; 

N3HC03c2 = 0,035011014778; 

Kw = le-14; 

Kal =4,47e-7; 

Ka2 =5,62e-Il; 

,Wa(l) =HN03cl; 

'b(l) =0; 

Wa(2) =-NaHC03c2; 

Wb(2) =NaHC03c2; 

Wa(3) =-Na0Hc3-NaHC03c3; 

Wb(3) =NaHC03c3; 

., .number of States fo the system 
Neq =3; 

^«State L^date 
y =[yl.k;y2_k]; 

ym «[yml_l;ym2J]; 

7«,,Initiali2ation 
u = [ul k;u2 k]; 

X = [yl.k;x2_lgx3_k]; 

...Setpoints vector 
ystp = [stpl;stp2]; 

Q * diâg(Q); 

R = diagíR); 

S =áag(S); 

%.. Cálculo de DeltaU para Umin < u < Umax 

Umax = Umax’; 

umin = UmIn'; 

deltumax = deltUmax'; 

deltumin = -deltumax; 

,. .linearízabon around nominal point 
h = yl k) 

WA = x2_k; 

WB = x3 k; 

Ql = ul k; 

Q3 = u2_k; 

A = [-(cv*p*h^(p-l))/At 0 0;-((Ql’^(Wa(l)-WA)+Q2*(Wa(2)-WA)+ 


Figura A. 16 - Destaque para bloco do controlador - IHMPC, condição ‘Trae”. 
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Figura A. 17 - Destaque para bloeo do eontrolador - IHMPC, eondição “False'\ 










































































Figura A.18 - Destaque para a eondição '"False’" para excitação de setpoint. 
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Figura A. 19 - Painel Frontal da VI implementado para os eontroladores RMPC. 
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